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PREFACE 

This report describes part of a comprehensive and continuing program of re- 
search concerned with advancing the state-of-the-art in remote sensing of the en- 
vironment from aircraft and satellites. The research is being carried out for the 
National Aeronautics and Space Administration (NASA), Lyndon B. Johnson Space 
Center, Houston, Texas, by the Environmental Research Institute of Michigan 
(ERIM), formerly the Willow Run Laboratories of The University of Michigan. The 
basic objective of this multi-disciplinary program is to develop remote sensing as 
a practical tool to provide the planner and decision-maker with extensive informa- 
tion quickly and economically. 

Timely information obtained by remote sensing can be important to such people 
as the farmer, the city planner, the conservationist, and others concerned with prob- 
lems such as crop yield and disease, urban land studies and development, water 
pollution, and forest management. The scope of our program includes (1) extending 
the understanding of basic processes , (2) discovering new applications, develop- 
ing advanced remote sensing systems, and improving automatic data processing to 
extract information in a useful form, and (3) assisting in data collection, process- 
ing, and analysis as well as in data verification via ground truth. 

The research described herein was performed under NASA Contract NAS 9- 
9784, Task II, and covers the period from 1 February 1973 through 31 October 
1973. Dr. Andrew Potter was Technical Monitor. The program was directed by 
R. R. Legault, Vice President of ERIM, J. D. Erickson, Principal Investigator 
and Head of the Information Systems and Analysis Department, and by 
R. F. Nalepka, Head of the Multi spectral Analysis Section. The ERIM number 
for this report is 190100-23-T. 

The authors wish to acknowledge their appreciation to Richard F. Nalepka for 
the many valuable suggestions he offered during this study and to Thomas J. 
Marshall for his contributions in the preparation of this report. 
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SYSTEMATIC MONITORING AND EVALUATION OF M7 
SCANNER PERFORMANCE AND DATA QUALITY 

1 

INTRODUCTION AND SUMMARY 


1.1 INTRODUCTION 

The purpose of this study is to establish and monitor the performance of the multispectral 
scanner operated by ERIM as an aid in the processing of data gathered as part of NASA's Earth 
Resources Program. Our goal is not simply to perform radiance calibrations and provide quan- 
titative measurements by which users can extract radiance information. Rather, it is to inves- 
tigate the temporal variations and overall uncertainties in scanner performance characteristics 
in order to meet data processing and interpretation requirements. Those areas potentially most 
helpful in improving data processing capabilities are emphasized. Efforts are also being made 
to improve radiance calibration techniques (including more efficient data reduction). 

1.2 SUMMARY 

Many users of remote sensing data gathered by the M7 scanner system are interested in 
extending spectral classification and mapping capabilities from one flightline to another, even 
when such lines are flown days apart. Such a capability obviously offers tremendous savings 
in manpower (including ground truthing teams) as well as in data processing time. In most 
cases, however, attempts to recognize objects in one flightline using signatures and decision 
rules generated from another line (even with flightlines only several minutes apart) have 
been characterized by significant decreases in the percentages of objects correctly classified. 

From the viewpoint of the data processor, it would be very helpful to know the magnitude 
of the scanner-associated signal variations one might expect during the course of a three- or 
four-hour mission. Information on longer-term signal variations covering three or four months 
would also be helpful. By knowing all the kinds of signal variations and their characteristics to 
be encountered over long periods of time— including scanner, atmospheric, and seasonal changes 
— the processor has a better chance of extending object classification capabilities from one 
flightline to another. 

Toward answering the question of short- and long-term scanner signal variations, one of 
the approaches we used was to analyze the so-called "standardization" tests conducted period- 
ically in the course of normal data collection over the past year and one-half. These tests were 
devised originally to indicate changes in system response. Forty tests, covering the periods 
5 May 1972 to 25 January 1973 and 22 June to 12 August 1973, were analyzed with a program 
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called CALIB., which was written expressly for the purpose of calculating a set of statistics con- 
cerning the solar reflective reference sources of the M7 scanner. General characteristics of 
the M7 scanner are documented in [1] . Some of these statistics, plotted by the computer as 
graphs of mean signal voltage versus test date, showed that under normal circumstances, stan- 
dard deviations of the mean lamp signal amounting to not more than 10% could be expected over 
several months and that in most cases the figure was less than 5%. A few instances of much 
larger changes were noted but these were under unusual circumstances, as we will discuss 
later. 

Since these standardization tests are planned to serve as continuing indicators of scanner 
performance, a system for rapid, computerized data reduction was designed. Once implement- 
ed, this system will enable one-day retrieval of desired data on the performance characteristics 
as determined in the standardization tests. 

Also completed during this study were three other programs designed to help determine the 
impact of scanner signal variation and calibration on data processing. One of these programs, 
THERML., provides a set of statistics on thermal reference sources as a function of time — much 
as the program CALIB . does in the solar reflective. NEATs* for the hot, cold, and ambient 
reference plates were calculated and compared as a function of time to qualitatively assess 
thermal data performance. In addition, a new technique for evaluating the accuracy of tempera- 
ture calibration was perfected using THERML . output. Also, we made a quantitative analysis of 
the frequency content of the thermal channel by using the output of another program called 
DGNSTC . The results of this study indicated that a line-by-line dynamic temperature calibra- 
tion was necessary to minimize errors resulting from noise. DGNSTC- was also used to deter- 
mine the darkest region in the scanner’s field-of-view and thus permit uniform baseline correc- 
tions for all processors of M7 data. 

The last program, called L AMR AD., was written to automatically calculate M7 scanner 
radiance calibration constants for each spectral channel. New methods employed in this sub- 
routine greatly improved the signal-to-noise ratio of the scanner signals as compared to old 
techniques. NEALs* (radiance) were also calculated by LAMRAD. 


*Noise Equivalent Temperature Difference, defined later in Section 5.3. 
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2 

M7 SCANNER NOISE CHARACTERISTICS (SOLAR REFLECTIVE REGION) 

Noise figures for the M7 multispectral scanner (and, for that matter, any scanner) have much 
more value than just indicating system performance. It is true that noise level monitoring does 
help indicate overall system performance; but to data processing personnel, noise figures 
especially the frequency content of the noise — can provide enlightening information as to which 
noise in the data is of sensor origin, and which results from other physical phenomena. Such 
significant discrimination must be achieved if optimal use is to be made of remotely sensed 
data. It is questionable at this stage in technique development, especially in data collected 
over agricultural areas, as to whether some of the signal variation problems encountered stem 
primarily from the scanner or from changes in the scene. In all probability, it is some com- 
bination of the two. The remainder of this section is aimed at classifying the different types of 
noise to be expected from the scanner. The noise is evaluated over different time intervals to 
help the data processor analyze specific problem areas, such as apparent effects during one 
scanline (high frequency noise) and during signature extensions (low frequency noise). 

2.1 LOW FREQUENCY (LONG-TERM SIGNAL STABILITY) 

Low frequency noise, as analyzed in this section, means signal variations in the M7 over 
a period on the order of days and weeks. During this length of time (days to even months), 
some changes in the optics, reference lamp spectral characteristics, and detector response are 
to be expected. Our main concern in observing such signal variations was to determine the 
frequency and amplitude of the changes as an indicator of system performance. Of secondary 
interest was the possibility that signature extension of data from one flight to the next several 
weeks later may be attempted in the near future. 

2.1.1 PROGRAM CALIB. 

The purpose of program CALIB. is to permit quick data handling and analysis of the stan- 
dardization test performed routinely on most data-collection missions (see next section). 

This program calculates the following set of statistics: 

(a) mean and standard deviation of the dark level (in absolute volts) 

(b) mean and standard deviation of the lamp pulse level (in volts relative to dark level) 

(c) mean and standard deviation of the middle reselm level in the lamp pulse (in volts 
relative to dark level) 

(d) mean and standard deviation of the position of the lamp pulse (average of 50% points 
in dark-level -corrected lamp pulses) 
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(e) mean and standard deviation of the integral of the solar reference pulse (in volt-reselms 
relative to the dark level). 

All voltages will have been initially adjusted to correct for any gain changes. A table of 
all five sets of statistics is printed out, as well as a plot of the mean voltage and standard de- 
viation of (a), (b), and (c) above versus standardization run number. Since these standardiza- 
tion runs occur days or weeks apart, the plots effectively indicate the scanner response versus 
time relationship. One of the desirable features of this plot routine is the capability of storing 
previously calculated standardization data on a special tape data base for updating as new data 
are gathered. Hence, whenever CALIB. is run on a standardization data set (or sets), output 
from that last set (or sets) is plotted along with all previously acquired data. An updated history of 
scanner response over any desired length of time is available for immediate analysis — subject, 
of course, to keeping the standardization tests current. A procedure for achieving rapid data 
reduction and thus fast turn-around time for these tests is discussed later. 

Other innovative features of the CALIB. program include new methods to calculate lamp and 
sky reference values* A common method in the past was to obtain a histogram of the lamp 
pulse over all the scanlines in the data set. The reselm corresponding to the peak lamp voltage 
from the histogram was then chosen as the point on the pulse to use on all channels. Unfortun- 
ately, when skewing errors are present, a single reselm cannot describe the position of the 
pulse peak on all channels simultaneously. Also, if tape jitter is present, the position of the 
lamp peak will shift from line to line by nearly the amount of the jitter. For jitter of more than 
1 or 2 reselms, a single reselm will not describe exactly the peak lamp voltage. One of the 
techniques being used in this program to determine the lamp pulse voltage level is given below 
(also see Fig. 1): 

(a) first, the maximum (MAX) magnitude of the lamp signal is found (this may, unfortunate- 
ly, include noise spikes) 

(b) next, the reselms closest to MAX/2 are found (points X and Y defining the 50% points) 

X + Y 

(c) then, the middle reselm — ^ — is calculated 

(d) bounds on the area of the lamp plateau used for the level determination are then set on 
either side of ^ ^ ^ any °dd num ber of reselms (called TOP) can be used; the plateau 
bounds are ±(TOP-l)/2; TOP must be an odd number to permit symmetric bounds. 

For purposes of our test, TOP was set at 5 reselms (although any reasonable width could have 
been used). The default case was also set at 5 reselms. 


*A typical scanline is depicted in Section 3. 
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FIGURE 1. LAMP-PULSE VOLTAGE VERSUS RESELM POSITION 
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An advantage of the technique outlined above is that it is independent of any lamp jitter 
present in the data. Also, it can reveal the amount of jitter present by calculating the standard 
deviation of the middle lamp reselm (X + Y)/2 position. There are several disadvantages, how- 
ever. It is subject to possible error resulting from large noise spikes in the flat area of the 
pulse. Calculations of the 50% points, if based on a noise spike maximum, would be wrong. 
Moreover, only a few data points per line (typically five) are available for statistical calcula- 
tions. Hence the standard deviation is larger than need be. In addition, both of these disad- 
vantages are compounded with the present scheme for determining lamp voltage levels. Such 
problems may be overcome, however, by means of a new technique using integration; it is out- 
lined below for the sky reference, and will be analyzed in more depth in the section on program 
LAMRAD. The advantage of integrating is that many more points per line are available for 
statistical averaging. This greatly improves signal -to -noise, and noise spikes do not cause the 
problems they do in the present system. 

The integration technique mentioned above was initiated primarily because of the geometri- 
cal constraints imposed on the present sky reference. The pulse shape from the scanner is tri- 
angular instead of nearly square — as the lamp pulse is. Hence, the maximum value of the 
signal is defined by at most one or two points; this makes it almost impossible to obtain an accu- 
rate measurement. With integration, the value of each and every dark-level -corrected reselm 
making up the pulse is summed. The resulting sum is a number proportional to the maximum 
or peak value (since the width of the pulse is fixed by the geometry of the sensor). Actually, a 
number close to the real maximum value can be obtained simply by dividing the integrated value 
by the half -width (or some other representative width). 

In practice, the region integrated is much larger than the pulse itself. This allows for pulse 
movements by reason of tape or scanner jitter. Since all the values are dark -level corrected, 
any reselms falling outside the pulse but inside the integration region will contribute nothing to 
the integral. Note that tape and scanner jitter do affect the integral approach since tape speed 
variations cause some fluctuations in pulse width. However, the magnitude of error from jitter 
is considerably smaller than those errors associated with current methods of determining lamp 
pulse voltages. A complete, detailed description of the program CALIB (including a program 
listing) may be found in Appendix A. 

2.1.2 STANDARDIZATION RUNS 

During the time of this study, 40 standardization data sets, covering three distinct periods 
from 5 May 1972 to 25 January 1973, and a fourth from 22 June to 12 August 1973 were digitized 
and the data analyzed by the program CALIB. No data were digitized from 26 January through 
21 June because of limited time and funds. 
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These standardization data sets were digitized so we would have a basis for comparing 
scanner response from one mission to the next. To achieve this, we adopted a standardized 
procedure for recording reference lamp voltage approximately once every mission. The pro- 
cedure uses a fixed lamp current setting and fixed amplifier gains. Assuming no amplifier 
gain changes and no spectral or ir radiance changes in the reference lamp (reasonable assump- 
tions over several months) , a shift in lamp voltage must indicate a change in overall system 
response. These changes, then, are basically the quantities being analyzed by the program 
CAUB. 

As it turned out, fixed amplifier gains were not maintained over the time periods of interest. 
Changes in some system components required that different gains be used. For this reason, 
the gains were compensated for by CALIB. Shown in Fig. 2 are plots of lamp response versus 
time for each spectral channel available (except thermal) covering only the first three time 
periods mentioned above. The computer plots shown in Figs. 3 (a - n) give a more detailed look 
at response for each channel. The overall time frame spanned in these detailed plots includes 
these three time periods plus a fourth which covers the most recent months of scanner operation. 
It should be noted that these response values are plotted against mission number. The date of 
each mission is indicated in Table 1. As mentioned, the data covers four main periods of 
scanner operation. 

Thefirst period covers the months of May and June 1972 which were bracketed between major 
scanner maintenance during the winter and performance tests the following July. As canbe seen 
from the plots, absolute system response does vary somewhat during this period, although the 
standard deviation of the mean response goes from a maximum of 10% in channel CIO (0.46- 
0.49 jum) to only 3% in C6 (0.55-0.60 iim). It should be noted, however, that even though the 
overall response may have varied as much as 10% over a two-month interval, the variations 
observed over one -day intervals (by two missions flown on the same day or, at most, on the 
next day) showed variations typically on the order of 2% or less. Hence, evidence from past 
data indicates a high degree of short-term (and here we mean one day or less) stability. 

Note from Fig. 2 that some degree of total spectral response shifts is indicated by the 
graph. For instance, all channels except CIO show a decrease in response between 4 May 1972 
and 5 May 1972; likewise, a decrease between 11 May and 12 May 1972 is noted for most of 
the channels. It is obvious, however, that not all channels change In the same direction all the 
time. A major exception is channel C8, which had a PM tube replaced midway during this 
period— the change there is obvious. 

The second period covered by this data set occurs during the months of July and August, 

1973, at a time when the scanner was undergoing a series of bench performance checks. We 
discovered in these checks that the lamp housing containing the secondary radiance standard 
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TABLE 1. MISSION NUMBERS AND DATES 


Mission 

Number 

Date 

Mission 

Number 

Date 

1 

4 May 72 

21 

5 September 72 

2 

5 May 72 

22 

14 September 72 

3 

11 May 72 

23 

30 September 72 

4 

11 May 72 

24 

1 October 72 

5 

12 May 72 

25 

17 October 72 

6 

12 May 72 

26 

16 November 72 

7 

19 May 72 

27 

17 November 72 

8 

not available 

28 

2 January 73 

9 

25 May 72 

29 

10 January 73 

10 

25 May 72 

30 

12 January 73 

11 

5 June 72 

31 

25 January 73 

12 

6 June 72 

48 

22 June 73 

13 

23 July 72 

49 

25 June 73 

14 

24 July 72 

50 

3 July 73 

15 

18 August 72 

51 

5 July 73 

16 

25 August 72 

52 

6 July 73 

17 

28 August 72 

53 

7 July 73 

18 

29 August 72 

54 

4 August 7 3 

19 

29 August 72 

55 

5 August 73 

20 

30 August 72 

56 

10 August 73 



57 

12 August 73 
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(quartz -iodine lamp) was subject to jarring and subsequent movement during scanner installa- 
tion. This fact is obvious from the graphs; not only does the magnitude of signal change, but 
also the changes in overall system response correlate almost 100% between channels. 

An effort was then made to prevent the lamp housing from being again inadvertently jarred. 
Subsequently, we decreased both the frequency of scanner installation and removal. This brings 
us into a third period covering the months of September 1972 through January 1973. Here the 
overall variation in system response is better than that observed in May and June, and the high 
correlation in changes between channels is evident. It is apparent from this signal correlation 
that the variations noted are primarily due to response changes such as dust collecting on the 
scanner mirror or secondary mirrors and lens. This also explains why the short-time response 
(one day or less) is usually much more stable than the long-time response. 

Because a lack of time prevented processing all the data, a period of five months was skipped 
to permit the data from the most recent mission to be compared to data acquired earlier. 
During April 1973 the lamp housing was modified substantially to permit better viewing by the 
scanner mirror; also, a 45-watt bulb was installed in place of the original 100-watt bulb. 

Because of the spectral differences in the smaller lamp, some change in absolute magnitude 
of the signals is to be expected. In general, the variations in lamp voltages are similar to those 
observed with the 100-watt bulb. Except for channel CIO which, for some unknown reason, in- 
dicated excessive variation, scanner signals over a two-month span were stable to within an 
average of 10%. Day-to-day stability tends to be typically better than 2%. 

2.1.3 ESTABLISHMENT OF PROCEDURE FOR WEEKLY MONITORING 

To facilitate a quick computer turn-around of this CALIB. output for the standardization 
tests, a system was initiated for analog -to -digital conversion of the data and operation of the 
CALIB. program on the data. CALIB. channel equivalents appear in Table 2. The idea was to 
minimize the time and cost of constructing the stability history of the M7 scanner. Forms 
were designed which, when properly filled out, will allow for complete data reduction. Table 3 
represents the two forms that, when filled in, will permit complete automatic data processing 
by the program CALIB. Since the standardization tapes are digitized in a set procedure, no 
operator intuition is required. Directions for filling out the forms are detailed in Table 4. A 
sample form, filled out with comments explaining what was done, is shown in Table 5. 

2.2 INTERMEDIATE FREQUENCY SCANNER SIGNAL NOISE (MINUTES TO HOURS) 

2.2.1 PROGRAM DGNSTC* 

The purpose of this program is to perform a more detailed analysis of possible noise prob- 
lems with the M7 scanner system — specifically, higher frequency problems. To facilitate this, 
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TABLE 2. S-CHANNEL EQUIVALENTS FOR PROGRAM CALIB. 


(jgn) 

0.41 - 0.48 
0.46 - 0.49 
0.48 - 0.52 
0.50 - 0.54 
0.52 - 0.57 
0.55 - 0.60 
0.58 - 0.64 
0.62 - 0.70 
0.67 - 0.94 

1.0 - 1.4 
1.5 - 1.8 

2.0 - 2.6 

1.0 - 1.4 

2.0 - 2.6 
0.33 - 0.38 
0.71 - 0.73 


Detector 

Notation 

C12 

CIO 

C9 

C8 

C7 

C6 

C5 

C4 

C2 

D3 

D2 

D1 

B3 

B2 

El 

FI 


S- Channel 
(for "CALIB.") 

16 

15 

14 

13 

12 

11 

10 

9 

8 

7 

6 

5 

4 

3 

2 

1 


*Note: These wavelengths are currently valid for the given 
detector notation. In general, however, care must be taken in 
associating a particular band with a detector notation. 
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TABLE 3. DATA FORM FOR PROGRAM CALEB. 

FORTRAN Coding Form 


GX2B-7327-6 U/M 050* 
Printed in U.SA 


programmer 


caro Electro number 
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TABLE 3. DATA FORM FOR PROGRAM CALIB. 

GX28-7327-6 U/M 050' ‘ 

FORTRAN Coding Farm „ Printed In U.S.A. 


PROGRAM 
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TABLE 4. DIRECTIONS FOR COMPLETING DATA FORMS FOR PROGRAM CALIB, 

Card 1: Supply only if beginning a second or following file set (See last instruction) 

Card 2: Supply FILMAX (see Table 5) 

Card 3: Supply INBIN, UNIT, FILE , NS A (see Table 5) 

Card 4: Supply F DATE, MISSION #, (see Table 5) 

Card 5: Supply COMMENT or TITLE (see Table 5) 

Card 6-18; Fill out 

(1) Multiplex channel => (MPXMUX, order of channels on digital tape, usually 
consecutively numbered 1->X) 

(2) SCHAN => see chart #1 

(3) AIRCHAN => Match previously numbered MPXMUX to original channels 
(’tape channels") 

Fill out as many of the 18 channel blanks as needed, then cross out 
remaining blanks! 

Card 19: Supply as In line 3 above (only FILE &NSA, see Table 5) 

Card 20: Supply new FILE and NS A # (Table 5) 

Card 21: Supply FILE DATE and MISSION # (Table 5) 

Card 22: COMMENT (Table 5) 

Card 23-35: Fill out as in 6-18 above 

Card 36: Fill out as in 19 above 

At this point repeat 19-36 above as needed to cover all files to be entered (flights). 
After last file (last line 19 completed) has been entered — proceed 

Last Card: MODE 1 = $ FINAL $ * (Enter as shown) 

For new file, set start at CARD #1 and supply MODE information as indicated. 
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TABLE 5. COMPLETED SAMPLE DATA FORM 


Card 1 

MODE(l) = $START$ 
♦(ASTERISK) 


Card 2 

CHANNL = 13 
DEFAULT = 13 


FILM AX = 8 
DEFAULT = 1 


LAMT0P = 5 
DEFAULT = 5 


TAPE = $NEW$ • 
DEFAULT = $OLD$ 


MISS0N =12 
DEFAULT = 0 


TBASE 1 = 799 

DEFAULT = 799 

TBASE 2 = 800 

DEFAULT = 800 
♦(ASTERISK) 

Card 3 

IN BIN =777 

DEFAULT = -1 

UNIT = 4 
DEFAULT = NONE 

FILE = 1 
DEFAULT = 0 


(Not used in this sample because there is only one file 
set.) Asterisk signifies End of Instruction; it is used 
after last entry on card. 

Maximum number of channels to be analyzed in any one 
file. (May process fewer than maximum) 

1 s CHANNL * 13 

Maximum number of files to be processed. FILMAX a 1 
(1 file/flight) 

(Plot drops files entered out of sequence. Ex: 1, 2, 4, 8, 
10, 7. 7 is dropped.) 

Number of data points at top of lamp pulse to be used in 
calculations 

0 < LAMT0P £ LAMPE - LAMPB + 1 
see card 4 

(If LAMT0P < 5, lamp bounds LAMPE AND LAMPB 
must be narrowed.) 

$NEW$ = No tape data exists. Present data are used as 
base and put on tape TBASE 1. 

$OLD$ = Tape base exists (by DEFAULT = 799), and is 
merged with processed data (data base compilation at 
799). 

$ PRINTS = Recall of precompiled data from last 
MISS0NX. MISS0N must be specified. Only 3 statements 
needed: 

PRINT = $PRINT$ 

MISS0N = X 
MODE = $ FINAL? 

Only for a preceding instruction of $ PRINT?. 

Number of mission up to which the accumulated data is 
plotted and tabulated (including above-specified mission 
data). 

Used for a preceding instruction of $NEW$ (BIN number 
of tape). 

BIN number of tape from another data base 
End of Instruction 


BIN number of input tape to be mounted 
(never = 0); if = -1, no mounts. 

Logical unit number of tape drive on which input tape is 
mounted (only 2, 3, 4, 8, 9) 

File on input tape to be processed; if < 1, no positioning. 
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TABLE 5. COMPLETED SAMPLE DATA FORM (Continued) 


NSA = 1 

- 225 
= 1 
= 1 
= 358 
= 1 

DEFAULT = 0 
* (ASTERISK) 


Starting line number 
Ending line number 
Line number increment 
Starting point number 
Ending point number 
Point number increment 

End of Instruction 


Card 4 

FDATE (1)= 05, 25, 73 
DEFAULT = 00, 00, 00 


MISS0N =12 
DEFAULT = 0 

LAMPB = 10 
DEFAULT = 10 


LAMPE = 60 
DEFAULT = 60 


DARKB = 185 
DEFAULT = 185 


DARKE = 225 
DEFAULT = 225 


S0LARB =310 
DEFAULT = 310 


S0LARE = 355 
DEFAULT = 355 

Card 5 

COMMNT (1) = $S9 t 
BLACK HILLS, 0745, $ 

♦(ASTERISK) 


Date of flight associated with file 
3 variables — ► FDATE (1) - mo - 0 < (1) s 12 

FDATE (2) - dy - 0 < (2) <31 

FDATE (3) - yr - 0 < (3) ^ 99 

Number of mission associated with file being 

entered. MISS0N * 0 

Number of reselm denoting beginning of lamp 
pulse (Note: DEFAULT ) 

LAMPB a NA (see "NSA,” card number 3) 

Number of reselm denoting the end of lamp pulse 
(Note: DEFAULT ) 

LAMPE ^ NB (see "NSA, 11 card number 3) 

Number of reselm denoting the beginning of dark 
signal (Note: DEFAULT ) 

DARKB * NA 

Number of reselm denoting the end of dark signal 
(Note: DEFAULT ) 

DARKE s NB 

Number of reselm denoting the beginning of solar 
signal (Note: DEFAULT ) 

S0LARB & NA 

Number of reselm denoting the end of the solar 
signal (Note: DEFAULT ) 

S0LARE ^ NB 

Any 4 8 ^character comment 

to be printed as given between $ signs 

End of Instruction. (Includes cards 4 & 5.) 


Card 6 

MPXMUX = 1 Number of multiplexer channel 

1 * MPXMUX s 13 

SCHAN = 8 Number of spectral channel MPXMUX represents 

(See chart ) 

0 - SCHAN ^ 16; if = 0, MPXMUX is ignored in processing* 
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TABLE 5. COMPLETED SAMPLE DATA FORM (Concluded) 


AIRCHN = 8 
♦(ASTERISK) 

Card 7 

MPXMUX = 2 
SCHAN - 9 
AIRCHN =13 
♦(ASTERISK) 
Card 8 
MPXMUX = 3 
SCHAN = 10 
AIRCHN = 6 
♦(ASTERISK) 

Card 9 

MPXMUX - 4 
SCHAN = 11 
AIRCHN = 5 
♦(ASTERISK) 

Card 10 
MPXMUX = 5 
SCHAN = 12 
AIRCHN = 4 
♦(ASTERISK) 

Card 11 
MPXMUX = 6 
SCHAN = 13 
AIRCHN = 3 
♦(ASTERISK) 

Card 12 
MPXMUX - 7 
SCHAN = 15 
AIRCHN = 2 
♦(ASTERISK) 

Card 13 
MPXMUX = 8 
SCHAN = 16 
AIRCHN = 1 
♦(ASTERISK) 

Cards 14-18 
Card 19 


Number of aircraft channel that MPXMUX represents 
End of Instruction 


End of Instruction 


End of Instruction 


End of Instruction 


End of Instruction 


End of Instruction 


End of Instruction 


End of Instruction 

For this example, MPXMUX cards 14 through 18 were 
not needed. 


File = 1 

NS A = 1, 225, 1, 1, 358, 1 NSA is same as card 3. 
♦(ASTERISK) End of Instruction 


Last Card 

MODE (1) = $ FINAL? f 


Close files. 


t NOTE: When a "Data Set" is closed (file set started with FlLMAX card and 
finished with MODEl = $ FINAL $), further program operation requires all cards 
and steps to be repeated. 

For a scries of files within a data set, you only need the following: 

(1) FILE, NSA 

(2) F DATE, MISS0N 

(3) COMMNT 

(4) MPXMUX, SCHAN, AIRCHN 
LAST, FILE, NSA 

Repeat 
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the program has been made as general as possible and designed for rapid execution. It provides 
a technique to measure the effect of varying scan rate or tape recorder speed on digitized data. 
Although the statistics do not point to the cause of the apparent speed change, they do show that 
a change has occurred in the system and how drastic it might be. 

The user is also able to determine the mean and standard deviation (1) of the signal level 
within any specified area, and (2) of the signal level of another area relative to the first. For 
example, when used with the lamp reference source, the user can measure changes in lamp ra- 
diance relative to a dark area within the scanner. In addition, plots of both these means and 
standard deviations versus scanline number may be obtained. With these plots, the user will be 
able to observe the detailed quality of changes in the signal level of an area in the scanline along 
the flight path. For instance, the user could look for cyclical fluctuations in an ambient area. 

A measure of speed changes in the system is accomplished by analyzing some fixed-position 
signal in the scanline such as the lamp reference pulse. The mean position of the middle of the 
lamp pulse is measured relative to the beginning sync pulse for each scanline. The standard 
deviation from this mean shows how much jitter there is in the digitized image of the lamp pulse. 

Differing only in output detail, the other two facilities of the program do virtually the same 
thing. One feature simply determines the mean and standard deviation of the signal level (either 
absolute or relative to another area) in a specified area of the scanline over the entire data set. 
The other feature calculates the mean and standard deviation of the signal level (either absolute 
or relative to another area) in a specified area of each scanline. Plots are made of these sta- 
tistics versus the scanline number from which they were calculated. Appendix B includes a 
complete description of the program DGNSTC. 

2.2.2 EVALUATION OF DATA SET 

Measurements were made on M7 scanner data to determine, quantitatively, the intermediate 
frequency noise (or drift) that one might encounter in processing multi spectral scanner data 
over short periods of time. The particular data set chosen for this study was obtained over 
Eaton County, Michigan, on 25 August 1972, as part of the ERTS underflight program. This 
mission was chosen specifically because it satisfied two requirements: (1) it had been obtained 
recently enough to qualify the scanner configuration and performance as being typical of the pre- 
sent system; and (2) it contained several long (12 -minute) data passes encompassing a total 
period of approximately one hour — this enabled us to obtain a near-continuous measurement of 
the short-term drift over a manageable time period. 

A total of five passes covering one hour and twelve minutes (see Table 6) were analyzed by 
DGNSTC. Since five full 12 -minute data runs would require a prohibitive amount of processing, 
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TABLE 6. LOG OF EATON COUNTY, MICHIGAN, FLIGHTLINES 


Line 

Date 

GMT 

Start 

GMT 

Stop 

Altitude 
(1000 ft) 

Aircraft 

Heading 

(deg) 

1 

25 August 72 

16:20:30 

16:30:00 

5 

360 

2 

25 August 72 

16:41:00 

16:51:00 

5 

180 

3 

25 August 72 

16:54:30 

17:06:00 

5 

360 

4 

25 August 72 

17:09:00 

17:18:30 

5 

180 

5 

25 August 72 

17:22:00 

17:37:30 

5 

360 
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five segments of data (equally spaced from beginning to end) were digitized for each run. Each 
segment required 15 seconds of actual data collection time out of the original 12 -minute run. 

Only a region of the dark level (base line) and the lamp was observed. Output from the program 
DGNSTC. was obtained on this data; the results are summarized in Tables 7, 8, and 9 
(X = 0.41-0,48, 0.55-0.60, and 2. 0-2. 6 pm, respectively). 

Only three spectral channels were analyzed since a considerable amount of computer time was 
still required for processing each channel. It should be noted here that the data were averaged 
(smoothed) over four lines before processing; hence, the noise figures (expressed as a standard 
deviation) listed in the tables do not represent the values on the original analog tape. It can be 
shown that smoothing n-lines has the effect of decreasing the noise by / n. Therefore, the actual 
noise values on tape are twice those indicated in the tables (smoothing over four lines). The 
values were left unchanged, though, because much of the computer processing being done uses 
the concept of smoothing to reduce noise. Hence, these values are directly related to current 
techniques. 

The mean value shown in Tables 7, 8, and 9 for the dark level and lamp pulse are the aver - 
age s of the means obtained from the region of interest in each line. The standard deviations are 
the deviations of the mean of each line about the average mean for the whole segment. For the 
dark level, the region digitized was midway between the reference lamp pulse and the sky refer- 
ence pulse (see Fig. 4). As the tables indicate, in the solar -reflective region the value of the 
dark level per line (averaged over 25 points) deviates only 0.1% from the mean dark level value 
for the whole 72 -minute data set. In the near-IR region where the dark level is adjusted man- 
ually before each run via dc-offset pots, the deviation is a remarkably small 0.3%. Note that the 
deviation over the whole 72 -minute set is no larger than that for the 12 -minute runs. 

The values shown for the reference lamp were calculated in a similar manner. The region 
used for averaging the lamp signal was 15 points wide and covered dark levels on both sides of 
the pulse. The average value obtained for the lamp is analogous to integration of the pulse and 
not to the peak value. Standard deviations of the lamp pulse varied from 2,3% of the mean in 
the blue channel (noisiest channel) to 1.3% in the near-IR. 

In general, then, when smoothing four lines and averaging 15-25 points per line, the stan- 
dard deviation of a signal from line-to-line (assuming a constant input) will generally not ex- 
ceed 2.5% (worst case) of the mean value over all lines collected in a period of time as long as 
one hour and twelve minutes. Of course, in many cases, smoothing between lines and averaging 
in a line is not possible. For this instance, statistical theory leads us to a standard deviation 
that will be Vi" x VTI larger than that obtained when smoothing over four lines and averaging 
over 15 points per line. Hence, the standard deviation from line to line with no averaging of any 
kind would be as follows: 
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TABLE 7. MEANS AND STANDARD DEVIATIONS OF 
DARK LEVEL AND LAMP FOR FIVE DATA SETS 
OVER EATON COUNTY. A = 0.41-0.48 (im 


Dark Level Lamp 


Line No. Seg 

ment No. 

X 

a 

X 

a 

1 

1 

140.18 

0.16 

13.45 

0.22 


2 

140.05 

0.15 

13.30 

0.20 


3 

140.07 

0.16 

13.34 

0.20 


4 

140.40 

0.22 

13.35 

0.50 


5 

140.32 

0.21 

13.47 

0.46 


AVE 

140.20 

0.18 

13.38 

0.34 

2 

1 

139.29 

0.19 

13.52 

0.37 


2 

139.50 

0.21 

13.28 

0.55 


3 

139.63 

0.20 

13.15 

0.50 


4 

139.80 

0.19 

13.10 

0.33 


5 

139.83 

0.18 

13.15 

0.25 


AVE 

139.61 

0.19 

13.24 

0.41 

3 

1 

140.82 

0.18 

13.15 

0.28 


2 

140.86 

0.18 

13.18 

0,26 


3 

140.90 

0.16 

13.21 

0.25 


4 

140.87 

0.17 

13.11 

0.23 


5 

140.91 

0.15 

13.15 

0.21 


AVE 

140.87 

0.17 

13.16 

0.25 

4 

1 

140.91 

0.16 

13.15 

0.22 


2 

140.92 

0.17 

13.15 

0.23 


3 

140.90 

0.14 

13.14 

0.17 


4 

145.02 

0.15 

13.34 

0.21 


5 

145.03 

0.14 

13.64 

0.18 


AVE 

142.56 

0.15 

13.28 

0.20 

5 

1 

145.21 

0.17 

13.46 

0.32 


2 

145.16 

0.19 

13.35 

0.29 


3 

145.17 

0.19 

13.40 

0.29 


4 

145.21 

0.18 

13.37 

0.34 


5 

145.22 

0.18 

13.35 

0.34 


AVE 

145.19 

0.18 

13.39 

0.32 

Mean Total 


141.69 

0.18 

13.29 

0.31 


Percent Deviation 0.1% 2.3% 
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TABLE 8. MEANS AND STANDARD DEVIATIONS OF 
DARK LEVEL AND LAMP FOR FIVE DATA SETS 
OVER EATON COUNTY. X = 0.55-0.60 fim 


Dark Level Lamp 


Line No. Se^ 

$ment No. 

X 

a 

X 

CT 

1 

1 

144.05 

0.15 

37.01 

0.19 


2 

143,90 

0.15 

36.74 

0,23 


3 

143.91 

0.15 

36.73 

0.21 


4 

143.92 

0.18 

36.82 

0.25 


5 

143.95 

0.18 

36.82 

0.25 


AVE 

143.95 

0.16 

36.82 

0.23 

2 

1 

144.05 

0.15 

36.74 

0.20 


2 

144.09 

0.16 

36.77 

0.23 


3 

144.07 

0.16 

36.63 

0.20 


4 

144.04 

0.15 

36.50 

0.20 


5 

144.01 

0.15 

36.61 

0.21 


AVE 

144.05 

0.15 

36.65 

0.21 

3 

1 

143.93 

0.15 

36.67 

0.20 


2 

143,91 

0.15 

36.58 

0,21 


3 

143.91 

0.16 

36.54 

0,25 


4 

143.86 

0.15 

36.60 

0.26 


5 

143.87 

0.15 

36.65 

0.27 


AVE 

143.90 

0.15 

36.61 

0.24 

4 

1 

143.93 

0.14 

36.62 

0.25 


2 

143.88 

0.16 

36.43 

0.25 


3 

143.83 

0.15 

36.36 

0.25 


4 

143.89 

0.16 

36.49 

0.44 


5 

143.90 

0.14 

36.77 

0.29 


AVE 

143.89 

0.15 

36.53 

0.31 

5 

1 

144.02 

0.13 

36.77 

0.22 


2 

143.97 

0.14 

36.43 

0.22 


3 

143.97 

0.14 

36.53 

0.22 


4 

144.01 

0.14 

36.62 

0.22 


5 

144.05 

0.15 

36.76 

0.26 


AVE 

144.00 

0.14 

36.62 

0.23 

Mean Total 


143.96 

0.15 

36.65 

0.24 


Percent Deviation 0.1% 0.7% 
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TABLE 9. MEANS AND STANDARD DEVIATIONS OF 
DARK LEVEL AND LAMP FOR FIVE DATA SETS 
OVER EATON COUNTY. \ = 2.0 -2. 6 jim 


Dark Level Lamp 


Line No. 

Segment No. 

X 

a 

X 

0 

1 

1 

160.24 

0.42 

47.89 

0.65 


2 

159.00 

0.40 

47.60 

0.57 


3 

159.74 

0.39 

47.46 

0.61 


4 

159.15 

0.39 

47.55 

0.58 


5 

160.48 

0.40 

47.03 

0.91 


AVE 

159.72 

0.40 

47.51 

0.66 

2 

1 

163.13 

0.43 

47.28 

0.58 


2 

163.53 

0.40 

47.37 

0.61 


3 

163.43 

0.46 

47.22 

0.57 


4 

163.77 

0.35 

47.07 

0.52 


5 

164.55 

0.35 

47.16 

0.59 


AVE 

163.68 

0.40 

47.22 

0.57 

3 

1 

163.87 

0.35 

47.18 

0.53 


2 

163.29 

0.37 

47.16 

0.54 


3 

163,64 

0.32 

47.17 

0.60 


4 

162.75 

0.37 

47.26 

0.54 


5 

161.93 

0.34 

47.21 

0.45 


AVE 

163.10 

0.35 

47.20 

0.53 

4 

1 

162.79 

0.44 

47.17 

0.59 


2 

163.62 

0.35 

47.03 

0.67 


3 

163.93 

0.38 

47.16 

0.57 


4 

165.28 

0.47 

47,00 

0.63 


5 

165.19 

0.37 

47.28 

0.52 


AVE 

164.16 

0.40 

47.13 

0.60 

5 

1 

164.81 

0.53 

47.15 

0.66 


2 

164.29 

0.61 

47.08 

0.83 


3 

163.51 

0.52 

47.02 

0.66 


4 

163.07 

0.37 

47.00 

0.55 


5 

162.92 

0.38 

47.13 

0.60 


AVE 

163.72 

0.49 

47.08 

0.67 

Mean Total 


162.88 

0.41 

47.23 

0.61 

Percent Deviation 

0.3% 


1.3% 
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standard deviation nQ average = ^4 x Vl5 x 2.3% (X - 0.41 0.48 pm) 


= 17.8% {worst case) 

This would be the variation, line-to-line, about the mean value over all lines expected in a 
scanner signal obtained from a single resolution element over a data run 12 to 72 minutes in 
length. 

2.3 HIGH FREQUENCY NOISE CHARACTERISTICS (SECONDS-TO- MINUTES) 

Thus far, noise measurements have been made on M7 scanner data to determine the signal 
variation to be expected over months, days, and hours. This kind of information is valuable in 
terms of analyzing signature extension capabilities over various time periods. Another aspect 
to the noise problem is the higher frequency noise that occurs within a single scanline. Knowl- 
edge of this noise (which is akin to NEAL) would permit the data processor, for example, to 
evaluate the potential for resolving the spectra of two similar objects. But, before discussing 
the actual measurement of this noise, several prefatory comments should be made about the 
M7 system in general. 

As in the thermal channel {discussed later), the noise observed in solar -reflective channels 
is almost always tape-recorder-limited; that is, the dynamic range of the tape recorder {Mincom 
PC -500) is almost always less than the dynamic range of signals encountered in viewing the 
terrain (including detector and amplifier noise). This recorder dynamic range is 39.5 dB, 
peak-to-peak, which results in a maximum signal -to-noise ratio (S/N) of 89:1. This assumes 
that the entire range of the recorder is being used and that the signal is indeed tape-recorder- 
limited (true most of the time). The one possible exception occurs when the highest gains are 
used over iow reflecting objects such as water or during low sun-angle situations. In this case, 
detector or preamplifier noise becomes the limiting factor. 

With this maximum signal-to-noise ratio in mind, data were analyzed over Run 1 of the 
same mission discussed in Section 2.2. The average dark level values (in bins) over 195 lines 
and 25 pts/line for each of the five segments are recorded in Table 10. Also tabulated is the 
average standard deviation per line j^ ne ) f° r eac h of the segments and their deviations (%-)• 

The same wavelengths discussed earlier were also studied here. The most important thing 
to note from these data is that in the solar -reflective region the noise in a single scanline (denoted 
by o” line ) is indeed equal to the tape recorder limit. 

To see this, one must be aware that in digitizing the data the maximum number of bins 
available to cover the full dynamic range of the recorder is 512. Since the average dark level 
(0.41-0.48 pm) is at 140.20, the number of bins available to the recorder for a maximum signal 
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TABLE 10. DARK LEVEL MEANS AND STANDARD DEVIATIONS 
FOR ONE SEGMENT OF DATA OVER EATON COUNTY 


Dark Level 


X 

Run No. 

Segment No. 

X Line 

G Line 

O’- 

er 

0.41-0.48 |im 

1 

1 

140.18 

0.39 

0.12 



2 

140.05 

0.35 

0.14 



3 

140.07 

0.39 

0.11 



4 

140.40 

0.46 

0.08 



5 

140.32 

0.45 

0,09 



AVE 

140.20 

0.41 

0.11 

0.55-0.60 4m 

1 

1 

144.05 

0.38 

0.13 



2 

143.90 

0.38 

0.12 



3 

143.91 

0.40 

0.13 



4 

143.92 

0.40 

0.11 



5 

143.95 

0,39 

0,13 



AVE 

143.95 

0.39 

0.12 

2.0 “2.6 4 m 

1 

1 

160.24 

1.43 

0.43 



2 

159.00 

1.36 

0.32 



3 

159.74 

1.37 

0.34 



4 

159.15 

1.37 

0.32 



5 

160.48 

1.36 

0.46 



AVE 

159.72 

1.38 

0.38 
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input is 512 - 140.20 = 371.80. Hence, the corresponding signal -to-noise ratio is 907.8 : 1 
(371.8 -7 0.41). In obtaining the data, however, 100 points were averaged to arrive at the fig- 
ures shown in Table 10. Therefore, the noise in each line (with no averaging) is actually VT00 
higher than indicated. This reduces the S/N to 90.8:1. The corresponding values for the other 
two channels are as follows: 

for A = 0.55-0.60 pm, S/N = 94.4 
A - 2 . 0-2 . 6 pm, S/N = 25 .5 

The near infrared channel does not appear to be tape-recorder -limited, although the cause of 
the noise is at this time unknown. 
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3 

DARK-LEVEL DETERMINATION 

One of the more seemingly insignificant decisions made by processors of multispectral 
data is the choice of a baseline region for dark -level correction. There are several regions 
of darkness viewed by the rotating scan mirror in one 360° swing, and in the past it has been 
more or less an arbitrary choice by the processor as to which region to use. One of the 
goals of the study was to determine if all these baseline regions were indeed equally dark 
and, if not, what region was optimal. 

Data were examined from a flight flown on 20 March 1973, at 15:30 GMT, over Washtenaw’ 
County specifically to locate an optimal dark level region. Only five spectral channels were 
digitized so that maximum spatial resolution could be obtained. The entire scanline was not 
digitized because of gate limitations on the analog -to 'digital equipment; nearly 1000 sample 
points out of a possible 1333 (at 4.7 mr/sample and 80K samples/sec) were obtained. The 
300 and some points not digitized are located in the region between the sky reference and the 
fixed sync (see Fig. 6). Included in this area are the hot plate and UV reference lamp. It is 
possible that this area contains a region of even darker base line than that determined from 
the remainder of the data; since the present system cannot digitize data in this area, however, 
it will be neglected here. 

Shown in Figure 5 is a plot of signal (in bins) versus reselm location for a single scanner 
line. The discrete level fluctuations in the data are the results of quantization into bins of 
voltage levels bn the analog tape during the analog -to -digital conversion process. Each bin 
is equal to 0.01953 volt. Tables 11 and 12 contain the tabulated results of a statistical 
analysis performed on nine regions of the base line data (see Fig, 6). Table 11 shows the 
mean (in bins) and standard deviation for all five channels and one scanline. Table 12 gives 
the same figures averaged over 100 scanlines. As the tables show, the best region for dark 
level correction is the area just before the beginning of the video (Region I). In terms of 
physical location inside the scanner housing, this region corresponds to a blackened area 80° 
from nadir. Unfortunately, under severe aircraft roll conditions, the roll -corrected sync 
pulse can move into this region. During turbulent weather, then, Region I should not be used. 

However, an alternative is to use Regions VIII or IX, which have values very near to that 
of Region I. These areas are part of one broad region located between the reference lamp and 
the sky reference. The blackened cold -plate thermal reference is included in this area. The 
pulse located at reselm number 715 is an SCR noise spike which comes from a solid state 
DC-to-AC inverter used to supply 110 V, 60 Hz to the equipment. This noise spike will al- 
ways be found somewhere in the data (at least for the present), although efforts are being 
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TABLE 11. SUMMARY OF DARK LEVEL AVERAGE VALUES FOR DIFFERENT SPATIAL AND SPECTRAL REGIONS IN A 
SINGLE SCAN PERIOD. Average values are in bins (1 bin = 0.01953 V). Numbers in parentheses are normalized values. 




X - 0.41-0.48 jim 

X = 0.50-0.54 

fim 

X = 0.58-0.64 

/am 

X = 0.67-0.94 

ji m 

X - 1. 0-1.4 

fim 

REGION 

LOCATION 

X 

a 

X 

a 

X 

cr 

X 

a 

X 

0 



(0.41-0.48) 


(0.50-0.54) 


(0.58-0.64) 


(0.67-0.94) 


(1. 0-1.4) 


I 

1-58 

153.31 

(1.000) 

0.50 

150.81 

(1.000) 

0.66 

149.05 

(1.000) 

0.47 

149.17 

(1.000) 

0.53 

154.60 

(1.000) 

1.22 

II 

530-570 i 

1 

154.83 

(1.010) 

0.49 

152.37 

(1.010) 

0.53 

150.80 

(1.012) 

0.40 

155.17 

(1.040) 

0.49 

156.54 

(1.013) 

0.50 

III 

600-620 

154.43 

(1.007) 

0.49 

151.86 

(1.007) 

0.47 

150.52 

(1.010) 

0.50 

154.29 

(1.034) 

0.45 

155.24 

(1.004) 

0.43 

IV 

640-660 

153.95 

(1.004) 

0.49 

151.67 

(1.006) 

0.56 

150.00 

(1.006) 

0 

153.28 

(1.028) 

0.45 

155.24 

(1.004) 

0.43 

V 

680-700 

153.81 

(1.003) 

0.39 

151.33 

(1.003) 

0.47 

149.95 

(1.006) 

0.21 

152.57 

(1.023) 

0.49 

155.05 

(1.003) 

0.21 

VI 

600-700 

154.01 

(1.005) 

0.52 

151.60 

(1.005) 

0.55 

150.16 

(1.007) 

0.39 

153.35 

(1.028) 

0.77 

155.22 

(1.004) 

0.44 

VII 

725 -750 

153.56 

(1.002) 

0.57 

151.22 

(1.003) 

0.42 

149.70 

(1.004) 

0.46 

150.44 

(1.009) 

0.50 

154.59 

(1.000) 

0.49 

VIII 

800-825 

153.38 

(1.000) 

0.49 

151.08 

(1.002) 

0.47 

149.58 

(1.004) 

0.49 

150.00 

(1.006) 

0.28 

154.42 

(1.000) 

0.49 

IX 

875-900 

153.42 

(1.001) 

0.49 

150.92 

(1.001) 

0.47 

149.29 

(1.002) 

0.45 

149.73 

(1.004) 

0.52 

154.19 

(0.997) 

0.39 
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TABLE 12. SUMMARY OF DARK LEVEL AVERAGE VALUES FOR DIFFERENT SPATIAL AND SPECTRAL REGIONS 
AVERAGED OVER 100 SCANLINES. Average values are in bins (1 bin = 0.01953 V). 

X = 0.41-0.48 fim X = 0.50-0.54 X = 0.58-0.64 /am X = 0.67 -0.94 (im X = 1.0 -1.4 jam 


REGION 

X 

cr 

X 

a 

X 

a 

X 

a 

X 

a 


(0.41-0.48) 


(0.50-0.54) 


(0.58-0.64) 


(0.67-0.94) 

l 

(1.0-1. 4) 


I 

152.47 

0.62 

150.01 

0.56 

148.36 

0.38 

149.07 

0,43 

153.77 

1.06 

II 

153.61 

0.53 

151.18 

0.49 

149.71 

0.33 

154.23 

0.40 

155.30 

0.43 

VIII 

152.58 

0.54 

150.13 

0,48 

148.60 

0.34 

149.60 

0.33 

153.54 

0.39 

IX 

152.50 

0.55 

150.06 

0.48 

148.52 

0.34 

149.44 

0.36 

153.68 

0.39 
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made to move it to an area of insignificance. For this reason, a standard region to use for 
dark level correcting is the one between regions VIII and IX which corresponds to the posi- 
tion of the cold plate. Since the SCR noise spike is always kept out of the hot and cold plate 
regions, there need be no worry about this pulse affecting the dark level value. From the 
standpoint of time, this dark level region occurs 1.326 msec before the sky reference. Any 
width up to approximately 30 reselms (4.7 mr/sample) can be used. It is expected that all data 
digitized in the future will use this region for dark level correction. 
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4 

M7 RADIANCE CALIBRATION 


4.1 RADIANCE CALIBRATION 

Eleven radiance calibration tapes have been generated since the M7 scanner was first 
put into service in June 1971. Their dates are listed in Table 13. Until now, only a portion 
of these had been digitally analyzed — using a variety of techniques. In order to understand 
more fully the changes that have occurred in these calibration data sets and to simplify the 
data reduction procedure, the program LAMRAD. was written. Designed to handle M7 scanner 
radiance tapes exclusively, this program calculates the radiance calibration constants (as a 
function of lamp current) used to provide radiance information for multispectral processing. 

A complete derivation and explanation of the procedure used to calibrate the M7 scanner is 
given in Appendix C. In summary, the radiance of a target on the ground is given by 



where L T (X) = radiance of target, assuming a Lambertian surface, for each channel 
= voltage of target on ground 
= peak lamp voltage during target data collection 
= peak lamp voltage during calibration 
Vpef = voltage off reference panel during calibration 
E M = spectral irradiance of quartz -iodine lamps 
P re j = directional reflectance of reference panel 


It is the quantity 


k cai>> ' V 


V L 'WV > > 


ref 


that is calculated by the program LAMRAD. 


4.2 TECHNIQUES FOR LAMP PULSE MEASUREMENTS 
4.2.1 SUMMARY OF TECHNIQUES 

Critical to the equation above for &q AI j W are the values of and V L , This is true be- 
cause Vj^and V L , of all the terms in that equation, are by far the most difficult to measure con- 
sistently. In general, V L is represented by a narrow pulse that on occasion can exhibit jitter equal 
to its pulse width. This makes it very difficult for the data processor to measure its peak value. 
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TABLE 13. DATES OF SPECTRAL 
RADIANCE CALIBRATION 


Calibration No. 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


Date 

20 August 71 

6 June 72 

7 June 72 
16 June 72 
12 July 72 

5 September 72 
24 January 73 
26 January 73 
3 April 73 

5 June 73 

6 August 73 
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One of the tasks of this study was to attempt to improve the technique for measuring this 
voltage. 

Thus far, two methods have been developed for calibrating the signal pulse of the lamp in 
the M-7 scanner using the program LAMRAD. Both methods result in a voltage value for the 
peak of the signal. One method, called the "mean level method’* and used exclusively up 
until this time, has certain drawbacks which a second method, called the ''integral method, " 
will hopefully eliminate. 

The mean level method of determining signal height (voltage) is calculated in the follow- 
ing way (refer to Fig. 7), The program will first select the point on the pulse .curve which 
has the curve's peak voltage "P M and then determine at which reselms the half -value of P 
(P/ 2) intersects the curve. There will always be only two intersections: one at the pulse 
rise and one at the pulse drop. This is dictated simply by the shape of the curve. The two 
reselm locations mentioned will be called x and y, respectively. The program will then aver- 
age these values and determine another reselm value (x 4- y)/2. At this point, LAMRAD. will 
determine the voltage of the curve at the point (x + y)/2 and take the average voltage of it and 
other consecutively preselected reselms on either side of (x + y)/2. Figure 7, for example, 
shows only two reselms selected on either side of (x + y)/2 and averaged to obtain the mean 
level. Thus far, this has been the configuration used in LAMRAD. The curve values taken at 
the reselm points and then divided by (N + 1) gives the mean level. The general equation for 
the mean level is: 



where Z may have only integral values and N can only be even. This is to insure that no frac- 
tional reselm values will exist and that one may select equal numbers of reselm samples on 
either side of the (x + y)/2 reselm, thus allowing for a balanced mean. This method of calcu- 
lating a mean voltage height of the pulse has some drawbacks, one of which is illustrated by the 
abnormally skewed curve in Fig. 7. If the curve should happen to be shaped oddly, the value 
for the height could be off true as shown. Another problem comes up in determining how many 
reselms to average on either side of (x + y)/2. We want to average as much of the pulse top as 
possible without going off the sides (the drop and rise usually occur within a reselm or so). 
Once this value is calculated, signal shift cannot be tolerated because a shift of a reselm or 
two would give values at the bottom of the curve, resulting in a "height" far below what actually 
exists. 

The second means of determining signal level is the "integral method.” As the title sug- 
gests, this method involves determining the area under the pulse in calculating the pulse 
height. With this method, instead of trying to pinpoint where the pulse and its peak will appear 
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Pulse is Abnormally Skewed for Illustrative Purposes 


p/2 Reselm Values 
(= (X + y/2)) 


FIGURE 7. MEAN-LEVEL METHOD OF PULSE-VOLTAGE CALCULATION 
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(reselm limits), one may integrate the curve within wide limits (to allow for reselm position 
shift) and then divide the area value by an experimentally predetermined pulse "width con- 
stant" to obtain the mean value of the signal. The "width constant," which is a value deter- 
mined over a period of time, is continually updated to offer accurate calibration values at 
any time. The "width constant" is simply determined by dividing the integrated area by the 
height of the pulse, as calculated in the first method (see Fig. 8). 

Certain conditions must exist if this system of lamp calibration is to offer a more accu- 
rate way of determining signal value. If the shape of the pulse radically changes, the "width 
constant" will be inaccurate, producing incorrect lamp voltage values (although restrictions 
on pulse shape change are certainly no more imposing than in the "mean level method"). 
Generally, the only time the signal can change shape is after some type of movement has 
occurred within the lamp or its housing (observed experimentally). In this case, the width 
constant must be re-calculated to give accurate pulse voltages again. The "width constant" 
must be a time-weighted mean to avoid the accuracy problems encountered in determining the 
signal voltage by the "mean level method." Finally, the integral limits must be set far enough 
apart so that the integral will accurately tolerate the pulse shift that the "mean level method" 
was incapable of handling. The integral limits should be set such that the curve closely 
approaches zero at those points. 

4.2.2 EXPERIMENTAL DETERMINATION OF WIDTH CONSTANT 

LAMRAD. has been executed for several calibration runs. The "Lamp Mean" and the 
"Lamp Integral Mean" were among the data calculated by the program. By taking the values 
of the lamp integral mean and dividing by their respective lamp mean values, we obtained a 
width constant. After looking at the width constants resulting from some seven calibration 
runs, we then came to several conclusions. 

First, the wavelength band in which the signal was observed had an effect on the width 
constant. For all thermal values (1.0pm and above), the width constant varied so radically 
and the deviations were so great that no real attempt was made to obtain a width constant for 
them. The remaining width constants calculated for the various detector bandwidths had such 
small deviations between themselves that we could divide them into just two groups. One 
group used the signals emitted from all detectors operating in the 0.41-0.52 jam range, while 
the other group contained all signals in the 0.50-0.94 pm. As Table 14 reveals, deviations of 
width constants within these two groups for any one run are on the order of 0.8%. Therefore, 
two width constants should be calculated, one for each of the two groups. These two values 
are calculated at the bottom of Table 14. 
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FIGURE 8. DETERMINING PULSE-WIDTH CONSTANTS (K w ) 
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TABLE 14. WIDTH CONSTANTS FORA SINGLE RADIANCE CALIBRATION. (Calculated with 

4.5 -6.0 amp lamp current only) 


Detector 

Bandwidths 

(M m ) 


0.41 - 0.48 

16.0837 

CM J 

| GROUP 2a 

0.46 - 0.49 

16.3078 

§ Mean = 16.2554 ± 0.1244 

] GROUP 2b 

0.48 - 0.52 

16.3747 

ps 

o 

i 

0.50 - 0.54 

16.5327 



0.52 - 0.57 

16.5803 


Mean = 16.4269 ± 0.1537 

0.55 - 0.60 

16.5558 

■h 

(Notice: worse deviation than 
either Group 1 or 2 alone) 

0.58 - 0.64 

16.5238 

£> Mean = 16.5127 ± 0.0750 


0.62 - 0.70 

16.3501 

- O 
£ 


0.67 - 0.94 

— 

16.5335 ' 

j 

1 


1.0 - 1.4 

22.7585 

Mean = 19.92 ± 2.8384 

i 

Disregard (deviation too 
) great) 

1.5 - 1.8 

17.0816 

l 

j 




Final Mean Width Constants (Runs 4-8) for Groups 1 and 2 
using currents of 4.5 A and greater 

GROUP 1 
GROUP 2 

0.41 - 0.48 
0.46 - 0.49 
0.48 - 0.52 

0.50 - 0.54 ^ 
0.52 - 0.57 
0.55 - 0.60 
0.58 - 0.64 
0.62-0.70 
0.67 - 0.94 J 


Mean = 16.3458 ± 0.0445 
Mean = 16.5719 ± 0.0554 
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Secondly, the width constant does not really remain very constant until calculated for 
lamp signals using 4*5 amps or more of current. Width constants calculated for signals using 
only 3.5 or 4.0 amps varied greatly. One way these signals might be made more useful is by 
increasing the gain at the time of calibration, instead of leaving the gain constant as is 
usually done. An increased signal gain may be able to bring the pulse level far enough out of 
the noise to allow more accurate calibrations to be made. Thus far, all experimental width 
constant data and their deviations have been calculated using lamp currents of only 4,5 amps 
and greater. 

Lastly* as previously mentioned, changes in signal shape will affect the value of the width 
constant (the ratio of lamp integral mean to lamp mean will change). Also, the only way in 
which the signal shape can change is if there is some type of alteration made in lamp position 
(such as might result in housing movement). Between calibrations #8 and //9, a change was 
made affecting the position of the lamp’s housing; the expected changes in width constants 
also occurred and are compared in Table 15. (Only signals in three spectral bands were 
analyzed.) New width constants must be calculated for calibration runs after #8. Calibration 
ii 9 gives an indication of what these values might be, but more data are needed to determine 
a good mean width constant for these new calibrations. 

4.3 PROGRAM LAMRAD. 

As mentioned earlier, the primary purpose of this program is to calculate a set of statis- 
tics from the calibration data set on the M7 scanner. These statistics are used to determine 
the radiance of the lamp reference source at various current levels* The statistics calculated 
by this program are; 

(1) the mean peak levels of the lamp and panel pulses (see Appendix C) relative to some 
specified dark area 

(2) the mean integrals of the lamp and panel pulses after dark -level normalization and 
pulse width adjustment 

(3) the standard deviations of the adjusted integrals of the lamp and panel pulses 

(4) the mean of the ratio of the adjusted lamp integral to the adjusted panel integral 
multiplied by a unique factor for each spectral channel 

(5) the standard deviation of the integral ratio 

(6) the mean integral of the lamp and panel pulses after dark -level normalization 

(7) the mean pulse widths for the lamp and panel pulses 
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TABLE 15. WIDTH CONSTANTS FOR SEVERAL RADIANCE CALIBRATIONS. 
(Means for currents of 4.5 -6.5 amps, 3 bandwidths only) 


Detector 

Bandwidths 

(pm) 

CALIB.4-8 

CALIB.9 

Comment 

0.41 - 0.48 

16.2830 ± 0.1975 
(1.21%) 

17.0864 ± 0.2040 
(1.94%) 

Note greater width constants 
in all cases for CALIB. No. 9 
(across -the -line change in 
pulse shape) 

0.55 - 0.60 

16.5530 ± 0.0885 
(0.535%) 

16.7218 ± 0.048 
(0.287%) 


0.67 - 0.94 

16.6314 ± 0.1435 
(0.863%) 

16.8995 ± 0.081 
(0.479%) 
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Before the lamp or panel pulses are used, the mean level of some specified dark area 
from the previous scanline is subtracted from each reselm level. This removes the effect of 
a varying dark level on pulse levels and allows comparisons between runs. It is assumed that 
the dark level does not vary significantly between any two consecutive scanlines. This as- 
sumption has been valid for all the data sots examined to date. 

The mean peak level of the lamp pulse and panel pulse is determined by selecting a few 
reselms at the top of each in every requested scanline. Then the sets of lamp pulse and panel 
pulse reselms are averaged. These two means are saved for later use. The reselms for the 
lamp pulse are selected by determining the two reselm numbers of the reselm. levels closest 
to half the maximum height of the pulse. These reselm numbers are averaged and two reselms 
on either side of this middle reselm are selected. The reselms for the panel pulse are se- 
lected in a similar manner, except that the two reselms nearest 80% of the maximum height 
of the pulse are used. (This 80% value was used instead of the 50% value because the panel 
pulse is flatter than the lamp pulse and does not drop below 50% of its maximum height.) 

The mean integral levels of the lamp and panel pulses are calculated by performing an 
integration (using Simpson's 1/3 Rule) on the lamp pulse and panel pulse in each scanline 
after normalization to dark level. The lamp integral is then divided by a unique pulse width 
(in reselms) for each spectral channel. The number that results is no longer an area measure- 
ment (in quasi-volt -reselms); instead, it is a level measurement (in quasi-volts). The panel 
pulse is subjected to a similar process. All the adjusted integrals for the lamp pulse and 
panel pulse are averaged to obtain their respective means. The standard deviations of the two 
pulses from their respective means are then calculated. The mean integral ratio is calculated 
by dividing the adjusted mean lamp integral by the adjusted mean panel integral and then mul- 
tiplying by the appropriate factor for each spectral channel. The program then goes back and 
calculates the integral ratio for each scanline and determines the standard deviation from the 
mean integral ratio for the data set. 

The program now performs pulsewidth calculations. The normalized lamp integral is ob- 
tained by multiplying the adjusted mean lamp integral by the pulsewidth constant used to adjust 
it originally. This returns the level measurement to an area measurement. This normalized 
lamp -pulse integral is then divided by the mean lamp level determined earlier. In this man- 
ner, the area measurement is converted to an equivalent -square -wave width (in reselms). 

The adjusted mean panel integral is transformed in a similar manner. 

The advantages of performing pulse -width calculations by integration instead of by the 
mean level method are summarized below: 

(1) The integration process usually results in lower standard deviations than those 
achieved with mean levels. 
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(2) The integration process is faster and simpler than determining 50% or 80% points, 
selecting reselms, and calculating the mean level. 

(3) Changes in calculated pulsewidth provide a way to determine if changes had occurred 
in the shape of the pulse and how drastic they might be (assuming the specified 
pulsewidth is a valid representation of the expected pulsewidth). 

A detailed description of the program LAMRAD. is included in Appendix D. 

4.4 SUMMARY OF CALIBRATION DATA 

All eleven radiance calibration sets were digitized and processed by LAMRAD. Although 
the volume of data generated is too large to be included in this report, a sample of the infor- 
mation obtained is listed in Table 16 for channel CIO (0.46-0.49 4m). As can be seen, the 

calibration constants K~ A (\) have undergone two major, discrete changes associated with 
CAL 

modifications to the system. These modifications include new amplifiers and lamp current 
controls (between CALIB. Nos. 1 and 2) and a different wattage lamp (between CALIB. Nos. 8 
and 9). Even without these modifications, variations of up to 50% in K CAL W are observed. 

It is apparent that more frequent radiance calibrations are needed if the changes observed 
between measurements are to be lowered. 

The NEALs calculated for this channel were obtained using the panel noise and the follow- 
ing relation: 

L x 2a 

NEAL = V1 L 
V L 

where L = radiance at lamp current i 

ct t = standard deviation at panel signal in bins (includes factor of 2 because of smooth - 

ing over 4 lines) 

V' = panel voltage in bins 
Li 

From these values, it is seen that the high frequency noise measured on the reference 
panel also varies somewhat from one calibration date to the next. However, there is no ap- 
parent relationship between the changes in the value of and noise. In fact, two of the 

least noisy signals (CALIB. Nos. 9 and 10) show the greatest change in 

One possible explanation for the change observed in the evaluation of is the 

aging of the reference lamps. A change in spectral character, which has been noted in aging 
quartz -iodine lamps, could indeed cause the differences in calibration constants. To measure 
this, the constants K CAL (M from several calibration data sets were normalized to the value 
at channel C2 (0.67 -0.94 /xm) and the results plotted (see Fig. 9), As can be seen, with the 
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TABLE 16. SUMMARY OF RADIANCE CALIBRATION INFORMATION FOR CHANNEL CIO (X = 0.46-0.49 pm, 10% points) 


CALIB. NO. 

Date 

Mean 

Dark Level 
(bins) 

Mean 

Panel Level 
(bins) 

K Ca! (X) ® 6 - 0A 
(pW-cm-2-ster-l-pm' 1 ) 

% Deviation 

NEAL 

(uW-cm- 2 -ster-l-4m-l ) 

Lamp Width 
(reselms) 

Panel Width 
(reselms) 

1 

20 August 71 

147.84 

11.87 

12,945 

3.4 

54.9 

13.87 

7.82 

2 

6 June 72 

147.66 

26.62 

5071 

2.1 

37.2 

16.50 

7.92 

3 

7 June 72 

148.73 

24.44 

5116 

2.1 

61.2 

16.32 

7.98 

4 

16 June 72 

153.44 

22.64 

5572 

1.9 

33.4 

16.59 

7.95 

5 

12 July 72 

147.99 

16.43 

6867 

2.5 

53.2 


7.40 

6 

No 

CIO 

values 






7 

No 

CIO 

values 






8 

26 January 73 

161.33 

18.61 

4880 

3.0 

63.1 

16.73 

7.66 

9* 

3 April 73 

148.25 

26.58 

2982 

1-6 

25.0 

16.83 

8.09 

10* 

5 June 73 

147.49 

27.54 - 

2299 

2.6 

35.6 

7.68 r 

4- 

7.76 

11* 

6 August 73 

149.05 

19.33 

2111 

2.4 

54.0 

8.10 * 

8.08 


*45 -Watt reference lamp. 

^ Due to slightly different A to O procedure. 
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FIGURE 9. RELATIVE SPECTRAL RESPONSE OF AGED LAMPS. 
Lamps set at 6.0 amps. 
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exception of calibration No* 1, the spectral character of the system has remained nominally 
constant. However, changes of 10% in some spectral bands are present, noticeably in chan- 
nels CIO (0.46-0.49 jim) and C9 (0.48-0.52 iim). Calibration No. 1 does not bear any resem- 
blance to the rest of the data except in the red region, even though the lamp was identical to 
those used in later runs. Since this data set was obtained with the old M5 scanner electronics 
(tlie new electronics were installed six months after the scanner itself), it is possible that the 
observed curve is attributable to the different spectral throughputs rather than changes in 
the lamp. 

In still another evaluation of spectral system response, data from the panel measure- 
ments for all eleven calibration runs were used to calculate a relative response factor given 
by the following; 

L x X X G 

Response, factor = -^r= 

panel, X 

where L = radiance on panel at current i and wavelength X 

X 

A = bandpass of spectral channel 
x 

G = electronic gain 

Vpanel X ~ volta £ e °bt aine d from panel 

A list of these values is given in Table 17, In almost all the solar -reflective channels, 
the scanner spectral response has varied less than 10% during the course of 1 and 1/2 years 
of radiance calibration. The two exceptions in the solar -reflective region are channel C12 
(0.41-0,48 /im), which shows a variation of 22%, and channel C8 (0.50-0.54 fira), which has a 
16% variation. Barring one strange value which occurred during calibration No. 5, channel 
C12 would have only a 7% variation. CALIB. No. 1 is excluded from this analysis as it ob- 
viously exhibits characteristics quite different from the rest. 

The variations observed in the near-IR region tend to be larger than in the shorter wave- 
length region, averaging approximately 15 to 20%, Considering the nature of these detectors 
and their associated preamps, however, this larger deviation is not surprising. In general, 
the response factor shows that the overall system has not significantly changed its response 
in almost two years of operation* 
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TABLE 17. RELATIVE SPECTRAL RESPONSE FACTOR 


Response = 


L X 4 XG . 

A A nun 

V 

panel A 


Calibration No. 


* 


Channel 

1 

2 

012 

— 

0.086 

CIO 

0.169 

0.046 

C9 

0.163 

0.028 

C8 

0,104 

0.061 

C7 

0.063 

0.022 

C6 

— 

0.042 

C5 

0.215 

0.094 

C4 

0.498 

— 

C2 

1.0 

1.0 

D3 

346.7 

0.0160 

D2 

131.2 

0.0076 

D1 

— 

0.0055 

B3 

— 

— 

B2 

— 

— 


r 


0.088 

0.08B ! 

| 0.035 1 

0.051 

0.053 

0.050 

0.026 

0.025 

0.026 

0.060 

0.060 

0.070 

0.022 

0.022 

0,022 

0.042 

0,042 

0.043 

0.097 

0.095 

0.097 

— 

0.143 

0.141 

1.0 

1.0 

1.0 

0.0163 

— 

0.0147 

— 

0.0075 

0.0051 

0.0081 

0.0053 

— 


6 

7 

8 

0.089 

0.102 

0.098 

— 

— 

0.053 

0.025 

0.030 

0.028 

0.078 

0.095 

0.085 

0.021 

0.023 

0.022 

0.041 

0.050 

0.047 

0.094 

0.102 

0.100 

0.137 

0.151 

0.145 

1.0 

1.0 

1.0 

0.0157 

— 

— 

0.0061 

_ 



— 

0.058 

0.055 

— 

0.0050 

0.0073 


9 

10 

11 

0.095 

0.096 

0.109 

0.044 

0.047 

0.053 

0.025 

0.022 

0.023 

0.077 

0.084 

0.091 

0.021 

0.020 

0.020 

0.045 

0.043 

0.056 

— 

0.096 

0.097 

0.143 

0,148 

— 

1.0 

1.0 

1.0 

1 0.055 

t 0.0183 


0.0056 0.0066 

0.0057 

0.0062 0.0071 

0.0065 


Deviation 

(percent) 

Mean (2-11) 

22 

0.0886 (0.0191) 

7 

0.0954 (0.0068) 

7 

0.0496 (0.0033) 

9 

0.0258 (0.0023) 

16 

0.0761 (0.0123) 

4 

0.0215 (0.0009) 

10 

0.0451 (0.0045) 

3 

0,0969 (0.0025) 

3 

0,1440 (0.0042) 


1.0 

78 

0.0332 (0.0259) 

7 

0.0162 (0.0012) 

16 

0.0064 (0.0010) 

14 

0.0064 (0.0009) 

3 

0.0565 (0.0015) 

19 

0.0062 (0,0012) 


3 4 5 


*Note: Measurements made with M-5 scanner electronics. 
^Note: Appear to b£ erroneous. 
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5 

MONITORING OF THERMAL DATA QUALITY 

Users of thermal data from a scanner such as the M7 usually are interested more in the 
accuracy (deviation from true value) of the temperature measurements than in some of the 
other performance parameters such as NEAT or stability. This is because many investigators 
(especially at this stage in technique development) are attempting to extend limited ground- 
truth information on temperatures to much larger areas by the use of imagery obtained by 
aircraft or satellite. When measurement differences occur between these sources of data, 
not only does the extension of temperatures become difficult but the researcher must do a 
considerable amount of investigation to identify (or justify) the causes of the differences. 

This is not to say that NEAT and long-term stability are not important parameters. 
Obviously, a system that cannot resolve, say, a 5°C temperature difference would not be suit- 
able for most applications involving small temperature ranges. Also, knowledge of NEAT is 
a good indicator of system performance in the thermal region. In most cases, though, the re- 
searcher will tolerate poorer temperature resolution (as long as it is not excessively poor) to 
gain exact temperature calibration. 

It is not the objective of this study, however, to improve the method of thermal calibra- 
tion. That is, while some new techniques concerning thermal data quality were analyzed, 
none specifically involved the method of calibrating the reference plates themselves. This 
study was addressed to the problems of maintaining some degree of confidence for a given 
thermal calibration and determining how often these thermal calibrations are needed. The 
remainder of this section presents the results of this investigation. 

5.1 PROGRAM THERML. 

As an aid in analyzing thermal data, the program THERML. was written. This program 
calculates a set of statistics for the hot reference plate, cold reference plate, and ambient 
plate similar to that which CALIB. calculates for the solar reflective data channels. The main 
difference is that no data-base tape forgathering a history of thermal response is available, nor 
can more than one channel of thermal data be processed at one time. Since the hot plate and 
cold plate are represented by fairly flat pulses, calculations to find the mid-point need not be 
performed on the pulses. Any width of flat area on the plates can be selected by specifying 
the beginning and ending reselm numbers of each of the plates. A default case does exist, 
however, for routine data reduction. 

To implement this program, the user simply specifies the location of the cold plate, the 
hot plate, and ambient area within each scanline of a scene. The program then calculates the 
average level and standard deviation of the cold plate. The same thing is done for the hot 
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plate and the ambient area, except that the mean level and standard deviation are relative to 
the mean level of the cold plate from the previous scanline. This assumes that the mean level 
of the cold plate does not change significantly between any two scanlines of a scene. (This 
proved to be a valid assumption throughout all the data sets we analyzed.) These statistics 
are calculated for as many data sets as required. 

After all the statistics have been calculated for the requested data sets, the program 
plots the means of the cold plate, the hot plate, and the ambient areas versus the data set 
numbers. The same thing is done for the standard deviations. These plots are followed by a 
table of all the statistics calculated. Together, these plots and the table provide the user 
with comparisons between data sets so that trends and discontinuities can be determined. 

Once the user has found that a trend has developed or a discontinuity has occurred, the causes 
can be analyzed more thoroughly via other techniques available to him. Detailed documenta- 
tion of Program THERML. is given in Appendix E. 

5.2 PROCEDURE FOR CHECKING THERMAL DATA QUALITY 

The M7 scanner has two thermo -electrically heated (or cooled)-plates which serve as 
temperature references for the determination of apparent temperatures at the entrance aper- 
ture of the scanner. These plates are checked periodically to ensure that the plate tempera- 
ture (apparent temperature) versus thermistor resistance (thermistors are embedded in the 
plates) has not changed. Though it has been determined that in general these plates do not 
change their radiative characteristics much from one calibration to the next, individual plate 
failure has often occurred. Once repaired, these plates exhibit temperature characteristics 
identical to those observed before failure. It is primarily during these periods of plate 
failure, then, that uncertainties exist concerning plate calibration. 

Plate failures are caused by changes in operating characteristics of the thermoelectric 
devices used to control plate temperature. Although no solution has yet been worked out to 
prevent plate failure, a technique for monitoring plate performance has been found. Shown in 
Table IB is a list of calculated and actual ambient surface temperatures located inside the 
scanner housing for a number of missions during 1972 and 1973. The surface being viewed 
for this experiment was a plate (identical to the reference plates) installed in an unused por- 
tion of the scanner housing during the summer of 1972. Although this plate was unpowered by 
thermoelectric devices and left to drift with ambient temperatures, a precision thermistor 
was embedded in its center to permit temperature monitoring. 

The data in Table 18 provide comparison between the ambient surface temperature of 


this plate as determined by the thermistor and the apparent temperature calculated by using 
the hot and cold reference plates. Although some difference can be expected since surface 




TABLE 18. CALCULATED PLATE TEMPERATURE COMPARED TO MEASURED 

PLATE TEMPERATURE 



Calculated 

Actual 



Ambient Temperature 

Surface Temperature 

Temperature 

Date 

(°C) 

(°C) 

Difference 

23 July 72 

14.20 

14.13 

-0.07 

24 July 72 

18.16 

16.38 

-1.78 

18 August 72 

19.77 

19.77 

0 

25 August 7 2 

25.80 

27.09 

+1.29 

28 August 72 

20.13 

20.59 

+0.46 

29 August 72 

23.39 

23.97 

+0.58 

29 August 72 

20.54 

21.07 

+0.53 

30 August 72 

24.71 

25.42 

+0.71 

5 September 72 

23.50 

23.94 

+0.44 

14 September 72 

14.57 

14.77 

+0.20 

17 October 72 

7.57 

6.85 

-0.72 

16 November 72 

18.42 

18.75 

+0.33 

17 November 72 

19.33 

20.51 

+1.18 

2 January 73 

6.54 

5.25 

-1.29 

10 January 73 

-2.37 

-5 .35 

-2.98 

12 January 73* 

10.02 

1.96 

-8.06 

25 January 73 

16.05 

15.67 

-0.38 


♦Reference plate failure 12 January 73. 
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temperatures are being compared to apparent temperatures, a major discrepancy would in- 
dicate a possible reference plate malfunction- As noted in the third column (containing tem-. 
perature differences), a significant error was observed on 12 January 1973. This coincides 
exactly with a scanner operator’s observation of a plate failure on that date. As noted for the 
three previous dates, there is strong evidence that some degree of malfunction, which went 
unobserved by scanner operators, was present even then. 

The data shown in Table 18, however, do not answer all the questions about plate per- 
formance- Although the principle is sound, only one failure was observed (for analysis pur- 
poses) during the period covered in the table. Other failures were observed during the spring 
and summer of 1973, but the ambient plate used in the test had by then been removed to make 
room for a second reference lamp for the solar reflective region. Also, the dates indicated 
in the table do not include all flights made with the M7scanner — the need for standardization 
tests (including the thermal regions) on every flight had not yet been demonstrated. 

Finally, two other entries (24 July and 25 August 1972) in the table indicate a difference 
significantly in excess of 1°C between measurements, but they have not been associated with 
a plate failure. The reason behind the larger discrepancies which occurred on these two 
dates is unknown, but it is obvious that some quantitative measure of maximum difference 
must be generated before the technique will become usable. 

During the winter of 1973, a new ambient plate w T ill be constructed for use during stan- 
dardization data runs. With these data runs being made every flight, some quantitative judg- 
ment of thermal reference plate operating performance (and, indirectly, of temperature cali- 
bration) should be available for each mission. 

5.3 NOISE EQUIVALENT AT (NEAT) 

In addition to monitoring thermal plate performance (which determines temperature 
accuracy), it is also very informative to monitor the NEAT of the thermal channel. By doing 
this, a check of the system’s temperature resolution is possible. Before going on, however, 
we should be careful to distinguish the difference between the NEAT of the system and the 
NEAT of the infrared detector. This distinction is important because the NEAT observed by 
an investigator who is digitally processing recorded data includes noise from the optics, de- 
tector, amplifier, and tape-recorder components of the system. Hence, a particularly poor 
NEAT could be the result of a predominance of noise from any one of these components. And, 
in fact, tape-recorded thermal data from the M7 is usually limited by the noise of the tape re- 
corder, as will be shown later. 

The infrared detector itself, however, may have a significantly better NEAT than that ob- 
served from processed data. In the ideal system, the noise of the system would be detector- 

limited; hence, ultimate temperature sensitivity depends on the detector. Of the four 
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components mentioned above, the detector is the least rugged and most likely to fail. Thus, 
it seems reasonable to monitor the NEAT of the detector in conjunction with total system 
NEAT to ensure proper system performance. This study, however, will concern itself only 
with the system NEAT measurements since most investigators are interested in temperature 
resolution primarily from the data processing standpoint — i.e., what results can be obtained 
from the recorded signals. 

As mentioned above, the tape-recorded thermal data from the M7 scanner is usually lim- 
ited by the noise of the tape recorder. The Mincom PC-500 tape machine used in the M7 system 
has a dynamic range of 39.5dB (peak-to-peak), which means that 89 temperature levels are 
available between the minimum and maximum temperatures recorded (this assumes that the 
full range of the tape machine is being utilized, which is seldom the case). For a typical tem- 
perature range of 10°C, the system NEAT would be 0.11°C (assuming a tape-recorder-limited 
system). It should be obvious, here, that the NEAT in this system is a function of the terrain 
temperature range encountered — the larger the range, the larger the NEAT. Since the nominal 
NEAT of the HgCdTe infrared detector used in the M7 is 0.03°C (which is independent of AT), 
it is apparent that, in most cases, the noise of the tape recorder is the limiting factor in tem- 
perature resolution. 

In order to quantitatively assess the noise performance of the thermal channel, data ob- 
tained from the standardization runs in 1972 and 1973 were analyzed for system NEAT. Table 
19 is a list of the NEAT s calculated from that data. The theoretical tape recorder limit is 
plotted (assuming 39.5dB dynamic range) along with data generated from the noise levels ob- 
tained from the cold plate, hot plate, and ambient plate. These NEAT s were calculated ac- 
cording to the relationship 


NEAT = 


AT X 2(r 
AV 


where AT = temperature difference between reference plates 

o - rms noise calculated from program THERML. — the factor of 2 appears because 
the data were averaged over 4 lines 
AV - voltage difference between reference plates 

In viewing these data, however, certain things must be kept in mind. As mentioned above, 
the theoretical recorder limit applies to the full dynamic range of the tape machine. In prac- 
tice, this full range is not used. Hence, the recorder limit may be lower than the theoretical 
limit by as much as 33%. Secondly, the cold-plate noise figure is based on absolute voltages 
and includes a low-frequency drift inherent in the preamplifiers for the HgCdTe detector. 
Since the hot plate and ambient plate are measured relative to the cold plate, this drift does 
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TABLE 19. SUMMARY OF NEAT VALUES OVER TIME 


NEAT 




Recorder 

Cold 

Hot 



Date 

AT 

Limit 

Plate 

Plate 

Ambient 

Comments 

14 May 72 

7.9 

0.09 

0.35 

1.09 

0.20 


5 May 72 

6.5 

0.07 

0.69 

0.71 

0.48 


11 May 72 

9.9 

0.11 

0.48 

0.65 

0.36 


11 May 62 

9.9 

0.11 

2.12 

3.20 

1.71 


12 May 72 

9.8 

0.11 

0.89 

0.49 

0.76 


12 May 72 

10.8 

0.12 

0.53 

0.45 

2,12 


19 May 72 

6.5 

0.07 

2.61 

1.07 

0.96 


25 May 72 

6.5 

0.07 

0.56 

0.49 

0.37 


25 May 72 

6.5 

0.07 

0.41 

0.54 

0.39 


5 June 72 

9.7 

0.11 

0.51 

0.53 

0.39 


23 July 72 

10.3 

0,11 

1.07 

0.77 

0.61 


24 July 72 

7.3 

0.08 

1.70 

0.71 

0.64 


18 August 72 

9.1 

0.10 

0.72 

2.12 

0.59 


25 August 72 

10.3 

0.12 

6.47 

2.61 

2,22 


28 August 72 

10.3 

0.12 

0.98 

0.65 

0.59 


29 August 72 

9.1 

0.10 

0.70 

0.20 

0.34 


29 August 72 

10,3 

0,12 

0.48 

0.41 

0.23 


30 August 72 

10,3 

0.12 

0.65 

0.24 

0.23 


5 September 72 

9,1 

0.10 

0.81 

0.30 

0.32 


14 September 72 

15,7 

0,10 

0.61 

0.83 

0,32 


17 October 72 

15.7 

0.18 

0.95 

1.87 

0.34 


16 November 72 

5.6 

0.06 

0.35 

0.35 

6.65 


17 November 72 

5.6 

0.06 

3.33 

0.45 

13.00 


2 January 73 

15.7 

0.18 

0.94 

3.94 

0.67 


10 January 73 

15,7 

0.18 

0.87 

4.81 

12.91 


12 January 73 

15.7 

0.18 

1.25 

2.57 

1.08 

Plate failure on this date 

25 January 73 

15.7 

0.18 

0.53 

1.97 

7.71 


22 June 73 

15,7 

0.18 

0.52 

1.52 

— 


25 June 73 

15,7 

0.18 

0.44 

1,62 

— 


3 July 73 

15.7 

0.18 

1.63 

1.73 

— 


5 July 73 

15.7 

0.18 

0,63 

2.58 

— 


6 July 73 

15.7 

0,18 

0,96 

1.22 

— 


7 July 73 

15.7 

0.18 

0.58 

1.43 

— 


4 August 73 

15.7 

0.18 

0.70 

2.07 

— 


5 August 73 

15.7 

0.18 

4,04 

2.43 

— 


10 August 73 

15.7 

0.18 

0.43 

2.82 

— 


12 August 73 

15.7 

0.18 

0.52 

1,47 

— 

Plate failure on 8/17/73 
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not affect them. Finally, the ambient plate is unpowered (i.e., it has no thermoelectric de- 
vices); so, any noise associated with thermoelectric devices would appear only in the cold and 
hot plate values. 

In looking at the hot plate NEAT s, it is quite easy to see which missions appear to have 
excessively noisy temperature resolutions. An arbitrary limit is any value over three times 
the recorder limit. Equipment performance on most missions prior to January 1973 appears 
to be nominal. But note the NEAT measurements during January just prior to a reference 
plate failure; obviously the signals were becoming consistently noisy. Also interesting is the 
fact that all NEAT s calculated from the hot plate during the period June to August 1973 appear 
to be high (no ambient values are available because the ambient plate had been removed). On 
August 17, soon after the last mission listed in the table, a reference plate failure occurred 
and a degraded HgCdTe operating performance was noted. It would appear that one or both of 
these problems had been affecting the thermal channel performance for several weeks. How- 
ever, since the noise on the cold plate during this period is generally acceptable, more than 
likely it was the reference hot plate that caused the excessively high noise. 

Given more frequent standardization runs and the availability of a new ambient plate for 
additional data, these NEAT calculations should provide a quantitative evaluation of every 
mission, whether or not the thermal channel is performing within the ideal limit. 

5.4 HIGH AND LOW FREQUENCY THERMAL NOISE 

Until now, statistical comparisons of thermal data have been made between missions. One 
of the questions to be answered in this study was: What frequency variation in thermal ref- 
erence signals can be expected during the course of a 12-minute run? In other words: Is the 
dominant noise in the thermal channel mainly attributable to high-frequency (<16 msec) or low- 
frequency (<12 minute) components? The frequency distribution should tell us (1) whether static 
reference-plate calibration of the data is adequate or whether dynamic calibration is needed; 
and (2) if dynamic calibration is dictated, how often must it be done. 

Using the program DGNSTC., we analyzed for noise a 12-minute segment of data obtained 
on 7 April 1973 over Lower New York Bay. To obtain a low-frequency component, the mean 
and standard deviations of the hot and cold plate were calculated for all lines in the entire 12- 
minute segment: 

Cold Plate: mean plate voltage = 185.07 

standard deviation (between lines) = 92.41 (50.0%) 

Hot Plate: mean plate voltage = 33.27 

standard deviation (between lines) = 14.45 (43.5%) 
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Although the high deviation in cold-plate value can be attributed to drift (low-frequency 
noise) in the infrared detector preamplifiers, the high deviation in the hot plate value is in- 
explicable since it has been measured relative to the cold plate. In any case, because these 
significant deviations occur during the course of a 12 -minute data collection run, it is apparent 
that static dark level corrections are out of the question. 

To measure high-frequency thermal noise, the program DGNSTC.also calculated the 
standard deviation of each line. This would be the noise encountered if dynamic thermal cali- 
bration were done on a line-by-line basis. The calculations obtained are as follows: 

Cold Plate: mean plate voltage - 185.07 

average standard deviation per line = 3.49 (1.9%) 

Hot Plate: mean plate voltage = 33.27 

average standard deviation per line - 2.52 (7.6%) 

The average cold -plate standard deviation of each line over the entire 12-minute run is only 
1.9% of the mean value as compared to a 50.0% standard deviation of all lines. Line-by-line 
dynamic temperature calibration is obviously needed since it provides a significantly lower 
error than that observed statistically from all lines. 
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6 

CONCLUSIONS AND RECOMMENDATIONS 

It has been demonstrated thatvaluable information about scanner performance canbe gained 
from the calibration data ("standardization” runs) obtained during data acquisition missions. In 
the solar reflective region, noise information (as a function of frequency) obtained digitally 
from "standardization" runs can indicate possible amplifier, detector, or recorder problems. 

To permit up-to-date monitoring of these noise figures, it is recommended that standardization 
data be acquired at least once every flight. Rapid processing of this information is essential 
to its optimal use. 

A region of darkness viewed by the rotating scan mirror in one 360-degree swung was 
obtained for baseline correction. We recommend that all data processors of M7 multispectral 
data use this region, which is located beyond the lamp reference but before the sky reference. 

A new' technique for obtaining values for the reference lamp pulse was examined. This 
technique, which uses integration instead of peak pulse value, is advantageous in that more 
points are available for averaging with a resulting lower noise figure. In addition, the inte- 
gration is relatively insensitive to jitter in lamp pulse position. For optimal use of the refer- 
ence lamp pulse information, this technique should be applied to all software associated with 
the extraction of radiance information. 

Analyses of all past laboratory radiance calibrations on the M7 scanner have revealed 
large differences between "adjacent" measurements (currently scheduled approximately every 
two months). Such differences (as large as 50%) appear even in the absence of a physical 
alteration in the system — such as a new reference lamp, for example. If smaller changes 
between radiance calibration measurements are desired, these measurements will have to be 
made more frequently, most likely every 2 to 3 weeks. 

Finally, we tested methods for evaluating both thermal reference plate performance and 
overall detector, amplifier, and tape recorder performance in the thermal region. Plate failure 
often begins to occur long before it is observed by the scanner operator. But it is possible, by 
using an ambient plate as a third reference, to determine w'hether thermal reference plates 
are failing. Thus, inclusion of an ambient reference in the standardization measurements 
made every flight would be desirable in order to permit reference plate evaluation as often as 
possible. An overall thermal region performance indicator has been demonstrated by the use 
of NEAT calculations. These permit a theoretical tape recorder limit to be compared to actual 
measured values. Since in most cases the system is tape-recorder-limited, any NE AT f s 
larger than the theoretical limit can indicate some type of excessive noise in the system. 

In summary, it has been shown that systematic recording and digital analysis of M7 scanner 
calibration signals can be valuable in evaluating both scanner performance and data quality. 
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As we acquire better knowledge through continued analysis of these calibration signals, we can 
more easily satisfy future data requirements to meet new user goals. 
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Appendix A 
PROGRAM CALIB. 

The purpose of this program is to calculate a set of statistics concerning the reference 
sources of the M7 scanner. These statistics may be calculated for as many flight missions as 
desired (the only limitation is the physical amount of core storage in the computer). After 
calculating the statistics for the requested files and channels, the program can merge them 
with those on a special tape data base (initially produced by this program), produce plots and 
tables, and update a tape data base for each of 24 possible spectral channels. In addition to 
creating the initial tape data base, the program can be requested to simply produce plots and 
tables from a given tape data base previously generated. 

The statistics calculated by this program are the following: 

(1) The mean of the dark level 

(2) The mean of the normalized lamp pulse level 

(3) The mean of the normalized level of the middle reselm in the lamp pulse 

(4) The mean of the reselm number of the middle reselm in the lamp pulse 

(5) The mean of the integral of the normalized solar pulse 

(6) The standard deviation from each of the previous five calculated means 

All of these statistics are normalized to 100X gain except the solar statistics (normalized 
to 10X gain). The gain in each channel is read as input. The lamp pulse points used in the 
calculation of its mean are determined as follows: 
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(1) Maximum height is determined [MAX] 

(2) Reselms closest to MAX/ 2 are found [X & Y] 

(3) The mean of X & Y is found (X + Y)/2 

(4) The bounds used in calculations are TOP/2 reselms on either side of (X + Y)/2 

Normalization is achieved by subtracting the average dark level from each point of the 
lamp pulse and solar pulse before they are used. The mean dark level is taken from the pre- 
vious line. Consequently, in order to start this chain, the lamp pulse and solar pulse are 
ignored in the first line of any file so that the mean dark level can be calculated. We assume 
that the dark level does not change significantly from line to line. 

CAUTION: For the program to operate correctly, the lamp pulse and the solar pulse must be 
positive -going. 

STEP (1 )- PERFORMED 1ST TIME FUNCTION CALLED OR WHENEVER M0DE1 - $START$ 

The maximum number of channels to be referenced and the maximum number 
of files to be referenced are initialized. The beginning and ending reselms for 
the dark area, lamp pulse, and solar pulse are set and other pertinent variables 
are set. 

Data cards (if any) are read that may change the default settings for the 
type of data base (TAPE = $OLD$, $NEW$, OR $PRINT$), the maximum number of channels, 
the maximum number of files, the maximum number of missions to be plotted (effective 
only if TAPE=$PRINT$) , the tape numbers of the tapes to be used for the tape data base, 
and the number of reselms at the top of the lamp pulse to be used in calculation of 
the lamp statistics. These variables are now checked for reasonableness. The storage 
for the main-in-core data base is then allocated, and the totalling arrays are 
zeroed. Control returns to POINT. 

STEP (2)- PERFORMED BEFORE EACH FILE IS PROCESSED 

The point counters for the lamp and solar pulses, for the data, and for the 
lines are zeroed. A flag is set so that the 1st line of the file will not be 
processed by the lamp or solar pulse routines. 

If this is the first time through the file, data cards are read to set the 
data (in month, day, year form), the mission number, and any special comment. The be- 
ginning and ending reselms for the dark area, lamp pulse, solar pulse (if different 
from the default settings are then specified. At this point, as many cards are read 
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in as there are channels on the input tape. Each card specifies the multiplexor chan- 
nel number, the spectral channel it represents, the aircraft channel it represents, and 
the gain in that channel. A title describing the spectral channel is printed 
along with each card. The spectral channel number and the multiplexor channel 
number are used to set up an array to relate the multiplexor channels to the 
spectral channels. [If a particular multiplexor channel is not to be used, simply 
assign zero as its spectral channel #.] This array is then checked for reasonable- 
ness. Likewise, the dark, lamp, & solar area boundaries are checked as well as the 
date and the mission number. No more than 100 reselms are allowed for either 

the lamp pulse or the solar pulse. In addition, there must be- some dark area 
specified. Control returns to POINT. If this is the second pass through the file, 

nothing is done and control returns to POINT. 

STEP(3) - PERFORMED BEFORE EACH POINT OF A SCAN LINE IS PROCESSED. 

No action taken. 

STEP (4) - PERFORMED AFTER EACH LINE IN FILE IS PROCESSED 

If this is the first line processed (FLAG=0) , the line counter is incremented, 
and control returns to POINT. . After incrementing the line counter, the following is 

performed on each nonzeroed spectral channel. 

1. The maximum point in the lamp pulse is found. 

2. The reselm location of the two points that are 1/2 the maximum or just 
greater than 1/2 the maximum are found (these will correspond to the boundaries 
of the lamp pulse if only the top of the pulse is specified), 

3. The reselm half-way between the preceding two points is found. New 
boundaries are set on either side of this middle according to the number of reselms 
requested for use at the top of the lamp pulse (LAMTOP) [Right & Left]. 

4. The new subset of reselms is checked to make sure they are within the 
original boundaries for the lamp pulse,, (if not, the calculated bounds are 
adjusted accordingly). 

The data point counter is incremented by right-left +1. 

5a. If this is pass 1, the values of the reselm subset of the lamp pulse are 
accumulated, the value of the middle reselm is accumulated , and the reselm number 
of the middle reselm is accumulated. 
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5b. If this is pass 2 , the squared deviations are accumulated, for the 
reselm subset, for the middle reselm, and for the reselm number. 

6. The number of points in the solar pulse is checked to make sure that it 
is greater than 5 and is odd. 

7. A Simpson integration is then performed on the points in the solar pulse 
and its result accumulated if pass one. Otherwise the squared deviation is 
accumulated. 

After all the multiplexor channels requested are processed, the point counters 
for the lamp pulse and the solar pulse are zeroed. 


STEPS - PERFORMED AFTER AN ENTIRE FILE HAS BEEN PROCESSED 
If this is pass 1; 

1. The file counter is incremented, 

2. The means are calculated for each channel requested. 

3. The information is stored in the main data base after being normalized 
to gain 100* 

4. "PASS" is set to 2 and the totalling arrays are zeroed. 

If this is pass 2: 

1. The standard deviations are calculated for each channel requested and 
stored in the main data base after being normalized to gain 100. 

2. "PASS" is set to 1. 

3. The date, mission number, and special comments are saved. 

4. The totalling arrays are zeroed, 

STEP6 - PERFORMED WHENEVER MO DEI - $FINAL$ 

Whenever "TAPE = $0LD$ or $NEW$; 

1. A tag sort is performed on the data array. 

2. The tapes specified by XBASE 1 & TBASE2 are mounted on tape drives 
2 and 3. 

a. If "TAPE" = $NEW$, tape unit 2 is made the output unit and its 
generation number is set to 1. 

b. If "TAPE" = $0LD$, tape unit 2 is checked for a generation number 
of 1. If so, tape unit 2 is set for input and tape unit 3 is set for 


84 



FORMERLY WILLOW RUN LABORATORIES, THE UNIVERSITY OF MICHIGAN 

output. If the generation number is not 1, the tape unit with the greater 
generation number is set for input, and the tape unit with the lower 
generation number is set for output. 

3. The bin number and generation number are written on the output tape. 

4. A call is made to PL0T1. to set up the grid. The minimum and maximum 
mission numbers to appear on the plot are determined. (The maximum mission 
number of the last file processed. The minimum mission number is 1 or the 
maximum mission $ minus 61- whichever is larger J 

For each spectral channel the following are performed: 

1. If "TAPE” = $PRINT$ , the number of missions in the spectral band is read 
from the tape data base. If "TAPE” - $NEW$ , this number is always zero. 

2* The storage necessary for the X and Y coordinate matrices is allocated. 

3. The calculated statistics and the statistics from the tape data base are 
merged (ascending order based on date). 

4. The merged means are then plotted against their corresponding mission 
numbers FOR THE LAMP ONLY. 

5. The merged standard deviations are plotted against their corresponding 
mission numbers. 

If "TAPE 11 = $PRINT$, the tapes are mounted and checked as before. The input 
tape is found also. In this case, it is assumed that there are no calculated 
statistics to merge. Therefore, only the data on the tape data base is plotted 
and tabulated as before. No output tape is generated and the existing tape 
data base remains unchanged. The maximum mission number to be plotted along 
the X axiB is determined by the cards read in during STEP1. 

PTPROC . 

This routine calculates the average dark level for each line. This average 
is used to normalize the data points that are saved in the lamp and solar pulse 
arrays. In addition, the dark level is accumulated over the specified area. If 
m FLAG m = 0 (indicating that this is the 1st line of the file to be processed), 
the data points for the lamp pulse and the solar pulse are not saved. 



85 



Terim 


FORMERLY WILLOW RUN LABORATORIES, THE UNIVERSITY OF MICHIGAN 


SIMPSN, 

This routine performs a Simpson integration of the data points in the solar 
pulse. The formula used is: 

Q 22 = X x + 4X 2 + 2X 3 + 4X 4 +....+ 4X n _ 2 + 2X n _ 1 + X r 

n must be ODD 
n > 5 


INSERT. 

This routine inserts either the calculated 
data tape into the 3-dimension coordinate array 
base. 

STEP (1) CHANNL 13 

INTEGER 

STEP (1) FILMAX 1 

INTEGER 

STEP(l) LAMTOP 5 

INTEGER 


STEP (1) TAPE $OLD$ 

ALPHABETIC 


statistics or the statistics from the 
composing the main in-core data 

The maximum multiplexor channel number being 
referred to by any file. It may be less than 
the actual number of channels available. It 
must not be less than 1 or greater than 13. 
The maximum number of input files to be 
processed , plotted >and tabulated as a set 
(you may process fewer files but not more). 
This variable must not be less than 1. 

The number of data points at the top of the 
lamp pulse to be used In the statistical 
calculations. Fewer than this number may be 
used if the given bounds for the lamp pulse 
are smaller. The variable must not be 
less than 0 or greater than L AMP E- LAMP B+l . 
$OLD$ indicates that a tape data base 
exists and the data in it should be merged 
with the processed data for plotting and 
tabulation, 

$NEW$ indicates that no tape data base 
exists and that a data base should be 
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STEP (1) MISSON 0 

INTEGER 

STEP (1) TBASE1 799 

INTEGER 

STEP (1) TBASE2 800 

INTEGER 


generated from the present data. This new 
tape data base is written on tape TBASEL. 

$ PRINT $ indicates that no data is to be 
processed. The existing current tape data 
base is plotted and tabulated. No tapes are 
altered. MISSON" must also be assigned a 
value. The number assigned to MISSON will 
be the number of the last mission table 
plotted. All higher- numbered missions on 
the tape data base will not be plotted. All 
mission numbers less than MISSION-61 or 
1 (whichever is higher) will not be plotted 
also. In either case all missions will be 
tabulated. Whenever TAPE = $0LD$ or $NEW$- 
the mission number of the last file processed 
will be the last mission number to be plotted. 
Higher-numbered missions will be plotted 
and mission numbers less than MISSON-61 
or 1 (whichever is higher) will not be 
plotted. In either case, all missions will 
be tabulated. 

This is the number of the last mission to be 
plotted and tabulated if TAPE = $PRINT$ . 

This number must not be less than 1. 

The bin number of a tape to be used for 
part of the tape data base. 

The bin number of the other tape of the 
tape data base. 

NOTE: No error checking is done either 

On TBASE1 or TBASE2 before the 
mount commands are issued. 
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STEP (2) FDATE(l). . .FDATE(3) 0 

INTEGER 


STEP (2) MISSON 0 

INTEGER 


STEP (2) COMMNT (1) . . . BLANK 

COMMNT ( 8 ) ALPHANUMERIC 

STEP (2) LAMPB 10 

INTEGER 

STEP (2) LAMPE 60 

INTEGER 


STEP (2) DARKB 185 

INTEGER 

STEP (2) DARKE 225 

INTEGER 

STEP (2) SOLARB 310 

INTEGER 


Any date to be associated with the present 
file being processed* 

PDATE(l) — the month must not be less than 
0 or greater than 12. 

FDATE(2) — the day must not be less than 0 
or greater than 31. 

FDATE(3) — the year must not be less than 0 
or greater than 99. 

The number of the mission to be associated 
with the present file being processed. If 
this is the last file being processed, 
this is also the highest mission number that 
will be plotted. This variable must not be 
less than zero. 

Any 48 -character comment to be printed 
along with the tabulated results for this 
file. 

This is the number of the reselm indicating 
the beginning of the lamp pulse. This number 
must not be less than the beginning reselm 
number of the line. 

This is the number of the reselm indicating 
the ending of the lamp pulse. The variable 
must not be greater than the existing reselm 
number for the line. 

Same as LAMPB for beginning of the dark area. 
Same as LAMPE for ending of the dark area 
Same as LAMPB for beginning of the solar area. 
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STEP (2) SOL ARE 
STEP (2) MPXMUX 

STEP (2) SCHAN 


STEP<2) AIRCHN NO DEFAULT 
INTEGER 

STEP (2) GAIN NO DEFAULT 
INTEGER 


355 

INTEGER 

NO DEFAULT 
INTEGER 

NO DEFAULT 
INTEGER 


Same as LAMPE for ending of the solar area. 

The number of the multiplexor channel. Must 
not be less than one or greater than 13. 

The number of the spectral channel that 
MPXMUX represents. This variable must not 
be less than 0 or greater than 24. When 
SCHAN - 0, the MPXMUX channel will be 
ignored. 


SCHAN 

SPECTRAL BAND 

1 

FI, 

.71-, 73 

2 

El, 

.33-. 38 

3 

B2, 

2. 0-2. 6 

4 

B3, 

1.0-1. 4 

5 

Dl, 

2.0-2. 6 

6 

D2, 

1.5-1. 8 

7 

D3, 

1.0-1. 4 

8 

C2, 

.67-. 94 

9 

C4, 

.62-. 70 

10 

C5, 

.58-. 64 

11 

C6, 

.55-. 60 

12 

C7 , 

.52-. 57 

13 

C8, 

. 50-. 54 

14 

C9, 

.48-. 52 

15 

CIO, 

.46 *.49 

16 

C12 , 

, .41-. 48 

17-24 

Not Assigned 


This is the number of the aircraft channel MPXMUX represents. 
No error checking is performed. 

This is the gain to be associated with this multiplexor 
channel. 
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OUTPUT 

PLOTS 


Plots are made only for those spectral channels with data. The first plot of 
each pair is a plot of the mean LAMP LEVEL, 

Plotting Char. 


D 

L 

M 

S 


mean dark level 

mean lamp level 

mean level of middle of lamp 

mean integral of solar 


The second plot of each pair is a plot of the standard deviations. The same 
plotting characters are used as before. After each pair of plots follows a table of 
DATE date 

MISSION mission # 

DARK MEAN mean of dark level x 100 

LAMP MEAN mean of lamp level x 100 

LAMP MIDDLE MEAN mean level of middle of lamp x 100 


SOLAR MEAN 
DARK DEV 
LAMP DEV 


mean integral of solar pulse x 10 
" *’ tr " " x 100 

n ii „ „ „ x 100 


LAMP MIDDLE DEV \ standard deviations from their respective means X10 


SOLAR DEV 


J 


respective means x 10 


LAMP RESELM MEAN mean // of point of middle of lamp pulse (ADD LAMPB-1 
to get actual reselm ft) x 100 


LAMP RESELM DEV standard deviation of above x 100 
Lamp pulse and solar pulse are normalized by subtracting average dark level 
from each point before processing. 


DATA CARDS SETUP 

1. FILMAX = 3, TAPE = $NEW$, CHANNL = 5* 

2. INBIN = 706, UNIT = 4, FILE = 1, NSA - 1,100,1,1,359,1 * 

3. FDATE(l) = 01,11,72, MISSON = 4 

4. COMMNT (1) = $TH1S IS A SPECIAL COMMENT $* 

5. MPXMUX = 1, SCHAN = 4, AIRCHN = 7, GAIN - 2* 
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6 . 

MPXMUX = 2, SCHAN = 5, AIRCHN = 4 f 

GAIN =4* 

7. 

" - 3 M = 0 " =8 

" = 6* 

8 . 

" = 4 " =2 " = 1 

M « 8 * 

9. 

FILE - 1, NSA = 1,100,1,1,359,1 * 


10 . 

" 2 " 1,110,1,1,320,1* 


11 . 

FDATE(l) = 01,25,72, MISSON = 7* 


12 . 

MPXMUX = 1, SCHAN " 3, AIRCHN = 6 , 

GAIN = 10* 

13. 

» =2 " =4 " =5 

" = 12 * 

14. 

FILE = 2, NSA = 1,100,1,1,320,1 * 


15. 

" = 3 " = 10,60,1,1,320,1 * 


16. 

FDATE(l) = 02,13,72, MISSON - 18* 


17. 

MPXMUX = 1, SCHAN - 0, AIRCHN = 3, 

GAIN * 14* 

18. 

" = 2 " = 3 " =1 

" - 14* 

19. 

FILE = 3, NSA = 10,60,1,1,320,1* 


20 . 

MODE1 = $F1NAL$ * 


21 . 

MODE1 = $START$* 


22 . 

TAPE = $PRINT$, MISSON - 72* 



1. A maximum of these files are to be processed with a maximum multiplexor 
channel number of five. The tape data base is assumed to be tapes 799 and 
800. A data base will be started on tape 799. 

2. This is the input data required by PROCESS under POINT. Tape 706 will be 
mounted on unit 4, and file 1 will be accessed. Every line from 1-100 will 
be processed and every point from 1-359 will be used. This starts the calcula- 
tion of the means. 

3. & 4. The date, mission, and comment given will be associated with this file, 

5-8. Multiplexor channels 1, 2, & 4 correspond to spectral channels 4, 5, & 2. 

Multiplexor channel 3 is to be ignored. File one has only 4 channels and all 
four channels must be specified even though only three of the four channels 
are actually to be used. The gains associated with the channels are as listed. 
All gains on the printout will be normalized to gain 100. 
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9. This is the input data required by PROCESS under POINT* This causes the previous 
file to be reaccessed so that the file can be read a second time. After the 
first pass through the file is completed, all the means have been calculated. 

This second pass is required in order to calculate the standard deviations 

for each of the corresponding means. 

10. This is more input required by PROCESS to access the second file and start 
the first pass through the data. 

11. Again the date and mission specified will be associated with this file. 

Since COMMNT wasn't specified, it is assumed to be blank. 

12. & 13. In this case both of the channels in file 2 are being used. 

14. This caused the file to be rewound and the second pass started. 

15. This encourages PROCESS to access the third file and start pass 1. 

16. Same as 11. 

17 & 18. Only multiplexor channel 2 is being referenced, but as before, both of 
the channels on the tape must be specified. 

19. Pass 2 on the 3rd file. 

20. This tells POINT to perform STEP (6) in the program. That section of the 
program is responsible for plotting and tabulating the results as well as 
updating the data base. 

21. This causes POINT, to reinitialize CAL IE. 

22. This tells CALIB. that the tape data bases (indicated by TBASE1 & TEASE2) 

are to be simply plotted and tabulated with absolutely no processing to take 
place. After the plots & tables are finished, the program returns control to 
the executive system. Therefore no more data dards must ever follow this card. 

ERROR COMMENTS 
ERR0R1 . LAMPB .L. NA. 

The beginning reselm location of the lamp pulse is less than the beginning 
re b elm for the line (as specified for PROCESS). 

ERR0R2 . LAMPE .0, NB. 

The end reselm location of the lamp pulse is greater than the ending reselm 
for the line (as specified for PROCESS). 
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ERR0R3. DARKB .L. NA. 

Same as ERRORl. except in reference to the dark area. 

ERR0R4. DARKE .G. NB. 

Same as ERR0R2 except in reference to the dark area. 

ERRORS . SOLARB .L. NA. 

Same as ERRORl except in reference to the solar pulse. 

ERR0R6. SOLARE .G. NB. 

Same as ERR0R2 except in reference to the solar pulse. 

ERR0R7. 13 .L. CHANNL -L. 1 

The maximum # of channels to be used in computations is either greater than 
13 or less than 1. 

ERRORS. FILMAX .L. 1. 

The maximum // of files to be processed before plotting and tabulation is less 
than 1« 

ERROR 9 . 24 .L. SCHAN .L. 0 

A spectral channel # greater than 24 or less than 0 was to be associated with 
a multiplexor channel. 

ERRORl 0. MORE FILE THAN FILMAX. 

The input required by PROCESS under POINT specified another file, but this 
file is one more than the number specified by FILMAX originally. 

ERRORl 1. INVALID MONTH. 

The number in FDATE (1) refers to a nonexistent month (less fhan zero or 
greater than 12). 

ERRORl 2. INVALID DAY. 

The number in FDATE(2) refers to a non-existent day (less than zero or 
greater than 31). 

ERR0R13. INVALID YEAR. 

The number in FDATE (3) refers to a nonexistent year (less than zero or 
greater than 99) . 

ERRORl 4. INVALID MISSON NO. 

The number assigned to MISSON is less than zero. 
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ERR0R15. SOLAR REF .L. 5 RES ELMS WIDE. 

The value of S0LARE- SOLARB+1 is less than 4. There must be at least 4 reselms 
In the solar pulse in order to integrate it. 

ERR0R16. GAIN .L. 1. 

The gain specified for the preceding multiplexor CHANNL was less than 1. 

ERR0R17 » INCORRECT TAPE MOUNTED. 

The program was not able to read a correct header label on the tape just 
mounted. This can be caused by mounting the wrong tape or by specifying 
TAPE = $OLD$ (or else letting TAPE default to $OLD$) when the tape data base is 
non existent or by specifying the wrong tapes with TBASE1 & TBASE2 . 

ERR0R18. LAMP E- L AMPB+1 .G. 100. 

The value of LAMPE-LAMPB+1 is greater than 100. No more than 100 reselms 
can be handled in the lamp pulse. 

ERR0R19 . SOLARE- SOLARB+1 -G. 100. 

The value of SOLAR- S0LARB+1 is greater than 100. No more than 100 reselms can 
be handled in the solar pulse. 

ERR0R20. INVALID MPXMUX CHANNEL. 

The value of MPXMUX is greater than 13 or less than 1* It is assumed that there 
will be no more than 13 and no less than 1 multiplexor channel on a tape. 

ERR0R2 1 . INSUFFICIENT STORAGE. 

Since this program uses dynamic storage allocation for its main incore data 
base, the possibility exists that some-one will try to process more data and 
consequently compute more statistics than can physically fit into the machine. 

The main in-core data base is dimensioned FILMAX*CHANNL*ll+FILMAX*10+4. 

In addition, there are two more arrays, each possibly as long as FILMAX+ENTRYS (ENTRYS 
is the largest number of records in any file of the tape data base.) If storage 
limitations are a problem: 

1. Process fewer files in a set. 

2. Process fewer multiplexor channels. 

3. Don’t attempt to merge in a large, old tape data base. Instead, 
create a new tape data base. 

See storage requirements section of this memo. 
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ERR0R22. LAMTOP . L. 0. 

Less than 1 reselm is to be used for computations at the top of the 
lamp pulse. 

ERROR23. LAMTOP .G. LAMP WIDTH. 

More reselms are to be used in computations than exist within the de 
fined limits of the lamp pulse. 

ERR0R24. NO DARK AREA PROVIDED. 

The value of DARKE- DARKB+ 1 is less than 1. The program normalizes all the values 
of the points in the lamp pulse and solar pulse by subtracting the average level of 
the dark area from the previous line. If no dark level is specified, this 
normalization can T t take place. 

SUBROUTINES NEEDED 

Following is a list of external subroutines needed for the correct opera- 
tion of this program. 

ASSIGN. - VERSION 2.2 This performs all the dynamic storage allocation in the 
a core box. 

CANCEL. - Part of ASSIGN. 

ERROR. - System error processing routine 
GETSPO. - Part of ASSIGN. 

LINK. - Routine for providing a connection to the point processing routine 
PTPROC . (Part of POINT.) 

MOUNT. - Routine for handling all tape mounts. 

MOVER. - A routine used for flipping an array end-for-end before output by 

All four routines are used in setting up and printing the point 
plot of the means and standard deviations. 

The main program that repetitively calls CALIB. with different levels 
to progressively process the data (see POINT, writeup). 

A routine called from POINT, that handles all I/O with data tapes 
(that are in ERIM format). 

A routine for setting up the dimensioning and subscripting for the 
dynamically allocated arrays. 


PL0T1. 
PLOT 2. 
PLOT 3. 
PLOT 4 . 

POINT. 


PROCESS. 


SETDIM. 
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SETEOF. A routine to perform branching when an end-of^file is encountered 
on input. 

SQRT. This routine calculates square roots. 

ZERO. This routine inserts either a binary or floating point zero In any 

specified variables or arrays. 

TAPE DATA BASE FOR CAL IB. 

This data base is used for storing all of the major results calculated by 
CALIB. For every file processed by CABIE. , the following information is saved: 

DATE 

SPECIAL COMMENTS 
MISSION NUMBER 
MEAN DARK LEVEL 
MEAN LAMP LEVEL 
MEAN LEVEL OF MIDDLE OF LAMP 
MEAN RESELM NUMBER OF MIDDLE OF LAMP 
STANDARD DEVIATIONS OF PREVIOUS 5 STATISTICS 
At the beginning of both tapes is stored the bin number of the tape (which is 
read to make sure the proper tape is mounted) and the generation # of the tape. 

The generation // is set to 1 when the tape data base is first generated (when 
TAPE = $NEW$) . The generation number is incremented by 1 each time the data base 
is updated (when TAPE = $OLD$; the data base is never altered when TAPE=$PRINT$) . 
When the data base is being updated, the tape with the highest generation number 
is the input tape. The tape with the lower generation number then becomes the 
output tape. This output tape will be given a generation number, 1 higher than the 
input tape (except when the input tape has a generation number of 1; then the output 
tape will be given a generation number of 3j Generation number 2 is never used ). 
The results computed from the files just processed by CABIB. are merged with the 
input tape data base and written on the output tape. 

The organization of the tape data base is as follows: 

1. There are 24 files representing up to 24 different spectral channels. 

2. Each record represents the results calculated by CABIB. from 1 channel 
of 1 file. 
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3- Each file has a 1-word header record containing the number of 


records in the file excluding the header record. 

4. If there are no results from CALIB. or from a previous data base referencing 


a channel, the corresponding tape data base file is simply terminated by an EOF 
and no header record is written. 


5, The records are in ascending order of date. 


6. Each tape of the data base has a 2 -word label containing the bin number 


and a generation number. 
FILE RECORD WORDS 

11 1 

2 

2 1 

3 1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


DESCRIPTION 
BIN NO. 

GENERATION NO. 

NO, OF RECORDS IN 1ST FILE. IF NONE, THIS 

RECORD IS AN EOF 

DATE 

MISSION # 

MEAN OF DARK LEVEL 

MEAN OF NORMALIZED LAMP LEVEL 

MEAN OF NORMALIZED LEVEL OF THE MIDDLE 

RESELM IN THE LAMP PULSE 

RESELM IN THE LAMP PULSE (RELATIVE TO LAMPB- 1) 
MEAN OF THE INTEGRAL OF THE NORMALIZED 
SOLAR PULSE 

STANDARD DEVIATION OF WORD 3 

ir ii ll it ^ 

ir n II n ij 

it it tl n g 

MEAN OF THE RESELM NUMBER OF THE MIDDLE 
STANDARD DEVIATION OF WORD 11 


13-20 SPECIAL COMMENTS 

[All means & standard deviations are integers multiplied by 100, Solar 
pulse statistics are X10] . 
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FILE RECORD DESCRIPTION 

1 4 same as record 3 

, [The form of record 3 is repeated for as many 

records as the number indicated in record 2]. 

2 1 No, of records in 2nd file. If none, this record 

, . is an EOF, 

[The form of file 2 is the same as file 2 starting with 
record 2, This is repeated for 24 possible files (each 
representing a different spectral channel, see description 
of "SCHAN" in parameter definitions)], 

STORAGE REQUIREMENTS 

The CALIB. program requires 7292 decimal locations before any 
dynamic storage allocation takes place. The formula for determining how 
much additional storage will be allocated for a particular data set is 
as follows: 

MAIN- IN-CORE DATA BASE 
FILMAX*CHANNL*ll+FILMAX*ll+5 

In addition, more storage is allocated depending upon the size of the tape 
data base. The maximum amount that will be allocated in addition to the 
preceding quantity specified is as follows: 

(MAX[ENTRYS] + FILMAX) *21+3 

MAX[ENTRYS] = the maximum number of missions in any one file 
of the tape data base. 
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MAC (C5 AUG 1565 VERSION) PHCGR AN LISTING 


EXTERNAL FLNCTICN C At I 8 . ( ST A T , CA T E , N SS ICN ,CNMNT , 

L TAC,XX,YY,YYY) 

NORMAL MODE IS INTEGER 
REFERENCES CN 

ELCATING PCINT F1,F2,F3, MCU3), NEANCI13), ME AM 1 1 3 ) , ME ANLM ( 13 ) , 

1 MEANMim, MEANS ( 1 3 ) , SCRT., XX, Vyy, STAT, G A IN I l 3 ) , 

2 T0TCI13), TCTL113), TLTLNI13), TCTMI13I, TCTSI13), 

1 PN TL ( LCC*13) » PNTS(100#I3> , C22, CCEF, F DATUM I 24 ) 

C1MENSICN FARAMSIB) • FCATE13) , NSCALE(A) , BCC(A) ♦ CCMMNTI8I, 

1 CHNNLS (13), ERRUC09I, TITLE(32>, OCM(B), CUTI2C), 

2 P T ( 1 3 ) 

LCL 1 VALENCE [MPXMLX ,FCATt ) , (SCHAN.HCO) , ( A I R CE N , CCN KN T) , 

1 (DATUM .ECATLN > 



VECTOR 

VALLES 

BCD < 1 )= $C», iLt,*Mt,tSi 




VECTOR 

VALUES 

Ml = -1 





VECTOR 

VALLES 

N'T 2 * -2 





VECTOR 

VALUES 

TITLE(L) ^ $F 1 i 

.71-. 73$ * $ E 1 * 

.33-. 38$ 

1 



$ B2 i 2,0-2, 6$ , 

$63* 

1.0-1. 4$ 


2 



$01* 2. 0-2. 6$ * 

£L 2 * 

1.5- 1.8$ 


3 



$C3. i.Q-l.4t ♦ 

$C2 * 

.67-. 94$ 


4 



$C4 * .62-, 70$ * 

$C5, 

♦5B-.64$ 


5 



$C6* .55-. 60$ * 

$C 7 * 

.52-. 57$ 


6 



$C8* .50-. 54$ * 

$C9, 

.48-. 62$ 


7 

VECTOR 

VALLES 

$ClOt .46-. 49$ * 
ERR ( 0C1 J = 

$C 1 2 * 

.41-. 48$ 


l 


$H # 1*****ERRCR1. LAMPE 

.L. N A . 

» * 

$ 


VECTOR 

VALLES 

E R R ( 0 4 3 ) = 




2 


$H i 1 *****EkRCR2. LAMPE 

.G. NR. 

» * 

$ 


VEC TOR 

VALUE S 

ER R < 085 ) - 




3 


$H*i*****ERRCR3* DARKB 

• L • N A . 

* * 

$ 


VECTOR 

VALLES 

ERR 1127) = 




4 


tH * !***♦**: RHCR4. CARKE 

.G. N6. 

• * 

t 


VEC TOR 

VALLES 

ERR ( 169 ) - 




5 


$H • SCLARb 

.L. N A 


$ 


VECTOR 

VALLES 

ERR (211 > = 




6 


w l*****t«KCKb* SCLARE 

.G. NB 


$ 


VEC TOR 

VALLES 

ERP<253> = 




7 


* H i i*****frkcR7. 13 . L • 

CHANM 

.L . 1 . 1 * 

$ 


VECTOR 

VALLES 

ERRI295) = 




e 


$ H •!*****£ rrcR 8. FILMAX 

.L. i. 

■ * 

$ 


VECTOR 

VALUES 

E R R ( 3 3 7 ) = 




9 


$HU*****ErtRC H9. 24 . L • 

SCFAN 

.L. Q . •* 

$ 


VECTOR VALLES ERRI379) = 

C tH ' 1 ** ** *CRRCR 1C. MORE FILES THAN' FILMAX.'*S 

VECTOR VALLES ERRIA21) = 

1 *H , l*****EKRCKLl. INVALIC MONTE.*# t 

VECIOR VALLES ERR ( A63 ) = 

2 $H* 1*####ERRCR12. INVALIC CAT.*# $ 

VECTOR VALLES ERR ( 5C5 ) = 

3 *H*L#*##*ERRCR13. INVALIC YEAR.** $ 

VECIOR VALUES ERRC547) = 

A ih ' 1*****ERRCR1A . INVALID MISSION NUMBER. •♦$ 

VECTOR VALLES ERRI589) = 

5 iH'l*****ERRCRlS. SCLAR REF .L. 5 RESEIMS**$ 

VECTOR VALUES ERRI631) = 

t th* I*****ERRCRL6. GAIN .L. I.»* t 

VECIOR VALUES ERRI673) = 

7 $H* 1*#***ERRCR17. INCORRECT TAPE MCUNT E C . * *S 

VECTOR VALLES ERKI715) = 
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c 




8 $ H * l*****ERRCR18. LAM PE-LAMPB* 1 .G. 100.'* * 

VECTOR VALUES ERR(757) * 

9 $H'1*****ERR0R19. SCLARE-SCL ARB+1 .G. 100'** 

VECTOR VALUES ERR(799I = 

C *H'l*****ERRGR20. INVALID MPXMUX CHANNEL. • ** 

VECTOR VALLES ERR ( 8AI ) = 

1 sh'1*****ERRCR21. INSUFFICIENT STORAGE. ■* $ 

VECTOR VALUES ERR(883) * 

2 *H , l*****ERRGR22. LAM TCP .L. 0.'* * 

VECTOR VALUES ERR(925> = 

3 *H* l*****ERRCR23. LAMTGP .G. LAMP N I CTH. 1 * S 
VECTOR VALUES ERR ( 9f>7 ) = 

« *H'1*****ERR0R29. NG DARK AREA PROVIDED.'* $ 

STATEMENT LABEL STEP I 7 ) 

VECTOR VALUES MSG1(C> = JMEAN LEVEL* 

VECTOR VALUES MSG2I0>= tDEVI AT ICN LEVEL* 

ERASABLE JUNK ( 22 5 > , CUN I T .CREEL , CF I LE ,GL I NE , TWR I TE, OCH AN , 

1 CNWCS,PRCG,MODE1,MOOE2,UNIT,CALINE.CCLINE,NSA,N$B, 

2 KS,NA »NB,KP*ID(1.> , TFLAG , I PACK, RESE RV ( A ) , Q FACTR ( A9>, 

3 CTITL2(19},QTITLE(I9) , CL I ST ( 19 > , CS P ARE I A6 ) , CFL AG, 

A CDANGtCBANG,CRP.CNA,CNSS,CNCHAN,gMCDE,CRECA,ORECC. 

5 QFILE, CREEL, CNWRDS 

ERASABLE DA T A ( 42 3 ) , 1 TE ST .CSTART ,DSTART , NV , NX ,NC , l , I P, 

1 1CP, NOP, NEXT, START2,REA0T2,READL2,R£ACP2,LABEL2, 

2 EXTRAI2C) , DATUM (24) , I CCDE 1 2A ) , ICHANI2A), IMAGEI570 ) 

WHENEVER NEXT .G. 6, FUNCTION RETURN 

TRANSFER TC STEPINEXTI 
STEP ( 1 ) L INK. IPTPRQC. I 
L « L ♦ OCHAN 
CHANNL = 13 
FILMAX = 1 
LAPTOP = 5 
LAMPB * 1C 
LAMPE = 6C 
DARKB * 185 
DARKE « 225 
SCLARB = 3 1C 
SOLARE = 355 
PARAMSIC) = 0 
FILES = C 
TAPE * 10LD* 

PASS = 1 
MISSCN = 0 
TBASE1 = 799 
T6ASE2 = 8CC 

READ AND PRINT DATA TAPE, CHANNL, FILMAX, MISSCN, LAMTOP, TBASE1, T6ASE2 
WHENEVER TAPE .E. 1PR 1 NTS 

WHENEVER MISSCN .L. 1, TRANSFER TO ERR1A 
TRANSFER TC STEP (6) 

END OF CONDITIONAL 
WHENEVER FILMAX .L. 1 
TRANSFER TC ERR8 

CR WHENEVER CHANNL .G. 13 .CR. CHANNL .L. 1 
TRANSFER TC ERR7 
CR WHENEVER LAMTOP . L. 0 
TRANSFER TC ERR22 

CR WHENEVER LAMTOP .G. LAMPE - LAMPB + I 
TRANSFER TO ERR23 
END OF CONDITIONAL 
SIZE = FILMAX* CHANNL*L1 
KK = ASSIGN. (0,01 

WHENEVER KK .L. ( SI ZE + F I LMAX*10*A ) , TRANSFER TO ERR21 

KK = ASSIGN. (CALIB.,PARAMS(0) .SIZE. ..3, FILMAX, FILMAX, 02 
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STEPU ) 


LP IB 


STEP! 3 ) 
STFPIAI 


FILPAXA8...2 I 


SETCIP. (CPPNT,FILPAX,8) 

StlOTP. ( STAT ,FILPAX,CHANNL,11) 

ZERO. < TC TC ( 11 . .. TCTC (CHANNL I . TCTLI1I .. .TCTLICHANNL ) t 

ICTLP I L ) . . . TCI LP (CHANNL I * TCTS t l I • • • TC TS (CHANNL I * 

TOTP ( 1 ) . . . TCTP (CHANNL I » 

PC<1) ...PD (CHANNL! * STATU) .. - ST AT (S IZE) > 

ELNUTION KE ILRN 

ZEHC.(PTL»PJ$»Pf...PT(13>«LINES»FLAG) 

WHENEVER PASS .C. I ,, . 

ZERL.(FCATE(1).«»FCAT6(3)«PISSCN*CHNNLS(1)...CHNNLS(13!) 


1 


<KK*l,l,KK.G.8,CCPPNT(KK) = $ *1 

«EAC AND PRIM DATA FUATftllt PISSCNr CHNNLSllIt 

LAPPfl»lAPPE,DARKB»E.ARKE*SCLARB,SCLA«E 


CCPPNT t 1 ) » 


THRCUC.H LP16, FCR K=l»l» K .G. C ACM AM 

REAU AMI PRIM DATA PPXPUX, SCHAN, AIHCHN, GAIK 
WHENEVER SCHAN .G. 2A .CR. SCHAN .L. 0, TRANSFER TO ERR? 
WHENEVER PPXPLX .0. 13 .CR. PPXPUX .L. 1, TRANSFER TG ERR20 
WHENEVER GAIN .L. 1. TRANSFER TC ERRlb 
GAIfv(K) = GAIN 


CHNNL S ( PPXPUX ) = SCHAN 

WHENEVER SCHAN .E. 0. TRANSFER TC LP18 

PRINT FORPAT $ 1H+ , T6U . H 1 SPEC TK A L B ANC= *»2C6*t* 

TlTLE(SCHAN*2-lf ,TITLE(SCHAN*2) 


CCNT l ME 
CHANTP = C NO HAN 
WHENEVER LAPP8 .L. NA 
TRANSFER TC ERR 1 


OR WHENEVER LAPPE .G. MB 
TRANSFER TC ERR2 
CR WHENEVER DARKB .L. NA 
TRANSFER TC ERR3 
CR WHENEVER DARKE .G. NB 
TRANSFER TC ERRA 
DR WHENEVER SCLARH .L. NA 
TRANSFER TC ERRS 
OR WHENEVER SCIARE .G. N8 
TRANSFER TO ERR6 

OR WHENEVER FDATEtl) .G. 12 .CR. FCATEI1) .L. 0 
TRANSFER TC EftRll 

OR WHENEVER FDATE12) .C. 31 .CR. FCATEI2) .L. 0 
TRANSFER TO ERR12 


CR WHENEVER FCATE ( 3 ) «L« 0 

TRANSFER TG ERR13 
OR WHENEVER PIISSGN . L. 0 
TRANSFER TC ERR14 
OR WHENEVER L A P PE-L APP8* l .G. 100 
TRANSFER ID ERR18 

CR WHENEVER SC L A R E - SC L AR B+ 1 .G. 100 

TRANSFER TC ERR 19 

OR WHENEVER DARKE - UaRKB + 1 -t. 0 
TRANSFER TC ERR2A 
END CF CONDITIONAL 
ENC CF CONDI TICNAL 
FLNCTION RETURN 
FUNCTION RETURN 
LINES = LINES + 1 
WHENEVER FLAG .E. C 
FLAG = L 
FUNCTION RETURN 
ENC CF CCNDI TICNAL 

THROUGH L PS » FOR K=1,1,K .G. CHANTP 

WHENEVER CHNNLSIK) .E. 0, TRANSFER TC EPS 
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MAX ^ C 

THROUGH LP6, FCR KK= L * 1 * KK .G . PTL 

WHEN E VER PNTL ( KK , K ) «G • MAX t MAX = PNTl(KK t K) 


LP6 

CONTINUE 





THROUGH l P7 , FOR 

KK=1,1,KK .G. PTL 




WHENEVER MAX 

/2 .L. PNTL (KK, K ) » TRANSFER 

TC 

OUT 1 

L P 7 

CONTINUE 





KK = 1 




QUT1 

LEFT = KK 





THROUGH LP8, FOR 

KK=PTL»-1, KK .L. 1 




WHENEVER MAX 

/2 .L. PNTL ( KK, K ) » TRANSFER 

TC 

0UT2 

LP8 

CONTINUE 





KK - PTL 




OUT 2 

RIGHT = KK 





MIDDLE = (LEFT ♦ RIGHT)/2 
LEFT = MIDDLE - LAMTOP/2 
RIGHT = MIDDLE ♦ LAMTCP/2 
WHENEVER LEFT .L. 1# LEFT = L 
WHENEVER RIGHT .G. PTLt RIGHT = PTL 
P T ( K ) = PT { K ) ♦ RIGHT - LEFT + 1 
WHENEVER PASS . E . L 

(KK=L£FT t I tKK.GwRlOHT.TCTL <K ) = TCTL(K) + PNTL ( KK ♦ K ) ) 

TCTLM(K) = TOTLM(K) ♦ PNTL (MIC CLEfK) 

TCTM(K) - TCTM(K) + MIDDLE 
OTHERWISE 

(KK^LEFT, L f KK.G.RIGHT,TCTLtK ) = TCTUKl ♦ 

I ( PNTL ( KK * K } - ME ANL ( K ) ) • P *2 I 

TCTLM(K) = TOTLM(K) ♦ { PN T L l M I DC L E , K ) - 
1 MEANLMIK) )«F.2 

TCTM(K) = TOTM(K) ♦ (MIDDLE - MEANMIKU.P.2 
END OF CONDITIONAL 

WHENEVER PTS .L. 5 . AND* SCLARE - SCLARfi ♦ 1 * NE * 0, TRANSFER TO ERR 1 5 
EXECUTE SIMPSN. 

WHENEVER PASS .E. 1 

TCTS(K) = TCTS(K) + Q22 
OTHERWISE 

TCTS(K) = TCTSIKI + <C22 - M E AN S ( K ) ) . P . 2 
END OF CONDITIONAL 

LPS CONTINUE 

P1L = 0 
PTS = C 

FUNCTION RETURN 

STEP(5) WHENEVER PASS .E. 2 * TRANSFER TC STCCEV 

FILES = FILES ♦ 1 

WHENEVER FILES .G. FILMAX. TRANSFER TC E RR 10 
THROUGH LPS f FOR K=l,i* K .G. CHANTP 

WHENEVER CHNNLS(K) .E. O f TRANSFER TC LP9 
MEAND(K) - TOTD(K)/( < D ARKE-CARKe+1 ) *L 1NES ) 

MEANL(K) = TCTL(K)/PT(K) 

MEANLM(K) = TOTL M ( K ) / ( L I NES- l ) 

MEANS ( K ) = TGT$(K)/(LINE$“1) 

MEANM(K) = TOTM (K) / ILINES-l) 

STAT(FlLEStXtl) = MEAND(K)*100 
STAT(FILES*K»2) = ME ANL < K ) *10 0/G A IN ( K > 

STAT (FILES »K t 3) = ME ANLM ( K J *100/GA I N < K ) 

STAT(FILES*K»4) = ME AN$ ( K )* IO/GA IN ( K ) 

STAT(F1LES»K,9) = ME ANM l K )* 100 
LP9 CONTINUE 

PASS = 2 

OUT 3 ZERO. (TCTD( L ) . . . TC TD ( CHANNL ) , TCTL <U .. .TCT L ( CHANNL > * 

1 TGTLM( 1 )•«* TOT LM (CHANNL > , TOTS U> ** .TCTS < CHANNL ># FLAGt 

2 TCTM(i)...TOTK< CHANNL I ) 

FUNCTION RETURN 
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STCCEV THROUGH LPLC. FOR K=l.l. K .G. CHAMP 

STATIFILES.K.llI = CHNNLSIK) 

WHENEVER CHNNLS ( K ) .6. 0« TRANSFER TC LP10 
FI = TGTO (K I /< (CARKE-DARKB+l ) *L I NES I 
SIATlFILES,K f 3)=SCRT.(Fl)*100 
FL = TGTL(K)/PT(K) 

STAT(FlLES,Kt6l*SCRT.(FlI*100/GAlMKI 
FI = TCTlMKI/tLtNES-1) 

STATIFlLES.Kt?) =SGRT.(F1I*100/GAIN(K) 

FI = TCTSIK) /UINES-II 

$TAT<FILES»K,8>=SCRT.(F1)*10/GAIMK) 

FI = TOTR (K) /(LINES-1 > 

$IAT(F[LES,K,ICI = S CR I . ( F 1 I * 100 
LP 1C CONTINUE 

PASS = l 

UATfc l F I LE S ) = FOA TE ( 3 I *10000 ♦ F0ATE(1)*1CC + FCAT E ( 2 ) 
MSSICNIFILESI = PI SSCN 

(KK=l,l,KK.G.8»CPMM (FILES.KK) = CCPMNTIKK)) 

TRANSFER 1C CUT3 
STEP(E) PAHARStCl = A 

WHENEVER TAPE .E. SPRINTS. TRANSFER TC PNT 
KK = ASSIGN. (0, Cl 

WHENEVER KK .L. FILES+1, TRANSFER TC ERR21 
KK = ASSIGN. tCALlB. .PARAPS (01 .FILESI 
IKK = 1, l.KK.G. FILES, IAG1KK) = KK) 

THRCCGH L P 12 i FCR KK=l,l,KK .G. (FILES - 1) 

THROUGH LP13, FCR J J= [ KK+ l ) , 1 » J J .G. FILES 
WHENEVER OATEITAGIKKI) .GE. OA TE ( T AG ( J J I ) 

TEMP * TAG ( KK ) 

TAG(KK) = TAG(JJ) 

TAGUJ) = TEMP 

CND CF CCNCITICNAL 
LP 1 3 CONTINUE 

L P 1 2 CONTINUE 

PNT WHENEVER MT1 .NE. THASEl 

MOUNT, t 2 . T8ASE 1 ,tINS) 

MT1 x TBASEl 
PNC OF CCNOI IICNAL 
WHENEVER MT2 .NE. TBASE2 
MOUNT. (3, TBASE2, SINS) 

MT2 = TBASE2 
ENC OF CONDITIONAL 
WHENEVER TAPE .E. SNEWS 
IN = <5S9 
OUT = 2 
BINCUT = 7SS 
BiNCUT = TBASEl 
GENCUT = 1 
TRANSFER TC WRTLAB 
END OF CONDITIONAL 
SLTEOF. IERR17) 

READ BINARY TAPE 2. B1N2.GEN2 
WHENEVER GEN2 .E. 1 
GEN2 = GEN2 + 1 
TRANSFER TC IN2 
ENC OF CCNCITICNAL 
SETEOF. ( IN2I 

READ BINARY TAPE 3. BIN3.GEN3 
WHENEVER 8 1 N 3 .NE. TBASE2. TRANSFER TC IN2 
WHENEVER 0 IN2 .NE. TBASEl, TRANSFER TC ERR 17 
WHENEVER GEN2 .G. GEN3 
IN 2 IN * 2 

OUT = 3 
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WRTLAB 


ARND1 

ENTRYC 


SETUP 

ARNC2 


RDOID 

RDNEWi 


tilNIN = TBASE1 
BINCUT = TBASE2 
GENCUT * GEK2 * l 
C THERM SE 
IN = 3 
CUT * 2 

8ININ * TBASE2 
BINCUT = TBASEl 
GENCUT = GEN3 + 1 
ENG OF CCNOITICNAL 
REMIND TAPE CUT 
WHENEVER TAPE *E. SPRINT! 

Z E R C . (NSCALE (1) ...NSCALE(A) ) 

N SCALE ( 0 ) = i 

EXECUTE PLCTi. (NSC ALE, 5, 10,6,101 
F2 = l. 

WHENEVER PISSCN - 61 .G. 0# F2 * MISSCN - 61 
TRANSFER TG ARNC1 
ENC OF CCNCmCNAU 

PRINT ON LINE FCRttAT SH* AT ENC CF JCB, PLEASE REMOVE •, 

I H'RING FRCM TAPE *,13*$, BINCUT 

WRITE BINARY TAPE CUT, B I N GU T , G E NCUT 
ZERC . I NSC ALE (1)... N SCALE (4 ) ) 

NSCALEIC) * 1 

EXECUTE PLCT1. INSCALE, 5 ,10,6 ,10) 

F2 = i. 

WHENEVER I M S S I CN ( F I LE S ) - 61) * G * 0, F2 = M S S I CN ( F I L E S ) - 61 
BLOCK = GETNXT. (8 > 

LEVEL = GETNXT. (BLOCK+4 ) 

ACCRES = GETNXT. (BLCCK > 

THROUGH IP I A , FCR K=i,l, K .G. 24 
WHENEVER TAPE .E. SNEWS 
ENTRYS = 0 

WHENEVER TAPE .E. SPRINTS, TRANSFER TC LP14 
TRANSFER TC SETUP 
END OF CCNOITICNAL 
SETEOF . (ENTRYC) 

READ BINARY TAPE IN, ENTRYS 
WHENEVER TAPE ,E. SPRINTS 

WHENEVER ENTRYS .E. 0, TRANSFER TC LP14 
FI = MISSCN 
TRANSFER TC ARND2 
END OF CCNOITICNAL 

WHENEVER ENTRYS+FILES .E. 0, TRANSFER TC WRTECF 
FL = MSS1CN (FILES) 

PARAMS(C) * 5 
KK = ASSIGN. (C,0) 

WHENEVER KK .L. (< EN TR YS + F I LES ) *2 1 + 3 I , TRANSFER TO ERR2 1 
KK = ASSIGN. (CALI6. ,PARAMS (0) ,ENTRYS+FlLESt 
l (ENTRYS+FILES)*19...2,ENTRYS4FJLES) 

SETDIM. (YY, ENTRY S+ FILES, 19) 

II * C 
JJ = C 

WHENEVER TAPE .E. SNEWS .CR. ENTRYS .E. 0, TRANSFER TC CMPNEW 
SETEOF. (CMPNEW ) 

READ BINARY TAPE IN* CD A T E ,C K S$ , C KC ,CKl , CMLM , CMS , GC D , CCL , OCLM , CDS , OPM , 
1 COM ,OCM ( 1 ) . . . CCP ( 8 ) 

U = II + 1 

WHENEVER II .G. FILES, TRANSFER TC CMPCLC 
WHENEVER CCATE .L. DATE ( TAG ( II ) ) 

11=11-1 

EXECUTE INSERT. (SOLOS) 

TRANSFER TC RCCLD 
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CMPOLL 


PLOT 


LP16 


ARN03 


LP17 

ViRTEUF 
l P 1 4 


I 


1 

1 

2 

2 

3 


1 


\ 


OTHLKfcl St 

exeCLTf INSERT. ( > 
transfer tc RONEW 


GNU OF CCfSD I TICNAL 

WHENEVER II .G. FILES, TRANSFER TC PLCT 
EXECUTE INSERT. I SNEWS) 

TRANSFER TC OMPNEW 
EXECUTE INSERT. ISCLDS) 

READ G rti NARv' TAPE IN, CO ATE , CP SS , CMC , CML , CM.LM. , Cl' S , DC C , CCL. CCIM , 00 S , CP M , 
CDP ,CCM ( 1 I . . .OCM ( 8 I 


TRANSFER TC DMPCLC 

WHENEVER JJ -E. 0, TRANSFER TC WRTECF 
( KK=1 » 1 ,KK.G. (ENTRYS+FILES) ,YYY(KK)= -999.1 
[KK=l,l»KK.G.JJ,YYY(KK) = YY(KK,2)I 

EXECUTE PLCT2.IIMAGE,Fl,F2,YYlJJ,2)*2.,YYIJJ»2)/8. 

EXECUTE PLCT3.iaC0(2J,XXUI,YYY(n,JJ) T1TLE12*IU 

PRINT FORMAT SlHl ,T 33 ,H ■ CHANNEL ■ , 2C6* * ,T I T L E ( 2 *K 1>,T1TLE(2 

EXECUTE PLCT4. ( 10,MSG1 ) 

PRINT FORMAT UH0.T36,H'PISSICN NUMBER *S 
EXECUTE PLC T2. ( I PAGE ,F 1 »F2 ,20. ,0. ) 

(KK=l,l ,KK.G. (ENTRYS+FILES) ,YYY(KK)= -9 99. ) 

THROUGH LPlfc. FCR II=i,l, II ,G. a 

(KK = 1,1,KK.G.JJ,YYYIKK) = Y Y < K K , I 1) / 1C . 1 

EXECUTE PLD T'3. I8CDII 1-4 ) ,XX(1 ) ,YYYI I) ,JJ> 

PMNT N FORPAT tlHl,T33,H l CHANNEL *,2C6*$.T ITLEI 2*K-l 1 , T ITLE I 2*K ) 
EXECUTE PLCT4. I 15.PSG2 > 

PRINT FORMAT tlH0,T36,H»KISSltN NUPBER'** 

WHENEVER IAPE .E. SPRINTS, TRANSFER TC ARNC3 

PRINT FORPM IlHl, 132, H* LAPP* , T 6 0 , F * LAPP* ,T75,H 'LAMP LAPP'/IF 

h , LAPP',Sl,H , PlCClLE , ,S3»H , SCLAR , iS2)*T74»2(H , RESELM*» 

Sll/ 

F* GATE PISSICN • , 4 (H ' P E AN 1 , S3 ) » S 1 , 4 I F ' OEV • , S4 ) • 

T 75 ,H 1 PEA N DEV*//*t 
THROUGH LP 1 7 , FCR KK=l,l.KK .G. JJ 
UUTll ) = YY(KK,9> 

SUB 1 = Y Y I KK , 9 ) 

SUB I = ( SUB L - SU61/10000*10000)*100 + SUBl/lOCOO 


CUT ( 2 ) = XX (KK) 

(M=3,1,P.G.1C,CLT(M) = Y Y ( K K , P—2 ) ) 

(M=ll, 1,P.G. 12, CUT (M > = YY ( KK , P+7 ) ) 

<PI = I3,1,P.G.20,CUT(P) = YY(KK,P-3I) 

PRINT FCRPAT $ 1H , 2 ( 16 ,S1 ) , SI , 10 1 16 »S 1) I 1H ,S15,8Cfc*S, SUB1, 
G U 1 1 2 ) . • . CUT 1 20 1 

WHENEVER TAPE .E. SPRINT*, TRANSFER TC LP17 
EXECUTE PCVER. (CUT (l ) ,CUT (20) ) 

WRITE BINARY TAPE CUT, OUT (1) . . .CU T I 2C ) 

CUNTINUE 

PRINT FORMAT SH*-SCLAR MEAN ANC SCLAR CEV ARE X1C. 1 / 

H* ALL OTHER STATISTICS ARE XLOO.*** 

WHENEVER TAPE .E. SPRINTS, TRANSFER TC LP14 
END OF FILE TAPE CUT 
FREESP. I BLOCK, LEVEL, AC ORES) 


, 


CONTINUE 

WHENEVER TAPE .E. SPRINTS, SYSTEM. 
ENC OF FILE TAPE CUI 
ENC OF FILE TAPE CUT 

WHENEVER IN .6. 959, FUNCTION RETURN 
REWIND TAPE IN 
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InR I TE BINARY TAPE IN, BI 
RE k I NO TAPE IN 
REklND TAPE CUT 
FUNCTION RETURN 


ERR 1 

K = 1 



ERR2 

TRANSFER 
K a A3 

TC 

BLEH 

ERR 3 

TRANSFER 
K = 85 

TC 

BLEH 

ERRA 

TRANSFER 
K ^ 127 

TC 

BLEH 

ERR 5 

TRANSFER 
K = 169 

TC 

BLEH 

ERR 6 

TRANSFER 
K = 211 

TC 

BLEH 

ERR 7 

transfer 

K = 253 

TC 

BLEH 

ERR8 

TRANSFER 
K - 295 

TC 

BLEH 

ERRS 

TRANSFER 
K = 337 

TC 

BLEH 

ERRIO 

TRANSFER 
K = 379 

TC 

BLEH 

ERR 1 1 

TRANSFER 
K - A21 

TC 

BLEH 

ERR 1 2 

TRANSFER 
K = A63 

TC 

8LEH 

ERR13 

TRANSFER 
K = 505 

TC 

BLEH 

ERRIA 

TRANSFER 
K = 5 A 7 

TC 

BLEH 

ERR 1 5 

TRANSFER 
K = 589 

TC 

BLEH 

ERR 16 

TRANSFER 
K - 631 

TC 

BLEH 

ERR 1 7 

TRANSFER 
K = 673 

TC 

BLEH 

ERR 18 

TRANSFER 
K = 715 

TC 

BLEH 

ERR 19 

TRANSFER 
K = 7 57 

TC 

BLEH 

ERR20 

TRANSFER 
K - 799 

TC 

BLEH 

ERR2I 

transfer 

K = 8 A 1 

TC 

BLEH 

ERR22 

TRANSFER 
K = 883 

TC 

BLEH 

ERR 23 

TRANSFER 
K = 925 

TC 

BLEH 

ERR2A 

TRANSFER 
K * 967 

TC 

BLEH 


TRANSFER 

TC 

BLEH 

BLEH 

PRINT FORK AT ERR 
ERROR. 


N I N t GENCUT-1 
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INTERNAL FUNCTION PTPRCC. 

WHENEVER TAPE .E. iPRINTJ, FUNCTION RETURN 

WHENEVER CMCDC .E. 1* ( I » 1 , L , I . G.2A , K = DATUMU » , FCATUM U)=K ) 
WHENEVER IP .LE. CARKE .AND. IP -GE. CARKe 

WHENEVER IP .E. DA RKB , ( K= l , 1 , K .G .CH AN T P, MD ( K I = 0.) 

F 2 = DARKE - DARKB + 1 

<K=1,1,K.G.CHANTP, FI 3 C ATOM I K I , PC ( K ) = MOIKI ♦ F1/F2) 
WHENEVER PASS .t. 1 

(K 3 1,1,K.G.CHANTP,TCTCIK) = TCTCIK) ♦ CATUM(K)) 
OTHERWISE 

(K=l,l,K.G.CHANTP, TCTCIK! 3 TCTCIK! ♦ (OATUM(K) 

I - RE AND I K I ) .P.2 ) 

END OF CONDITIONAL 
FUNCTION RETURN 
OTHERWISE 

TRANSFER TO LAPP 
ENC OF CCND I T I CNAL 

LAMP WHENEVER FLAG .E. 1 .AND. IP -IE. LAMPE .AND. IP .GE. LAMP6 

PTL = PTL + 1 

tK=l,l,K.G.CHANTP,PNTL(PTL,K) = CATUMIK) - MCIKI! 

FUNCTION RETURN 
C THERWI SE 

TRANSFER TO SOLAR 
ENC OF CONDITIONAL 

SOLAR WHENEVER FLAG I .AND. IP .LE. SCLARE .AND. IP .GE. SOLARB 

PTS * PIS M 

(K=l,l ,K.G.CHANTP,PNTS(PTS,K) = CATUMIK! - MC(K)J 
FUNCTION RETURN 
CIHERWISE 

FUNCTION RETURN 
ENC OF CONDITIONAL 
ENC OF FUNCTION 


INTERNAL FUNCTION SIMPSN. 

WHENEVER PTS/2*Z .E. PTS, PTS 3 PTS - 1 
02 2 3 PNTSU ,K> 

COEF = 2. 

THRGUGH LP1L, FOR KK-2,1, KK .G. PTS-1 
•WHENEVER CCEF -E. A. 

CCEF = 2. 

OTHERWISE 
CCEF 3 A. 

END OF CONDITIONAL 

C22 = 022 ♦ CCEF APNTS I KK ,K ) 

LP1I CONTINUE 

022 = Q22 + PNTSIPTSfKI 
022 = 022/3. 

FUNCTION RETURN 
END OF FUNCTION 


INTERNAL FUNCTION INSERT. ( TYPE 1 

MTAG = TAG II II 

WHENEVER TYPF .£. J NEW t 

THROUGH LP19, FCR N= 1 , 1 , N -G. CHANNL 
WHENEVER STAT IMTAG,N»ll) .E. K 
JJ 3 JJ + 1 
XXIJJ! 3 MSS I CM FT AG I 
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LP 19 


(M=ltltK.G.8*YY(JJ»P) = STAT I PTAG»N t P ) ) 

YYUJ.9I = DATE IPT AG I 

(KK=10,1,KK.G.17,YY( JJ.KK) = CPPNTIPTAG ,KK-9II 
(M=16,1,P.G.19*YY( JJ.P) = $TAT lPTAG»N»P-9> ) 
FUNCTION RETURN 
ENO CF CONDITIONAL 
CONTINUE 

OR WHENEVER TYPE .E, SOLDI 
JJ = JJ ♦ 1 
XXIJJI = CPSS 
YYlJJtll = CPC 
Y Y ( J J » 2 > = CPL 
YYIJJ.3) * CP LP 
YY(JJ,«) = GPS 
YYIJJ.5I » CCD 
YY(JJ,6> = CCL 
Y Y( J J * 7 1 * CD LP 
YY(JJ,8> = CDS 
YYIJJ.9) = CD ATE 
EXECUTE POVER. (CCP (1 ) »0CP18) I 
(KK=1C»1,KK.G.17.YYI JJ.KKI = GCPIKK-91) 

YY(JJ»18> » CPP 
YYIJJ'I?) = CCP 
OTHERWISE 

FUNCTION RETURN 
ENC OF CONDITIONAL 
FUNCTION RETURN 
ENC CF FUNCTION 
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Appendix B 
PROGRAM DGNSTC. 

The purpose of this program is to calculate a few statistics from multispectral data tapes 
in ERIM format. These statistics can be used to observe the behavior of the data in certain 
specified areas of the scanline. This version has three options concerning the statistics ex- 
tracted from a scene. 

(1) The mean level and the standard deviation around this mean level of the data within a 
specified area relative to some defined dark area. 

(2) The mean reselm number and the standard deviation around this mean of the middle 
of some specified pulse. 

(3) A point plot of each absolute mean and standard deviation of a specified dark area 
plotted line by line and/or a point plot of each mean and standard deviation of a spe- 
cified area relative to the mean of the previously defined dark area plotted line by line. 

For example, with the use of option (l),the user can determine if there are large-scale changes 
(either absolute or relative to some other area) in an area parallel to the flight line. Option 
(2) gives the user the capability to determine how much a particular pulse is shifting back and 
forth along the scanline over a given set of scanlines. The plots provided by option (3) allow 
the user to see exactly how the mean of a specified area (either absolute or relative to some 
other area) changes from one scanline to the next over a given set of scanlines (trends or 
discontinuities). 

With all of the options, the user must specify the beginning and ending reselm numbers 
of two areas along the flight line to be analyzed. The first area must always be specified. 

The reselm numbers of the second area may be chosen such that it is zero if it is to be ignored. 
In addition, the maximum number of reselms within either pulse, the option type, and the maxi- 
mum number of scanlines must be specified for each set of files. Also, the user must provide 
the multiplexer channel number to be referenced and the spectral channel it represents along 
with the mission number, the date, and any special comment. 

Option (1) averages all the reselms within the boundaries of the first area for each line 
as well as over all of the selected lines. If the second pulse's width is non-zero, the average 
of the first area from the previous line is subtracted from each reselm of the second area. 

The result of this subtraction is averaged over all the selected lines. Then the program goes 
through the data set again and calculates the standard deviation of the reselm levels from the 
corresponding overall mean for the area. 
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Option (2) determines the reselm number of the middle of the pulse within the specified 
boundaries (see Fig. B.l). This is done for each line of the requested scene. After the aver- 
age of all these middle reselm numbers is determined, the program goes back and calculates 
the standard deviation of these reselm numbers from the overall mean. The method for de- 
termining the middle of the pulse is as follows: 

(1) FL = 0.5 

(2) The maximum height is determined [MAX] 

(3) The reselms closest to MAX * FL are found by first searching from left to right for 
the reselm whose level is just greater than MAX/2 and then performing a search 
from right to left. [X & Y] 

(4) If X & Y are equal to the boundaries of the pulse, FL = FL + 0.1. Repeat (3) above. 

(5) The mean of X & Y is found (X + Y)/2. This is the middle of the pulse. 

NOTE: If the boundaries given for the pulse don’t include reselms at or below the 
MAX * FL level, X & Y will end up being either of the boundaries given for the pulse. 
Consequently, the middle reselm number given for the pulse could be quite inaccurate. 

Option (3) determines the absolute average height and standard deviation from that aver- 
age of the first pulse for each line. In addition, the user may request the program to determine 
the average height and standard deviation of a second pulse relative to the first for each line. 
The program then plots. 

(1) The average height of the first pulse for each line 

(2) The standard deviation of the first pulse for each line 

(3) The average height of the second pulse relative to the first pulse for each line 

(4) The standard deviation of the second pulse for each line 

PROGRAM OPERATION 

STEP (1) PERFORMED THE 1ST TIME DGNSTC. IS CALLED OR WHENEVER PROCESS 
READS M0DE1 - $START$" 

This section of the program perforins the initialization necessary for 
POINT, and sets the default values for the input variables. Data cards 
are then read which set the beginning and ending reselm numbers of both 
pulses (only pulse 2 can have zero width; END(2) - BEGIN (2) +1=0), the 
maximum number of points in either pulse, the option type the program is 
to execute, and the maximum number of lines to be processed. The input 
variables are then checked for reasonableness, and the necessary storage 


is allocated. Control returns to point. 
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2 


FIGURE B.l. LAMP-PULSE VOLTAGE 
VERSUS RESELM POSITION, 
SIMPLIFIED 
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STEP (2) PERFORMED BEFORE EACH FILE IS PROCESSED. 

First the program zeros the point index array » the date array, the 

comment array, and the mission number. Then more input cards are read 

specifying the multiplexor channel, the spectral channel, the aircraft 

channel (optional), the mission number, the date, and any special comment. 

The input variables are also checked for reasonableness. The beginning and 

ending reselm numbers are also checked again for reasonableness. Control returns 
to point. 

STEP (3) PERFORMED BEFORE EACH POINT OF A LINE IS PROCESSED 
No action taken 

STEP (4) PERFORMED AFTER EACH LINE IN THE FILE IS PROCESSED 

The average is first calculated for pulse one. If this is the first 
line of the file, the program returns. This allows the program to calculate the 
relative height of pulse 2 in the following line. If this is not the first 
line of the file, the line counter is incremented; the middle reselm is determined 
for both pulses (if option 2) ; the average height of pulse 1 and the relative 
height of pulse 2 is calculated (if option 1). The mean and standard deviation of the 
height for pulse one and pulse two (if nonzero) are saved for later plotting 
(if option 3). Control returns to uotLnt. 

STEP (5) PERFORMED AFTER AN ENTIRE FILE HAS BEEN APPROVED 

If the first option is specified, the mean and standard deviation are 
printed out of pulse I’s height and pulse 2 f s relative height. 

If the second option is specified, the mean middle reselm number and 
their standard deviations are printed for pulse 1 and pulse 2 (if requested) . 

If the third option is specified, the point plot image is set up and the 
mean and standard deviation of the pulse l f s height is plotted versus its line 
number. The same is performed for pulse 2 T s height relative to pulse 1. The 
pulse height (overall mean pulse height +25) is plotted along the Y-axis* 

The line numbers are plotted along the X axis. One hundred lines are plotted 
per page and will continue for as many pages as required to plot all the lines 
for which pulse height calculations were made. (this may be fewer than the maximum 
number of lines declared in step (1). Control returns to point. 
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AVEDEV. 

Given an array of points and the number of points to be referenced, this 
routine calculates the mean and standard deviation of those points. 

P TP ROC. 

This routine saves the value of each reselm that is within the specified 
boundaries of each pulse and within the specified channel. These values 
are stored in an array according to the point 1 s position relative to the beginning 
of the pulse and according to the pulse number (2 possible). Pulse one is' used 
to dark-level -correct pulse two when option 1 or 3 has been requested. That 
is, the average level of pulse one in the previous line is subtracted from each 
point value of pulse two before the pulse two points are stored for later 
processing. 


BOUNDS. 

This routine determines the middle of any given pulse using the method 
specified in the introduction to this memo. 

WHEN 

SPECIFIED INPUT VARIABLES DEFAULT/MODE DESCRIPTION 


STEPU) BEGIN(l). . .BEGINC2) 0 

INTEGER 


END(l) . ♦ . END (2) 0 

INTEGER 

PTS 50 

INTEGER 


This is an array of 

beginning reselm numbers for each of 2 
possible pulses. These numbers must not 
be less than the beginning reselm number 
for the line. 

This is an array of ending reselm numbers 
for each of 2 possible pulses. These 
numbers must not be greater than t^e ending 
reselm number for the line. 

This is the maximum number of res elms in 
either of the pulses. This number must not 
be less than 1. 
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ir 


it 

STEP (2) 


STEP (2) 


t? 


TYPE 


LINMAX 


REJECT 

MPXMUX 


MISSON 


FDATE(l) . * . FDATE(3) 


$N0RMAL$ 

ALPHABETIC 


100 

INTEGER 


$YES$ 

ALPHANUMERIC 

NO DEFAULT 
INTEGER 


0 

INTEGER 


COMMNT (1) ♦ . .COMMNT(8) BLANKS 
ALPHANUMERIC 


This variable is used to specify which 
option the programmer wants to use. 

$N0RMAL$ - Option 1 
$JITTER$ - Option 2 
$LINPLT$ - Option 3 

This variable specifies the maximum number 
of lines that will be processed in any of 
the files between this card and the next 
MODE 1= $START$ or 
MODE 1- $RETURN$ card. 

This variable must not be less than one. 

This variable controls the listing 
of Points outside ±25 of the overall 
mean of the plot points . 

This variable is assigned the multiplexor 
channel number to be referenced for processing. 
This number must be less than or equal to 13 
and greater than or equal to 1. 

This is the mission number that is to be 
associated with this file. This number 
must not be less than zero. 

This is the date (in month, day, year form) 
that is to be associated with this file. 
FBATE(l) must not be less than 0 or greater 
than 12. 

FDATE(2) must not be less than 0 or greater 
than 31. 

FDATE(3) must not be less than 0 or greater 
than 99. 

This is any special comment that is to be 
associated with this file. No error checking 
is performed. 
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OUTPUT 


Pulse 1 


Pulse 2 


When TYPE = $NORMAL$ 

MEAN, DEV, MIN, AND MAX OF DARK LEVEL CORRECTED PULSE 


DATE = 

xx-xx-xx 



MISSION - 

xxxxxx 

"comment*' 

|'mean(1)- 

xxxxxxx.xx 



' DEV (1)= 

xxxxxxx.xx 



l MIN(1)~ 

xxxxxxx 

LINE 

NUMBER = YYYYYY 

f 

V_MAX (1 ) = 

xxxxxxx 

LINE 

NUMBER = YYYYYY 

fMEAN(2) = 
j 

xxxxxxx.xx 



{ DEV (2)= 

xxxxxxx.xx 



MIN(2)= 

xxxxxxx 

LINE 

NUMBER = YYYYYY 

1 

l_MAX(2) = 

xxxxxxx 

LINE 

NUMBER = YYYYYY 


The reader will notice that the minimum and maximum mean values are printed as 
well as the actual line number within the file where they occurred. 

When TYPE = $JITTER$ 

MEAN, DEV, MIN, AND MAX OF MIDDLE RESELM VALUE 


DATE = 

XX-XX-XX 




MISSION = 

XXXXXX 

"special comment" 

MEAN (1) = 

XXXXXXX.XX 




DEV (1)= 

xxxxxxx.xx 




MIN (1)= 

xxxxxxx 

LINE 

NUMBER 

= YYYYYY 

MAX (1)= 

xxxxxxx 

LINE 

NUMBER 

= YYYYYY 

MEAN (2) 

xxxxxxx.xx 




DEV (2) 

xxxxxxx.xx 




MIN (2) 

xxxxxxx 

LINE 

NUMBER 

= YYYYYY 

MAX (2) 

xxxxxxx 

LINE 

NUMBER 

- YYYYYY 
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In this case, the means and deviations are in terms of the middle reselm 
number for the respective pulses. In addition to the minimum and maximum 
reselm numbers, the line number location within the file where the extremes 
were found is also printed. 

When TYPE = $LINPLT$ 

First the means of pulse one are plotted versus their respective line numbers. 
One hundred lines are plotted per page; the plots will be continued on to as many 
pages as are necessary to plot the actual number of lines processed. This plot is 
followed by the overall mean, and deviation of all the lines as well as the 
minimum and maximum values with line number locations. If REJECT =$YES$ 
all those points which were not within +25 of the overall means are tabulated. 

This procedure is repeated for the standard deviation of pulse 1, the mean 
of pulse 2, and the standard deviation of pulse 2. Wo minimum and maximum 
values are printed for the standard deviation of either pulse. 

DATA CARDS SETUP 

1 . BEGIN(l) = 10,100, END(l) = 90,190, PTS - 91, 

TYPE = $LINPLT$, LINMAX - 500* 

2. INBIN * 700, UNIT = 4, FILE =5, NSA = 100,489,1,1,200,1 * 

3. MPXMUX = 5 , SCRAN = 7. MISS0N ~ 12, 

FDATE'fl) = 10,29,73, COMMNT(l) = $ABSURD$* 

4. FILE = 5, NSA = 100,599,1,1,200,1 * 

5. MPXMUX *= 7, SCHAN = 8, MISSON “12. 

FDATE(l) “ 10,29,73, C0MMNT(1) = $STRANGE$* 

6. M0DE1 = $START$ * 

7. BEGIN (1) = 90,10, END(l) ° 190,9, PTS - 91, 

TYPE - $JITTER$, LINMAX = 400 

8. FILE = 5, NSA = 100,499,1,1,200,1 * 

9. MPXMUX = 7, SCHAN = 8, MISSON = 12, 

FDATE(l) ^ 10,29,73, COMMNT(l) = $BR0KEN$ * 

10. M0DE1 * $START$ * 
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11. BEGIN (1) - 10,70, END(l) = 59,85, * 

12. TILE - 7, NSA - 10,109,1,1,86,1 * 

13. MPXKUX = 4, SCKAN = 8* 

14. M0DE1 = $RETURN$ * 

1-3 These cards specify that a point plot Is to be made of a maximum of 500 
lines under the control of option three. Multiplexor channel 5 (representing 
spectral channel 7) is to be referenced. A mission number of 12 and the date 
10/29/73 is to be associated with the file. Although a maximum of 500 lines 
are provided for, only 390 lines will actually be plotted since those are the 
only ones processed. 

4-5 These cards are a continuation of the previous processing. Therefore a 
point plot is still to be made for a maximum of 500 lines with the same pulse 
boundaries as were used before*. This time multiplexor channel 7 

(representing spectral channel 8) Is to be used. This time, the full 500 lines 
will be plotted. 

5_9 These cards specify that option 2 is to be performed on pulse 1 only (pulse 
2 had zero width) using multiplexor channel 7 of file 5. A maximum of 400 lines 
will be and are processed. 

10-13 These cards specify that option 1 is to be performed on pulse 1 and 2 using 
multiplexor channel 4 af file 7. The maximum width of either pulse defaults to 
50; the maximum number of lines to be processed defaults to 100. The mission 
number and date default to zero, and there is no special comment. 

14. This returns contool to the calling program. 

ERROR COMMENTS 

ERR0R1 . MPXMUX .G. QNCHAN 

The requested multiplexor channel number Is greater than the number of channels 
in the file. 

ERR0R2. PTS.L. 1 

This indicates that the maximum number of points in either pulse is less 
than 1. 
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ERR0R3. INVALID TYPE* 

The variable TYPE was given some value other than $NQRMAL$, $JITTER$ I or 
$LINPLT$. 


ERR0R4 . LINHAX.L.l. 

The maximum number of lines to be processed is less than 1. 


ERRORS . PTS.L. PULSE WIDTH .L.D. 

This indicates that the width of one of the pulses is either less than zero 
or greater than the maximum number of points allowed in either pulse. The number 
of points in a pulse 1 b 

END(I) - BEGIN (I) +1. I = 1 or 2 

ERR0R6 ♦ NO DARK AREA SPECIFIED. 

The width of pulse 1 ie zero. 

ERROR 7 * INSUFFICIENT STORAGE 

The particular combination of input values requires more storage for processing 
than is available. This problem can be solved by reducing LINMAX and/or PTS. 
ERR0R8. 13-L.MPXMUX.L. 1. 

This indicates that the multiplexor channel number is either less than 1 
or greater than 13. 

ERR0R9 . 24. L. SCHAN.L.l, 

This indicates that the spectral channel number is either less than 1 or 
greater than 24. 

ERROR10, MISSON .L.O. 

The variable MISSON has been assigned a value less than zero. 

ERR0R11. INVALID MONTH 

FDATE(l) has been given a value either greater than 12 or less than zero. 
ERR0R12. INVALID DAY 

FDATE (2) has been given a value either greater than 31 or less than zero. 
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ERR0R13. INVALID YEAR 

FDATE(3) has been given a value either greater than 99 or less than zero. 
ERR0R1 4. BEGIN (X) . L . NA. 

This indicates that pulse "X" has a reselm number lower than the beginning 
reselm number for the line. 

ERR0R15. END (X) .6. NB. 

This indicates that pulse "X" has a reselm number greater than the ending 
reselm number for the line. 

ERR0R16. LINES. G. L INMAX 

The program attempted to process more lines than the declared maximum. 
LINMAX may have been assigned an incorrect value or the NSA field of the PROCESS 
input card specified more lines than were supposed to be processed. 

SUBROUTINES NEEDED 

The following are a list of external subroutines needed for the correct 
operation of the program. 

ASSIGN. 

ERROR. 

LINK. 

PL0T2 . 

PL0T2 . 

PL0T3. 

PL0T4. 

POINT. 

PROCESS. 

SQRT. 

ZERO. 


This routine performs all of the dynamic storage 
allocation necessary in the A core box. 

System error processing routine 

This routine provides POINT with a connection to 
DGNSTC.'s point-processing routine PTPROC. 

All four of these routines are used in setting up and 
printing the point plots of DGNSTC option 3. 

This is the program that repetitively calls DGNSTC. with 
different levels to progressively process the data 
(see POINT writeup). 

This routine is called by POINT, and handles all I/O with 
data tapes (in ERIM format). 

This routine calculates square roots. 

This routine inserts either a binary or floating point 
zero in all of the requested variables or arrays. 
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In addition to these external functions, a define package must also be 
included. This package inserts 3 new operators into the MAD compiles and defines 
what machine instructions they each represent. These three operators are .LI., 
.LIF., & .SI. and are called the indirection operators. They are used in this 
program in conjunction with the dynamic storage allocation routines for accessing 
dynamically allocated arrays without using the usual subscription routines provided 

by the MAD compiler. Instead, they use the two internal functions S, & S22 
defined in the program for all subscription routines. 

S. is used for accessing a 2-dimensional array with maximum dimensions of 
(X, LL/B). LL/B - PTS*LINMAX (OPTION 1) or LINMAX (OPTION 2 or 3 \ X»2 when 

TYpE»$NORMALS$, or TYPE = $JITTER$. X * 4 when type= $LINPLT$, 

S22.is used for accessing a 2-dimensional array with maximum dimensions 
of (2 ,PTS) . (See writeup on Indirection operators). 

STORAGE REQUIREMENTS 

This program requires 2744 decimal locations before any dynamic storage 
allocation takes place. The following equations describe how much additional 
storage will be required for any given set of input cards. 

PTS*2 + LL + 2 

LL = PTS*LINM AX * 2 When type = $ NORMAL $ 

LL. = LINMAX * 2 when type = $JITTER$ 

LL = LINMAX * 4 when type - $LINPLT$ 
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MAC ( G9 AUG 1965 VERSION! PROGRAM LISTING 


EXTERNAL FUNCTION DGNSTC. 

REFERENCES ON 

* DEFINE PACKAGE FCR LOAD INC 1 RECT I . L I . , . U F . I ANC 
STORE INDIRECT!. SI.). THE VARIABLES BEING REFERENCED 
ML ST BE IN THE A COKE BOX. IF THt ACCP-ESS BEING REFERENCEC 
IS IN THE 0 CORE 8CX, THE VARIABLE CONTAINING THE ACTRESS 
NL ST HAVE A MODE OF 5. 

DEFINE UNARY OPERATOR .LI., PRECEDENCE SAME AS .ABS. 

MODE STRLCTLRE 1 = .LI. 1 


JMP 

*+5,AC,*+l 

JMP 

>+6»rc f *+r 

JMP 

*+7 f LA,*+l 

CLA* 

H 

LIT 

AC 

S ft 

T 

JKP 

♦ -3 

““ v 

T 

JMP 


$L k 

T 

JMP 

*-7 

tNC 

_ 


MODE STRLCTLRE C = .LI. 0, SAME SECUENCE AS .LI. 1 
"DEF I Nt" LN'A'R V '0'P!ER'ATGR'“VLT F.T PRECEDENCE SAME AS .ABS. 
KCL’t STRUCTURE 0 = .LIF. 1, S AM*E SECUENCE AS .LI. 1 
MODE STRUCTURE 1 = .LIF. C, SAME SECUENCE AS .LI. 1 
ObriNE BINARY OPERATOR .SI., PRECEDENCE SAME AS - 
MODE STRUCTURE 1 = 1 .SI. 1 


JMP 

<=+t,BT,* + A 

JMP 

<>+'3, AC, *<1 

JMP 

*+9,MC,*+l 

JMP 

*+lC»LA,*+3 

JMP 

* + 1 1 t AC ,* + l 

JMP 

* + J 2 ? L A * * + 1 

CL A 

B 

... STG *~- 

A 

CL T 

/ 

ST C* 

A 

JMP 

*-2 

STG* 

A 

JMP 

*'4 

SL V.*‘ 

A 

JMP 

£ — 6 

STC 

T - 

JMP 

*-ic 

SLV> 

T 

JMP 

*-12 

ENC 

— 


, MOCE STRUCTURE C = C .SI. 0, SAME 5ECUENCE AS 1 .SI. I 

MODE STRUCTURE 1 = 1 .SI. 0, SAME SECUENCE AS 1 .SI. 1 

KOCt STRUCILRE C = C .SI. 1, SAME SECUENCE AS I .SI. I 

NORMAL MODE IS INTEGER 

HL CAT. INC. POINT FI, ME AN , 06 V ,TCT , SCRT . , F C A TUM ( 2 A I , MX 1 , E E I C ) , X ( ICC ) 

DIMENSION BEGINI2) * END (2 ) »FARAMS(2 ) , FT (2! ,FCATE( 3) , C OMMNT ( 8 I , 

1 MIM2),MAX(2),NSCALE(A),LINMIM2),LINMAX(2),UT(3) 

ECU l VALENCE (FDA TUM, DATUM) , ( PNT L » F AM AMS I 1 ) ) , 

1 (PMF.PARAMSI2) ) , ( E t( 0 I , JUNK ( 0 ) ) , 

2 UFLAG,FDATE) , (PT.LCNT) , 

1 INC, END). , (MS, BEGIN) . 

VECTOR VALLES ERR(CLl) = IH 1 X * + * ** ERROR 2 . PTS .L. !.•* 
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VECTOR VALLES ERRI0A3I = $H» L** * ** ER9CR3 . INVALJC TYPE.'* I 

VCC TOP, VALUES ERR C 0 H ■> ) = f H * 1*** ** ERRCR4 . LINMAX .L. 1.'* t 

VECTOR VALUES ERRI127) = $H • 1*****£RKCR5. PTS.L. PULSE WlDTF .L.C.'*i 

VECTOR VALLES EKRCIF.9) = tH *1 * * * ** ERRCR6 . NO DARK AREA SPEC IF I EC. • * * 

VECTOR VALLES ERR1211) = *H» 1 * ** ** E RRCR 1 . INSUFFICIENT STORAGE. • * l 

VECTOR VALLES ERR ( 2 5 3 ) = $H»1 ** * ** ERRORS . 13 .L. MPXMIJX .L. l.»* S 

VECTOR VALUES ERR 1 2 3 ) = Hi' 1 ** * ** t PRCR9 . 24 .L. SCHAN .L. 1.** i 

VECTOR VALLES ERR 1337) = tH • 1 ***** tKKGK 10 . MISSCN .L. O.'* 5 

VECTOR VALLES EHRI379I = th • 1 ***** ERROR 1 1 . INVALID MONTE.'* t 

VECTOR VALUES ERR1421) = tH • 1 ***** ERROR 12 . INVALID DAY.'* t 

VEDIOR VALLES ERR (463) = $H ' L*****ERRCR 1 3 . INVALID YEAR.'* $ 

VECTOk VALLES ERR I 505 ) - tH • l** *** ERROR LA . BEG I N { • , 1 1 ,H • ) .L. NA.'* * 

VECTOR VALLES ERR15A7) = JH'l ***** ERRCR IS . £NC( • , 11, H. • ) .G. NP.»* * 

VECTOR VALLES r»R<*>e < 9» = tH • 1 ***** ERROR ID . LINES .G. LINN AX.'* t 

VECTOR VALLES ERRI631) = tH • l ** * ** ERRCR 1 . MPXMUX .0. CNCHAN . • * % 


VECTOR VALLES FSGIO) = tMEAN LEVEL i 
VECTOR VALLES M SG t 2 ) = tOEvlATICN LEVEL* 

VECTOR VALLES NCHAR * 10,0,15 
VECTOR VALLES PCD = t*l 

VECTOR VALLES TIT = tS T ANDARD CE V t » tME AN $ 

STATEMENT L/.ITEL STEP I 7 ) 

ERASABLE JUNK ( 225 ) , CUN I T . CREEL ,CF I EE . CLINE ,TWR I T E . CICH AN , 

1 CNWDS,PKCG,MU0E1 .MOOE2 ,IJMT , C AL I N E , CDL I NE , NS A, NSB . 

~"kS,NA,M 3,KP»[0(1) ,TFLAG, l PACK,RESERV< A ) ,CFACTR<49 ), 

0 T 1 TL2I1S) ,CTI T LEI 19) .CUSTt 19) , CS P AR E I 46 ) , CFL AC , 

Cl' AWG «CBANG , CRP » CNA , LNSS , CNCH AN , CMCCE , CR EC A, QR ELC , 

OF 1LE, CREEL, CNViRCS 

ERASABLE DATA(4?3),ITEST,CSTAKT,CSTART,NV,NX,NC,L,IP, 

1 Ii:P,NCP,NEXT,START2,READT2,R6ADL2»REACP2,LAHEL2, 

2 ■' EXT'<A (2C) ,UATUM2A) , ICCCK2A ), ICHAM2A ), IMAGE! 1040) 
INTERNAL FUNCTION S22.(Si,S2) = PTS*{Sl-l) + S2 

INTERNAL FUNCTION S.(S3,S4) =LL/B*(S3-1) + S4 
VHFNFVEK NEXT .G. 5, FUNCTION RETURN 
TRANSFER TC STEP (NEXT) 

STEP ( 1 ) L I NK . ( P TPRCC • ) 

L =" L + nCHAN 
LINMAX = ICC 
P IS = SC 
TYPE = tNCKMAL* 

' RFJECT = tYESS 

(R = l,l,K.r,.a,CCRMNT(K) = t $) 

ZfcKU. (BEG I N. . .BEG IN 12 ) »EN'C. . . bNCi2 ) ,P ARAMS. . .PAKAMSI 2 ) » 

1 N SCALE II 1...NSCALEI3)) 

READ AND PRINT DATA 8EG1MU, ENO I 1 ) , PTS, TYPE, LINMAX 
M SCALE = 1 
NSCALE ( 2 ) = -1 
N SC ALL I A ) = -1 

nFiENEVER' PTS .Li 1 

TRANSFER TC ERR2 

CR WHENEVER TYPE .NE. tNORMALt .AND. TYPE .NE. tJITTERt .AND. 

I TYPE .NE. tLINPLT* 

TRANSFER TO E«R3 

OR WHENEVER LINMAX .L. 1 . 

, ... _ TRANSFER TO ERRa’ 

END OF CONDITICNAL 

THROUGH Ll’l , FOR 1 = 1,1, '"i 7G. 2 

WHENEVER cNOin-HEGINl IJ+l .L. 0 .CR. END ( I J-BEG INII1 ♦ 1 .G. PTS, 
l TRANSFER TC ERRS 

LPl CONTINUE . 

E) = 2 

FI = ENDIl) - BEG INI 1) + 1 
WHENEVER FI .L. 1. .AND. 

1 (TYPE .E. tNCRMAL* .OR. TYPE .E. tLINPLTt), TRANSFER TO ERR6 
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LL“ PTS*lIt\/ , 'AX*2 

WHENEVER TYPE ,F- $JITT£K$» LL= LiNMX*2 

~ WHE NEVE RETYPE *LI‘NPLT$ 

LL = LINMAX*4 
It = A 

END UF CONDITIONAL 
KK = ASSIGN. (C,C) 

WHENEVER KK , L . LL+ PTS*2 + 2» TRANSFER TC ERR7 

KK * ASSIGN, (0,PARAPStPTS*2tlL) 

FUNCTION RETURN 

STEP! 2 > ZikC. CP T • ,,PT (2 > *FDATF.. . .FCAFE<3) , F L AG t L I N E S , M I S SCN ) 

SCI- AN - 24 

HEAD AND PRINT DATA MPXMUX, SCHAN, MlSSCN,FCATEm,COMKNT Cl), REJECT 

WHENEVER MPXMLX - L - i .CR. MPXMUX .G * 13 

TRANSFER' TC EKK8 

Ck WHENEVER SCHAN . L. 1 -CR. SCHAN .G. 24 
TRANSFER TC E R 
OR WHENEVER M l S SON ,1* 0 
TRANSFER TC'FHRLO 

GR WHENEVER FDATECl) . L. 0 ,CR. FGATC(l) »G. 12 

CK WHENEVER F 0 A T fc ( 2 ) .1. 0 .OR. K4HM2) *t. 3t 
■ ’ TRANSFER TC ERR12 

cr< WHENEVER F 0 A T C ( 3 > ,L. 0 .Lrt* FUAT E ( 3 ) .G. 99 
TRANSFER TC FRR13 
CK WHENEVER MPXMUX ,G* CNLHAN 

• - - ’ TRANSFER TC CRkl 

CNU CJF CONDITIONAL 
MIN (1) = 9S9 

HM2) - 999 

f'AKU) r ~S99 

f i A X 1 2 ) = -9 9 9 _ _ ...... 

• * ’ THROUGH L P 2 , FOR T =1 , I • T • G • 2 

WHENEVER (3 E G I iN { l ) , L • N A 
TRANSFER TC ERR14 
OR WHENEVER FNU (I) * G . NO 

TRANSFER TC ERR 1 5 
END OF CONDITIONAL 

LP2 “CONTINUE 

FUNCTION RETURN 
STEP ( 3 ) FUNCTION RETURN 

STEP<4) IFLAG = l 

EXECUTE AVEDEV . C JUNK ( 7 7777K-PNTL + B22 , ( 1 , 0 I ) , 

1 ENDtD-BEGINU J+l) 

-J = >' F A N + ,5 

WHENEVER FLAG ,E. 0 
FLAG = i 
FUNCTION RETURN 
fc\C OF CONDI TICK AL 
L INFS •= LINES + 1 

"WH.LNEVER LINES ,G. LlNMAXf TRANSFER TC ERR16 
THROUGH I.P3, FCR 1=1,1* I ,G. 2 

WHENEVER ENnm - BEGINU) + 1 ,le. o, transfer TO UP 3 
WHENEVER TYPE .£. $J1TTFR$ 

EXECUTE BOUNDS, 

10NTF-S, ( I, LINES)) .SI* MIDDLE 

WHENEVER MIDDLE + BEGIN! i ) - 1 ,L, MIMI) 

HMD = MIDDLE 4 BFG IN ( I) - l 
LlKMIM l ) = LINES + NSA - 1 
END CF CONDITIONAL 

WHENEVER NICOLE + BEGIMI) - 1 ,G. MAXU) 

MAX (I) = MIDDLE + BEGIN! I) - 1 
L INMAXI l J = LINES + NSA - 1 
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SMI CF CONDITIONAL 
TRANSFER TO L F3 
OR WHENEVER TYPE .E. tNORFAL* 

IK=1 ,1,K.G. lENOll I-0EGINI I>*1>» PT(I> = Pill) + l, 

IPNTF-S.ILrPTIll)) .SI. . L I -. (PNTL-S22.I I » K ) I ) 

TNO OF' CCND I T I CNAL 

I FLAG = l 

hH r NFVER l .E. 2, EXECUTE A VEOE V . ( JUNK ( 77 77 7K-PNTL + S2 2 . I 2 , C )) * 

ENC12)-BEGIM2)+1) 

WHtNEVFK FEAN .L. FIMI) 

HINT I ) = FEAN + .5 
' LINFIMI ) = LINES + NSA - 1 
EMJ OF CONDITIONAL 
WHENEVER FEAN . 0. FAXIII 
FAX ( 1) = FFAN ♦ .5 
LINFAX(l) = LINES ♦ NSA - 1 
END OF CONDITIONAL 
UFiENL XlirK TYPE .E. FLlNPLTi 

IPNTF-S. (1*2-1, LINES)) -SI. FEAN 
IPNTF-S. ( 1*2, LINES) > .SI. DEV 

iM 


LP3 CCMiNUE 

FUNCTION RETURN 

STEP(E) WHENEVER TYPE . E . INQRFAL* 

PRINT FGRFAT $H • 1FEAN , CEV, FIN, AND FAX CF C ARK LEVEL CORRECTED 1 , 
I H* PULSE.'///** 

CK WHENEVER TYPE .E. tJITTEKt 

PRINT FCRFAT IH’IPEAN, CEV, FIN, AND FAX OF FIICOLE RESELF VALUE.'/ 


1 //** 

END UP CONDITIONAL 

PRINT FCRFAT FH 1 -CATE = • , 1 2 , S 1 , 12, SI, I2,S10,H* FISSION * ',16, SIC, 

1 HC fc* i * FCATfc(l).. . FDATE 13) ,FISSCN» 

2 CCFFNT { 1 ) . • .COFFNT 1 8 I 

WHENEVER TYPE «E. JLINPLTt, TRANSFER TC LINPLT 
T HRUUGH LPb, FOR 1=1,1, I .G. 2 

• WHENEVER END ( 1 1 —13 CO I N ( I I + 1 ..LE. 0, TRANSFER TC- LP5 

WHENEVER TYPE .E. iNORFAL* 

EXECUTE AVECEV. (JUNK! I77/7K-PNTF+S. ( 1,01 ),PT( I ) I 
OTHERkl SE 

EXECUTE AVFDtV. I JUNK I77777K-PNT F+S. II ,0 ) ), LINES) 

FEAN = FEAN + BEGIN I I ) - 1- 

F.NO UF ClnOI TICNAL ' ” 

PRINT FCRFAT pH 1 PE AN I ' » 1 1 ,H ' I = • , S 10 , F 10 . 2/H • DEVt'.Il.F') = ', 

1 S1C,F1C.2/ 

I H* FINISH, H'J = SSIO.I 7,S5,H*LINE NUFBER = ',16/ 

1 h' F AN t ' , I L ,H* I = SSIG , 17 ,Sb,h 1 L IKE NUFeER - ',16*1. 

2_ _I,FEAN,t,DEV,I,FIN|I),LINFIN(I),l,FAX(I) , LINHAX II) 

LPS “ 'CONfiNCe 

LINPLT WHENEVER TYPE .P . iLINPLTt 

EXECUTE PLOT 1 • IN SCALE, 5, 10, 5, 20) 

I FLAG = C 

THROUGH LP1C, FCR 1=1,1, I .G. A 

PS = 2 • 

WHENEVER 1/2*2 .ME. I , PS = C 

J = I 1 + 1 )/Z 

WHENEVER END! J) - BEGINIJ) + l .LE. 0, TRANSFER TC LP10 
EXECUTE AVECEV. (0, LINES) 

THUCLGH LPU, FCR LCN f = 10l , ICG, LCN'T .GE. ILINES+101) 

(K = ILCM-10D) ,l,K.G. (LCNT-1) , X I K-LCNT+ 10 1 ) = K> 

FXECLTE PLOT 2. ( l PAGE » LCNT- 1 » , LCM- 101 . » F £ AN + 25. * 

1 FEAN-25.) 


NO = ICC 

WHENEVER LINES .L. LCNT+1, NC = LINES - LCNT + 101 
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LP11 


IP4 

LP1C 

BLfcKl 
ERR 1 
ERR 2 
ERR 3 
ERK4 
£Rk5 
ERP.6 
ERR f 
ERR B 
ERR9 
ERR 10 
tRRll 
ERR 1 2 
ERR 13 
L R R 1 4 
ERR 15 


1 

1 


i 


1 

2 


tXECtTf PLCT3. IHCC.Xtl) .EE (77777K - PNTf- ♦ 

S. ( 1 ,ICM-L00> ) ,NC) 

~ PKlM FC'PPAT tlHl .552.H* PULSE 1 , I 1 , S 1 » 2C6* i . J * 
1 IT I (1-1/2*21*2), TIT <<I-I/?*2M2+1) 

FXECLTF PLCTA. (IvCHAMMS) , PSG(KS) ) 

PRIM FCKPAT tlHS . S55 ,K 1 L 1 A E KUKEER'** 


CENT INLE 

PRINT FCPP6! iH'Lf-tiN = ■ * S5 . F 10. 2/ /h ' EEV = 
‘ ' PhAKtUtV 

WHFNF'VFR 1/2*2 .K£. 1 

PRIM FCUf'AT *H* M I N = * » St> * I 7 tS5*H’HNE 
H> MAX = • ,S5,17,S5,MLII\6 NUMBER = 
MNtJ) » L 1 NK 1 N( J ) iPAXIJ)*LINi w AXIJ) 


1 ,S5,F1C.2//*S» 


NUMEER = '.It// 

• » 16 *$. 


FNH CF C CM) I T 1CNAL 

V»HhKV Vt:K REJECT . £ . tNCt, TRANSFER TC LP10 


K K = 0 

PRIM FORMAT $lHliS50tH'PClNT$ CUT L F RANCEVIH « S A 5 , H • L I ME 1 t 
1 S6pH f PEAK 1 , SS.H^STMIRU DEV' //** 

THkCUUH LP4 ? FCR LCM = lfl» LCM . C * LINES 

rtHFWEVtR EFC/7777K-PNTP + S. UfLCM) ) .L. PFAN+25. • A N K . 

y~ Cl : I77777K-PNTF+S. I I'.LCNT ) ) .0, RE AN-25 • ♦ TRANSFER TG LPA 

K K = KK ♦ l 

PRINT FORMAT $IH « S44 ♦ I 6 * S2 • F 1 0 « 2 i S 2 * F 1 0 * 2 * % * LCNT t 

1 Ec 177777 K-PNTP+S- < J *2 -It LCNT) ) » 

2 EE(77777K-PNTF + S.(J*2tLCM)> 

CONTINUE 

pr<j NTT'CRKAT "$H 1 -TOTAL "'NUMBER' CF POINTS OUT OF RANGE - KK 


CONTINUE 

FNB or CCNDITIGNAL 
PR [NT r GR t* A T % LH1*$ 
FUNCTION RETURN 
PRINT FORMAT EKKOOtl 

TRRGR 

K = 6 3 i 

TRANSFER TC BLFH 
K = 1 


TRANSFER TC BLFH 


l< = 4 3 

TRANSFER TC BLEH 
K = B 5 

TRANSFER TC RLEH 
K * 127 

TRANSFER TC BL EH 
K = 169 

TK AN SF ER T G BLEH 
K * 211 

TRANSFER TC BLEH 
K = 253 

TRANSFER TC BLFH 
K = 295 

"TRANSFER TC BLEH 


K a 33 7 

TRANSFER TC BLFH 
K = 379 

TRANSFER TC BLFH 


K - 4 7 1 __ 

TF ANSI LK TC BLEH " 

K - Ah3 

TRANSFER TC BLEH 
K - 5 T 5 

TRANSFER TC BLtHl 
K = 547 

fk AN S PE iV TC BLEH I 
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ERR 16 
BLEH 


LP6 


LP 7 

LPH 
DUU 
LP 1 2 


K = 609 

transfer to hlfh 

PRINT FCRFAT ERRIK) 


INTERNAL FUNCTION AVECEV- (FAT t P MS ) 

WHENEVER IFLAG .E. L 

HHAN = { TOT * C* iK = l 1 1 * K . G. PNTS» TGT + PNT<K>) / PNTS 
JOT * nCT«0,,K»lil.K f G.PMS* TCT 4 [ PNT ( K ) -WE AM . P . 2 ) 

0 h V s S tfi V. ( 7 C T / P N T S > 

CTHERWi SE 

FFAN ^ ( T C T = C * |K-I tl »K.G.PNTS# 
l TCT + .LIF. ( PNTF-$. U,K) H/PNTS 

TUt = ITCT-C.*K-l*l#K.G.PNTS* 

1 TL T + C • L I F • IFMF-S. tI,K) ) - FEAN). P.2) 

DEV = SORT. ( TCT /PNTS ) 

END OF CONDITIONAL 
FUNCTION RETURN 
END CF FUNCTION 


INTERNAL FUNCTION PTPROC. 

THROUGH L P 6 . FCR [*1 ♦ 1,1 -G. 2 

WHENEVER IP’.LE- PM) III .ANO. IP . GE • BEGINU) 

WHENEVER CFCDE . E • i 

K ■= FDATUfMPPXPUX) 

D ATU M ( R PXHUX ) = K 
EM) CF C CMJ I T I CN AL 

WHENEVER TYPE • E * i J 1 T T E R $ .CR. II .E. l • ANJJ ♦ 

1 (TYPE .E. $N OR MAL$ -OR- TYPE .6. iLINPLTS)) 

(PM L-S22.II,lIP-f3EOIN(l) + l))) .SI. C ATUM { MPXtfUX ) 

OTFtRMSt " 

(PNTL-S22.il v (IP-bEGlN(I )+l))) .SU 
1 DA TUP I FPXFUX ) - J 

END CF CONDITIONAL 


END OF CCNDITICNAL 
CONTINUE 
F LNCtl CN ’ RE TO RN 
ENG OF FUNCTION 


INTERNAL' FUNCTION BOUNDS. 

MXl = 0. 

THROUGH LP 7 1 FCR KK=l t l* KK «G. (EKQ(T)-REGIN(I)+l) 

WHENEVER . LI • ( PN f L-S22. U «KK )) *G. RX1, 
l MXl = .U . (PNTL-S22. (I'tKKl ) 

CONTINUE 

L C F T = f ~ ' 

THROUGH LP12, FCR FI = .5,.l, LEFT . NE . 1 

THROUGH LPH, FCR KK = l » L *KK .G. UNC < f l-EEG t N< I ) + l > 

WHENEVER F X 1 *F 1 .LE. • L I • ( PML-S 22 . < I f KK 1 ) « TRANSFER TO 0 UT 1 
CONTINUE 
LEFT = KK 
CONTINUE 

THROUGH LP9, FCR KK- ( END < ll-BEG IN I I) ♦ l > *- 1 1 KK .L. 1 

WHENEVER FX1*F1 .LG. .LI- (PNTL-S22.I I»KK) ) « TRANSFER TO CUT2 
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LPS 

0UT2 


COM I NUE 
lUGFT = KK 

MCULE = (LEFT + RIGHT 
" FI NCT I CN RETURN 
F.fvR OF FUNCTION 
END OF FUNCTION 


+ l )/2 
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Appendix C 

RADIANCE CALIBRATION 

The scanner calibration lamp signal is used basically as a transfer standard to obtain the 
apparent spectral radiance of terrain objects. Its main purpose is to account for responsivity 
changes in the scanner system. In general, these are almost impossible to monitor conven- 
iently. The scanner components themselves are designed to respond linearly with increased 
radiance from the ground terrain (i.e., an increase in radiance gives a corresponding increase 
in detector signal). Hence, the scanner detector voltage for each channel can be represented 
by the following general equation: 

V T+p (X) = KfL T (X) + L p W] (C.l) 

p t e t W 

L T ' 7T 

where V T+p (\) = scanner signal voltages recorded by aircraft recorder for each channel 

K = constant dependent on system factors such as responsivity of detectors and 
photomultipliers, system transmission losses or gains, etc. 

L t (\) - radiance of target (assuming a Lambertian surface) for each channel 

Lp(M = radiance produced for each channel by scattering of radiation by molecular 
and aerosol particles in the atmosphere 
P r = object directional reflectance 

At present, the quantity L^(\) is impractical to measure conveniently and simultaneously 
with airborne scanner imagery. Considerable effort is being spent on modeling the atmosphere 
so that the appropriate corrections can be made in the data during processing. 

In terms of scanner radiance calibration, however, it is not important to consider the 
individual terms of Eq, C.l, but rather the total radiance received at the scanner’s aperture. 
Hence, the apparent target radiance (which includes the contribution of path radiance) is a 
more appropriate quantity to measure. Simplifying Eq. C.l, then, we get 

V T+p = KL T+p^ (C.2) 

where = apparent radiance of target [L,p(A) + L^(X)] for each channel 

Were it not for the fact that the constant K in Eq. C.2 actually does vary somewhat be- 
cause of changes in detector and responsivity, the radiance calibration of a scanner would be 
easy. A quick look at a radiance source (such as a reflectance standard illuminated by cali- 
brated irradiance standards) would suffice. Since it does vary, however, a means of removing 
K from Eq. C.2 is necessary. 
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The prime purpose of the lamp reference source, as stated earlier, is to remove the 
factor K. This is done by dividing the voltage obtained from the target by the lamp voltage. 
Hence, if the target voltage is given by Eq. C.2 and the lamp voltage is given by 

V l = KL l (\) (C.3) 


then, dividing C.2 by C.3 we get 


V 


T+p 


V L L W X) 

l l< x > 


(C.4) 


Radiance calibration is performed in a similar manner except that the target, instead of 
being a terrain object, is an object of known radiance. If the calibration target is 3M white 
paint (a well-measured and stable reflectance standard) illuminated by quartz-iodine lamps, 
then the voltage generated by the radiance standard is 


v 3M (A) = K ' L 


A) 


where L gM (X) = 


3M V 
P 3M E QI ( ^ 


(C.5) 


E^, = spectral irradiance of quartz-iodine lamps 
vl 

P 3M “ directional reflectance of 3M white paint 
K 1 = system constant which has varied from K in previous equation 


Dividing C.5 by the lamp voltage 
V L ,(X) = K'L l (X) 


(C.6) 


we get 


3M 


V L ,P 3M E QI^ 


(C.7) 


Substituting Eq. C.7 into Eq. C.4 and rearranging, we have 


L T+p (X> 





(C-8) 


Equation C.8 gives the apparent radiance of the target in terms of known or measured 
quantities. It is obvious from Eq. C.8 that the accuracy of radiance calibration depends upon 
two factors: (1) how well and E ^ are known, and (2) L l (A) remaining spectrally constant. 
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Appendix D 
PROGRAM LAM RAD. 

The purpose of this program is to calculate a set of statistics from a set of special test 
runs on the M-7 scanner. These statistics are used to determine the radiance of the lamp 
reference source at various current levels. The test runs are made with the aircraft on the 
ground hooked up to laboratory power with a special test panel sealed underneath the scanner. 

The user must specify the reselm numbers of the beginning and ending of the dark area, 
the lamp pulse and panel pulse (whose average height is to be determined), and the lamp pulse 
and panel pulse (whose integral is to be determined). Before the lamp pulse or panel pulse is 
used, the mean of the dark level from the previous line is subtracted from each reselm of the 
lamp and panel pulses. We assume that the dark level does not change significantly from 
line to line. After this, a given number of reselms are selected at the top of the panel pulse 
and the lamp pulse. The level of these reselms is averaged over the entire set of lines re- 
quested for the lamp pulse and the panel pulse. At the same time, the integral of the entire 
lamp pulse and the panel pulse can be calculated for each line and averaged over the requested 
lines. This is, in turn, divided by a specified pulse width resulting in a number equivalent to 
the mean height. After this, a level ratio (ratio of the mean lamp level to the mean panel level) 
and an integral ratio (ratio of the mean integral of the lamp pulse adjusted by the lamp pulse 
width to the mean integral of the panel pulse adjusted by the panel pulse width) are calculated. 
These two ratios along with the mean lamp level, mean panel level, mean lamp integral, and 
mean panel integral are used to calculate the standard deviation of the respective means. 

The level ratio standard deviation is calculated by determining the level ratio for each 
line and determining the standard deviation of these level ratios from the mean level ratio for 
the file as a whole. A similar process is used for calculating the standard deviation of the 
integral ratio. The standard deviation of the lamp level and panel level is calculated by de- 
termining the square root of the average squared deviation of each of the selected reselms 
from the corresponding mean lamp level or mean panel level. A similar process is followed 
to calculate the standard deviation of the mean integral level for the lamp pulse and the panel 
pulse. This is done for up to 13 possible channels. The procedure for selecting the reselms 
at the top of the panel pulse and the lamp pulse is as follows: 
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1. The maximum height is determined [MAX]. 

2. The reselms just above MAX/2 <4*MAX/5 if this is the panel pulse) are found [X&Y] . 

If the pulse boundaries are too close or the digitizing gates were too narrow, X&Y 
may correspond to these boundaries. 

3. The level of the X&Y reselms are compared. If the difference between these levels 
is greater than 10 units, a 2nd degree curve is fitted to the data points. The curve is 
then extrapolated to determine the X&Y boundaries such that they have equal levels at 
these new boundaries. 

4. The average of X&Y is found (X + Y)/2. 

5. The reselms selected are those (TOP-l)/2 reselms on either side of (X + Y)/2. 

The aforementioned calculations are performed for each file of the given set of files re- 
quested. After all the files have been processed, plots are made of the integral ratio vs. lamp 
current and then the level ratio vs. lamp current. These two plots are followed by a table of 
all the calculated statistics. The plots and table are made for all channels that were requested. 
Although the user must specify a certain area as a dark area to be used in correcting the other 
pulses, the user may also calculate: 

a. Mean and standard deviation levels of two pulses. The number of reselms selected 
for either pulse may be different. 

b. Mean and standard deviation integrals of two pulses not necessarily related to the 
previous pulses mentioned in (a). 

NOTE: In order for this program to work correctly, all pulses must be positive -going. 
Also, each set of lamp radiance tests must have the files in ascending order of lamp currents 
and the lamp currents must be evenly spaced. 
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STEP (1) — PERFORMED 1ST TIME FUNCTION CALLED OR WHENEVER "M0DE1" = $STARTS$ 

This step initiates the maximum number of channels to be referenced and the 
maximum number of files in the set of radiance test runs. Also the minimum 
and maximum lamp currents are initialized in addition to the maximum number of 
reselms in any pulse and the maximum number of reselms to be selected at the top 
of the pulses. The beginning and ending reselm numbers of the pulses are zeroed. 

Data cards are read in (if any) that may alter the table of constants used 
in calculating the radiance ratios, the maximum number of files in a set, the 
beginning and ending reselm numbers of the pulses, the maximum number of channels 
to be referenced, the minimum and maximum lamp current values, the number of 
reselms to be selected at the top of the pulses, and the maximum number of points 
in any pulse. These values are then checked for reasonableness and the storage 
for the data base is allocated. 

At this time, as many cards are read in as there are channels to be referenced 
Each card specifies the multiplexor channel number, the spectral channel it 
represents, the aircraft channel it represents and the lamp and panel pulse 
width constants to be used with this channel. Each card is checked for 
reasonableness and then a title is printed describing the spectral channel. After 
all cards are read in, the pass number is set to l,and control returns to POINT. 

STEP (2) PERFORMED BEFORE FILE IS PROCESSED 

The totalling arrays and the point counters are zeroed. If this is the first 
pass through the file, the actual number of channels in the file is checked to 
make sure that it isn't less than the number of channels to be referenced. 

Also, the beginning and ending reselm numbers of the pulses are checked to make 
sure that they are within the limits of the lines in the file. Control then 
returns to POINT. 

STEP ( 3) PERFORMED BEFORE EACH POINT OF A SCAN LINE IS PROCESSED. 

No action taken. . 

STEP (4) PERFORMED AFTER EACH LINE IN THE FILE IS PROCESSED. 

If this is the first line, the line counter is incremented and the average 
level of the dark area is determined. No other processing is done on the first 
line. Control returns to POINT. 
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For all succeeding lines during the first pass, the line counter is incre- 
mented and the following occurs . 

1. Reselms are selected at the top of the 2nd -and 3rd pulses as 
described in the introduction of this memo. The number of points actually- 
selected is accumulated in the point counter "PNTS n for later use. 

2. The value of each reselm is accumulated for the first 3 pulses. 

3. A simpson integration is performed on the 4th and 5th pulses and 
adjusted by their respective pulse width constants. These are also accumulated 
separately. 

4. The point counter "FT" j_ s zeroed. 

5. Control returns to point. 

For all succeeding lines during the 2nd pass, the line counter is incremented 
and the following is done: 

1. The squared deviations are accumulated for the 1st through 3rd pulses. 

2. The mean is determined for the 2nd and 3rd pulses. Then the 2nd pulse 
mean is divided by the 3rd pulse mean. This ratio is then multiplied by the 
appropriate factor from the table (the level ratio) 

3. The squared deviation of this ratio from the overall mean ratio is 
accumulated. 

4. The squared deviation of the integrals are accumulated for the 5th and 
6th pulses. 

5. The adjusted integral of the 5th pulse is divided by the adjusted 
integral of the 6th pulse. This ratio is then multiplied by the appropriate 
factors from the table (the integral ratio). 

6. The squared deviation of this ratio from the overall mean ratio is 

accumulated. 

7. The point counter "PT" is zeroed. 

8. Control returns to point. 

STEP (5) PERFORMED AFTER AN ENTIRE FILE HAS BEEN PROCESSED 

If this is pass 1, 

1. The file counter is incremented and checked to make sure it hasn T t 
exceeded the maximum number of files set in STEF 1. 
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2. The totals for pulaes 1 through 3 are divided by the actual total number 
of reselms referenced (stored in the point counter ,, PNTS tr ) * 

3. The totals for pulses 4 and 5 are divided by the value of the line 
counter minus 1. 

^ mean of pulse 2 is divided by the mean of pulse 3* This ratio is 

multiplied by the appropriate factor from the table. This new number is the 
overall mean of the level ratio. 

5. The mean of pulse 4 is divided by the mean of pulse 5. This ratio is 
multiplied by the appropriate factor from the table. This new number is the 
overall mean of the integral ratio. 

6. If TYPE = $PWIDTH$, the mean of the lamp level is saved for later 
calculation of the pulse width. 

7. The pass number is set to 2 and control returns to POINT. 

If this is pass 2, 

1. The totals for pulses 1 through 3 are divided by the actual total number 
of points used (PNTS). The square root of this number is called the standard 
deviation. 

2. The totals for pulses 4 and 5 are divided by the value of the line 
counter minus 1. The square roots of these numbers are called the standard 
deviation of the integrals. 

3. The totals of the level ratio and integral ratio are divided by the line 
counter minus 1. The square roots of these numbers are called the standard 
deviation of the level ratio and the standard deviation of the integral ratio, 
respectively. 

4. The pass number is set to 1 and control returns to POINT. 

STEP (6) PERFORMED WHENEVER M0DE1 = $FINAL$ 

This step is not performed if TYPE is not equal to PWTDTH or if any of 
the pulse widths are zero. 

All previous working storage is freed except for the arrays containing 
the results calculated in step 5* New storage is allocated for use by this section 
of the program. The plotting image array is setup and the X coordinates are 
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determined* These X coordinates correspond to the lamp current and are calculated 
by dividing the interval between the minimum and maximum current levels (set in 
STEP1) into as many intervals as there are files in this set. The plotting 
routine is then called three times, with three different sets of Y coordinates 

The 1st time is with the level ratios, the 2nd time with the level ratios multi- 
plied by 10, the 3rd time with the level ratios multiplied by 100. The plot 
is then made. 

The plotting routines are then setup all over again and the. integral ratios 
for each current are used. Again the ratios are plotted 3 times. The plot is 
then made. 

Finally, a table is produced of all the results: 

1. Mean of the dark level, the lamp level, the panel level, the level ratios, 
and the integral ratio. 

2. The standard deviations of all of the above. 

Control returns to POINT. 

PTPROC, 

This routine saves the value of each reselm that is within the specified 
boundaries of each pulse. These values are stored in an array according to the 
position relative to the beginning of the pulse (the point counter), according 
to the multiplexor channel number, and according to the pulse number (3 possible). 
Pulse number one is considered the dark level and its uncorrected value is saved. 
All remaining pulses are dark-level corrected before being stored. That is, the 
average dark level of the previous line Is subtracted from each point value 
before it is stored. 

BOUNDS . 

This routine selects the points in the 2nd and 3rd pulses that are to be 
used in calculations. It calculates the reselm numbers of the left and right 
boundaries relative to the beginning of the dpulse. (The selection procedure 
is described in the introduction to this memo.) In addition, this routine 
counts the number of points actually selected for each pulse of the first 
three pulses for each channel. These numbers are used later in calculating 
the mean and standard deviation of these pulses (see STEP5). 
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ADJUST. 

routine is called by BOUNDS, whenever 3 least^squares curve fitting 
is to be done. It takes the calculated coefficients from LSQAR. and solves the 
quadratic equation for the reselm number that will result in a symmetrical pulse. 
This reselm riutnueir is then used to determine the left and right boundaries of the 
set of reselms to be selected. 

SIMPSN. 

This routine performs a Simpson integration of pulsea 4 or 5. The formula 
used is: 


Q ZZ " X 1 + 4X 2 + 2X 3 + 4X 4 + *• - + 4X n-2 + 2X n-l + X n 


n must be odd 
n > 5 


INPUT VARIABLES 
WHERE USED VARIABLE 
STEP (1) FACTOR (1) 


DEFAULT 


DESCRIPTION 


(1) 

2725. 

(13) 

1244 

.3 

table of factors that ratio of 

(2) 

133.8 

(14) 

1017 

.3 

pulse 2 to pulse 3 and ratio of pulse 

(3) 

1016.1 

(15) 

840. 

65 

4 to pulse 5 is multiplied by. 

(4) 

3484.5 

(16) 

622. 

69 

No checking is performed on these 

(5) 

1016.1 

(17) 

0 . 


numbers. Once this table is 

(6) 

2166.1 

(18) 

0 . 


changed, the default values can 

(7) 

3484.5 

(19) 

0 . 


be reinstituted only by reloading 

(8) 

3078.5 

(20) 

0 . 


the program or by reading them 

(9) 

2459.7 

(21) 

0 . 


in off input cards. 

(10) 

2111.2 

(22) 

0 . 



(ID 

1712.7 

(23) 

0 . 



(12) 

1469.6 

(24) 

0 . 




FILKAX 


6 

>1 

INTEGER 


The maximum number of files to be 
considered as one set of lamp 
radiance test runs. This number 
need not be exact (just so it is 
not fewer than the actual number 
of files to be processed) since the 
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STEP (1) BEGIN CO. . . 

BEGIN (5) 


END(l) . . . 
END (5) 


actual number of files processed is 
used to determine the lamp current 
test points. 

0 This is an array of beginning reselm numbers for 

INTEGER each of 5 possible pulses. These must not be 

less than the beginning reselm number specified 
for the line. 

0 This is an array of ending reselm numbers for each 

INTEGER of 5 possible pulses. These must not be greater 

then the ending reselm number specified for the 
line. 

NOTE: END (K) -BEGIN (K) + 1 must be less than or equal to 

PTS and greater than or equal to 5 (for K = 4 & ,5) . 

The result may be zero if a particular pulse is 
to be ignored. This should 

be the same for K-4 and 5 for lamp radiance 
tests so that the integral ratio will be correct. 


CHANNL 


MIN 


MAX 


13 

INTEGER 


3.5 


FLOAT 


6.0 


FLOAT 


The maximum number of channels being referenced 
by any file* It may be less than the actual 
number of channels available, but it must not be 
greater than the number of channels in any of 
the files in the set. Also, it must be 
between 1 & 13 inclusive. 

This is the minimum lamp current that is used. 

It also corresponds to the current in the lamp of 
the first file to be processed in the set. It 
must not be less than 3.0. 

This is the maximum lamp current that is used. It 
also corresponds to the current in the lamp of 
the last file to be processed in the set. It must 
not be greater than 7.0. 
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NOTE: 

If MIN, MAX, & FILMAX default, the following is 

assumed : 


FILE 

6 FILES ACTUALLY PROCESSED 

5 FILES 

ACTUAL! 


LAMP CURRENT 

LAMP 

CURRENT 

1 

3.5 


3.5 

A. 

4.0 


4.125 

3 

4.5 


4.75 

4 

5.0 


5.375 

5 

5,5 


6.0 

6 

6.0 




STEP (1) TOP (1)... TOP (2) 5 

INTEGER 


PTS 


MPXMUX 


SCHAN 


AIRCHN 


TYPE 


» PFILE 


50 

INTEGER 

NO DEFAULT 
(N/D) 
INTEGER 

N/D 

INTEGER 

N/D 

INTEGER 

N/D 

INTEGER 

5 

INTEGER 


The number of data points at the top of pulses 
2 and 3 to be used in the calculations. Fewer 
than this number may be used if the pulses are 
strangely unsymmetric. This number must be less 
than or equal to the pulse width. 

This is the maximum number of reselms in any 
of the 5 pulses. This number must not be less 
than 1, 

This is the number of the multiplexor channel. 

This number must not be less than 1 or greater 
than 13. 

This is the number of the spectral channel that 
MPXMUX represents. This number must not be less 
than 0 or greater than 24, When SCHAN^O, 
the corresponding MPXMUX channel will be ignored. 
This is the number of the aircraft channel 
that MPXMUX represents. No checking is done on 
this number. 

This is a switch to inform the program that pulse 
width calculations are to be made (TYPE*$PWIDTH$) 
This variable specifies which file of the following 
set is to be used for pulse width calculations. 
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» P CONST 3 This is the pulse width constant that is to be used 

FLOATING to convert the lamp Integral into a lamp level 
for this associated multiplexor channel. 

« PANCON 7.9 This is the pulse width constant that is to be 

FLOATING used to convert the panel integral into a panel 
level for the associated multiplexor channel. 

OUTPUT 

PLOTS 

Plot b are made only for those spectral channels whose SCHAN was not equal to 
0. The first plot of each pair is a plot tf the level ratios: 

$ - the level ratio multiplied by 100 for each current level 

* the level ratio multiplied by 10 for each current level 

$8 the level ratio for each current level 

The second plot of each pair is a plot of the integral ratios for each current 
level. The same plotting characters as before are used. After each pair of plots 
follows a table of all the statistics calculated. 

CRRNT lamp current being referenced 

DARK MEAN mean of dark level 

LAMP MEAN mean of selected reselms at top of lamp pulse 

PANEL MEAN mean of selected reselms at top of panel pulse 

LAMP INTEGRL MEAN mean adjusted integral of lamp pulse 

PANEL INTEGRL MEAN mean adjusted integral of panel pulse 

LEVEL RATIO mean of the level ratio 

INTEGRL RATIO mean of the Integral ratio 

DARK DEV 
LAMP DEV 

PANEL DEV Standard deviation of the previous means 

LAMP INTEGRL DEV 

PANEL INTEGRL DEV 

LEVEL RATIO DEV 

INTEGRL RATIO DEV 

The lamp pulse and the panel pulse are both dark-level normalized by sub- 
tracting the mean dark level of the previous line from all the reselms to be 
used in calculations 
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If tYPE=$PWIDTH$, the lamp integral and the panel integral have been divided 
by the pulse width constant. This turns the LAMP INTEGKL & the PANEL INTEGRL 
columns into mean levels of the lamp and panel pulses. 

DATA CARDS SETUP 

1 BEGIN(l) - 10,60,80,60,80, END(1)=50, 70, 90, 70, 90, CHANNL=4 , PTS = 41* 

2 MPXMUX = 1, SCHAN - 2, AIRCHN = 3* 

3 MPXMUX = 2, SCHAN = 3, AIRCHN - 4* 

4 MPXMUX = 3, SCHAN = 4, AIRCHN = 5 * 

5 " MPXMUX = 4, SCHAN = 5, AIRCHN = 6* 

6 INBIN = 701, FILE = 1, UNIT = 4, NS A = 10,100,1,1,100,1 * 

7 FILE = 1, NSA = 10, 100, I, 1, 200, 1 * 

8 FILE = 2, NSA 10, 100, 1, 1, 200, 1 * 

9 FILE=2 , NSA - 10, 100, 1, 1, 200, 1 * 

10 FILE * 3, NSA = 10, 100, 1, 1, 200, 1 * 

11 FILE=3 , NSA - 10,- 100, 1, 1, 200, 1 * 

12 - FILE = 4, NSA = 10, 100, 1, 1, 200, 1 * 

13 FILE = 4, NSA = 10, 100, 1, 1, 200, 1* 

14 ■ FILE = 5, NSA - 10, 100, 1, 1, 200, 1, * 

15 FILE = 5, NSA = 10, 100, 1, 1, 200, 1 * 

16 FILE = 6, NSA - 10, 100, 1, 1, 200, 1 * 

17 FILE = 6, NSA = 10, 100, 1, 1, 200, 1 * 

18 MODE1 = $FINAL$* 

19 MODE1 - $START$ * 

20 FILMAX - A, BEGIN (1) = 10, 60, 80, 60, 80, END(l) = 50, 70, 90, 70, 90, 

CHANNL = 3, MIN = 4.0, MAX = 5.5, PTS = 41* 


21 

MPXMUX=1, 

SCHAN- 3, 

AIRCHN-5 

* 

22 

MPXMUX=3, 

S CHAN-0, 

AIRCHN=7 

* 

23 

MPXMUX=2 , 

SCHAN*=4, 

AIRCHN=6 

* 

24 

FILE=9, NSA-1 , 100 

'• i. 1, 100, 1, 

23 

FILE=9, NSA=1, 100,1 

,1,100,1 

* 

26 

FILE = 8, 

NSA=*1 , 100 

,1,1,100, 

,1 * 
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27 

FILE = 8 

, NSA = 1, 100, 

1 , 1 , 100 , 

1 * 

28 

FILE = 7 

, NSA * 1, 100, 

1 , 1 , 100 , 

1 * 

29 

FILE - 7 

, NSA - 1, 100, 

1 , 1 , 100 , 

1 * 

30 

HO DEI - $ 

FINAL $ * 




31 MODEL = $ STARTS* 

32 FILMAX=1 , BEGIN (1) =1C> 80,120,12,12, END(l)-70, 119,159, 11 ,11. CHANNL=3 , 

3 PTS«1 , TYPE = $PWIDTH$, PF1LE=4* 

33 MPXMUX = 1, SCHAN = 2* 

34 MPXMUX = 2, SCHAN = 3* 

35 MPXMUX = 5, SCHAN = 4* 

36 FILE = 7, NSA = 20,40,1,1,160,1* 

37 MODE! = $STAXT$ * 

38 F1LMAX “4, BEGIN(1)=10,12,12,B0,120, END(1)=70,11,11,119,159, CHANML-3 , PTS-61* 

39 MPXMUX=1, SCHAN” 2 , PCONST=3. 9* 

40 KPXMUX=2, SCHAN=3, FC0NST-17.6* 

4 ! MPXMUX=5 , SCHAN=4 , PCONST=14.0, PANCON=4.9* 


42 

FILE=4 , NSA=20, 120,1 

,1.160,1 * 


43 

FILE=4 , NSA=20 , 

120 , 

1 , 1 , 160, 

1 * 

44 

FILE =5, 

NSA=20, 

120 

, 1 , 1 , 160 

, 1 * 

45 

FILE = 5 

, NS A = 

= 20 , 

120 , 1 , 1 , 

160, 1 * 

46 

FILE = 6 

, NSA = 

20 , 

120 , 1 , 1 , 

160, 1 * 

47 

FILE = 6 

, NSA = 

20 , 

120 , 1 , 1 , 

160, 1 * 

48 

FILE = 7 

, NSA = 

20 , 

120 , 1 , 1 , 

160, 1 * 

49 

FILE - 7 

, NSA = 

20 , 

120 , 1 , 1 , 

160, 1 * 


50 MODE1 » $FINAL$ * 

51 M0DE1 - $RETURN$ * 


This card specifies that the dark area begins at reselm 10 and ends at 
reselm #30. The lamp pulse begins at reselm #60 and ends at reselm it 70. 

The panel ” #80 " 90. 

The lamp 11 (to be integrated) begins at reselm #60 & ends at #70- 

The panel " " ” " 11 80 ” ” ” #90 ' 

Only four channels are to be referenced. A maximum of 41 reselms are in any one 


pulse. The default FACTOR array will be used, A maximum of 6 files will be 
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referenced (default). The minimum current level is 3.5 and the maximum current 
level is 6.0 (default). A maximum of 5 points will be selected at the top of pulse 
2 and 3 (default). 

2-5 These cards specify that multiplexor channels 1, 2, 3, & 4 will correspond to 

spectral channels 2, 3, 4, & 5, and to aircraft channels 3, 4, 5, &6. 

6 This is the input data required by PROCESS under POINT. Tape 701 will be 

mounted on unit 4 and file 1 will be accessed. Every line from 1 through 100 will 
be processed and every point from 1 to 200 will be used. This starts the 
calculation of the means. 

7 This is more input data required by PROCESS under point* INBIN & UNIT do not have 

to be specified since the tape is already mounted. This starts the calculation 
of the standard deviations. 

8-17 These cards correspond to cards 6&7 (see above). It is implied by the 

ordering of the cards that file 1 references the lowest lamp current of 3.5 
amperes and that each successive file references a higher lamp current than the 
previous* It is implied by the number of files referenced that the lamp current 
levels are 3.5, 4.0, 4.5, 5.0, 5.5, & 6.0 amperes. 

IS This tells POINT to perform STEP(6). This section of the program is responsible 

for plotting & tabulating the data that has been calculated previously* 

19 This tells POINT to start all over again with STEP(l). 

20 These cards specify that: 

a maximum of 4 files are to be processed 
The dark area begins at reselm #10 & ends at reselm #50; 
the lamp pulse 11 60 " 70; 

the panel " 80 " 90. 

The lamp pulse (to be integrated) begins at 60 & ends at 70; 
the panel " 80 ” 90. 

Only 3 channels are to be referenced. 

The minimum current level is 4.0 amperes. 

The maximum 11 " 11 5.5 11 

A maximum of 41 reselms are in any 1 pulse. 

The default FACTOR array will be used (default). 
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21-23 


24-29 


30 

31 

32-37 

38-30 


51 


A maximum of 5 points will be selected at the top of pulse 2 & 3 (default). 
These cards specify that multiplexor channels 1 & 2 will correspond to spectral 

channels 3 & 4, and to aircraft channels 3 & 6. Multiplexor channel 3 will be 
ignored. 

These cards perform the same functions as cards 6-17. It is implied by the 

number of files referenced and the order they are referenced that file 9 represents 

a lamp current of 4.0 amperes,file 8 represents a lamp current of 4. 75 amperes, and 

file 7 represents a lamp current of 5.5 amperes, fou will notice that despite 

the original declaration of FILMAX=4, the lamp current levels are based on the 

actual number of files referenced. 

Same as card 18 
Same as card 19 

This requests the program to perform pulse width calculations on the 4th file of 
the next data set. The mean level of the lamp pulse and panel pulse are calculated 
and saved. 

This performs the integral ratio calculations of files 4 through 7. The ratio 
of the lamp integral to the panel integral is printed; the lamp and panel integral 
are printed divided by the pulse width constant for the appropriate channel. 

File 7 (the 4th file of the set) is used for pulse width calculations. 

Returns control to the program that called POINT, 

ERROR COMMENTS 

ERR0R1. 1 ,G. MPXMUX. G. 13. 

This indicates that MPXMUX is either less than 1 or greater than 13. 

ERR0R2 . 0. .G. SCHAN .G. 24. 

This indicates that SCHAN is either less than 0 or greater than 24. 

ERR0R3. FILMAX .L. 1. 

This indicates that the maximum number of files to be processed is less than 


1 . 


ERR0R4 . 1 ,G. CHANNL .G. 13. 

This indicates that the maximum number of channels to be referenced is 
either less than 1 or greater than 13. 

ERR0R5, CURRENT MIN INVALID 

This indicates that the minimum current level is less than 3.0 amperes or 
greater than 7.0 amperes. 
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ERR0R6. CURRENT MAX INVALID. 

This indicates? that the maximum current level is greater than 7.0 amperes 
or less than 3.0 amperes. 

ERR0R7. INSUFFICIENT STORAGE 

Since this program uses dynamic storage allocation for its incore data bases, 
the possibility exists that someone will try to process more data than can physically 
fit into the machine. The main incore data base requires FILMAX + CHANNL* 14+ 

PTS * (2 + CHANNL *5) + 5 locations, (see the section on storage requirements 
for more detail.) 

ERR0R8. TOP(X) ,L.O. 

This indicates that the number of points to be selected at the top of pulse 
X is less than zero. "X" is 1 when referencing the lamp pulse; "X" is 2 when 
referencing the panel pulse. 

ERROR9. TOP (X) ,G. P WIDTH 

This indicates that the number of points to be selected at the top of pulse 
X is larger than the number of reselms in the pulse. "X” is 1 when referencing 
the lamp pulse; M X ,f is 2 when referencing the panel pulse. 

ERRORIO. CHANNL .G. QNCHAN. 

This indicates that the specified maximum number of channels to be referenced 
is greater than the number of channels in the file. 

ERROR11. BEGIN (X) .L. NA. 

This indicates that the beginning reselm number of pulse "X" is less than 
the beginning reselm number for the lines of the file. 

ERR0R12. END(X) .G. NB. 

This indicates that the ending reselm number of pulse "X" is greater than 
the ending reselm number for the lines of the file. 

ERR0R13. PULSE WIDTH .G. PTS. 

This indicates that the width of one of the pulses is larger than the 
declared maximum number of reselms in any one pulse. 

ERR0R14. PULSE WIDTH .L. 5. 

This indicates that the width of one of the pulses is less than 5 reselms. 
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ERR0R15. MORE FILES THAN FILMAX. 

This indicates that the program attempted to process more files than 
originally declared with FILMAX. 

ERR0R16. NO DARK AREA PROVIDED. 

This indicates that END(1)-'REGIN(1)+1 was less than or equal to zero. 
ERR0R17 . DISCRIMINANT .L. 0. 

After the 2nd degree least-squares curve was calculated, the equation must 
be solved for the two zeros. Unfortunately, the discriminate was less than 0, 
indicating that the roots of the equation are imaginary. This might indicate a 
misshapen pulse where the left and right boundaries are not lower than the peak. 
Consequently, the boundaries cannot be determined at the top of the pulse. 
ERR0R18. PTS. *L* 1. 

This indicates that the declared maximum number of reselms in any one pulse 
is less than 1. 

ERR0R19. AVE HEIGHT NOT FOUND 

This indicates that a pulse had a width of 2 ero and TYPE-$PWIDTH$ was not 
specified. 

ERROR20. PFILE. .L. 1. 

This indicates that the relative number of the file relative to the first 
file of the next data set was less than zero. 

SUBROUTINES NEEDED 

The following are a list of external subroutines needed for the correct 
operation of this program: 

ASSIGN. this performs all the dynamic storage allocation necessary in 

the B core box. 

CANCEL. Part of ASSIGN 

ERROR System error processing routine 

FREESP. This routine frees as much of the dynamically allocated space 

as necessary. (Part of ASSIGN ) 

GETNXT. Part of ASSIGN 
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LINK. 

LSQAR. 

PL0T1 . 
PL0T2 . 
PLOT 3 . t 
PLOT 4 . J 

POINT. 


PROCESS. 

SETDIM. 

SQRT. 

ZERO. 

PULSE WIDTH CALCULATIONS 

The purpose for performing pulse width calculations In lieu of just 
calculating mean levels and the level ratio is four-foldj 

1. We found that integration process usually results in a lower standard 
deviation than can be achieved with mean levels. 

2. Calculating integrals goes a lot faster since boundaries don r t have to 
be determined and no curve fitting is ever done. 

3. Changes in the calculated pulse width provide a way to determine if 
changes in the shape of the pulse have occurred and how drastic these 
changes might be. 

4. Through the use of pulse width constants, the results can be messaged by 
the user to account for errors in the digitization process or for errors 

in the data. 

The general method for pulse width calculation is as follows : 

1. The program calculates the mean level of the lamp and panel from a 


This routine provides POINT with a correction to the LAMRAD* s 
point-processing routine PTPROC. (Part of POINT). 

This routine performs a least squares curve fit to a given set 
of data. 

All four of these rcutlTies are used in setting up and printing 
the point plots of the lamp radiance curves. 

This is the program that repetitively calls LAMRAD. with 
different levels to progressively process the data. (See POINT, 
writeup) . 

This is a routine called by POINT, that handles all I/O with 
data tapes (in ERIM format). 

This routine sets up the dimensioning and subscriptions informa- 
tion necessary for the dynamically allocated arrays. 

This routine calculates square roots. 

This routine inserts either a binary or floating point zero in 
all of the requested variables or arrays. 
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file chosen by the user for pulse width calculations. These means (for each 
channel) are saved for later use. 

2. The user then instructs the program to calculate the mean integral of 
the lamp and panel and provides pulse width constants to divide into these 
integrals. Consequently, the integral is no longer an area measurement, but 
a measure of the average height of the pulse. 

3. These adjusted integrals are then used for the integral ratio. Standard 

also 

deviations are/calculated for these adjusted integrals. 

4. This process (# 2 & #3) is repeated for each file of the lamp radiance 
calibration set. 

5. After all the files are processed, the following is printed: 

a. mean lamp height from step $1. 

b. mean panel height from step ill 

c. integral of the lamp. This is the adjusted integral of the lamp multiplied 
by the pulse width constant. 

d. integral of the panel (similar calculations as in c) . 

e. calculated pu]se width for lamp. Integral of lamp (c.) divided by the 
mean lamp height from step #1. 

f. calculated pulse width for panel( similar calculations as in e).v 

g. pulse width constant used for the lamp 

h. pulse width constant used for the panel 

STORAGE REQUIREMENTS 

The LAMRAD program requires 4345 decimal locations before any dynamic storage 
allocation takes place. The formula for determining how much additional storage 
will be allocated for a particular data set is as follows: 

AFTER M0DE1 = $ START $ & 

BEFORE M0DE1 = $FINAL$ 

FILMAX*CHANNL*14+PTS*(CHANNL*5+2) -1-5 
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AFTER M0DE1 • $FINAL$ & 
BEFORE M0DE1 - $START$ 

FILMAX * ( CHANNL * 1 4+2 ) +4 
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PAC <C5 AUG 1565 VERSION) PRCGRAF LIFTING 
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4 
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2 
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1 

2 

2 

4 

c 

t 

7 


1 

2 


C 

1 

2 

2 

4 


t 

1 

e 

9 

C 

1 

2 


EXTERNAL rlM l i l in 
REFERENCES CN 

ncrml NCCE IS integer 

FLOATING FC1NT TCF(2), 
TCT2 (13) , TCT3U3) • 


TCF(2) * FACTCR<24), MN. 

r r nv U C A K CEV» X 

CCEF 


(*AX, TCT(t3*6) f 
Y r F I 2 J f 
Ml, ( 3 J • 


TCT 1 (2 ) t 


CTJU3U TC T 3 t 1 3 ) , SCRT.i ^ E A N , CEV, X, Y, FI2) f 
22, Fit F 2 , F 3 , FCATLM24), MX, CCEF, Ml, A(3>, 
,*X,YY,SAVLAF (13) ,FUILL,PCCNSTI 13*2 ) » S A V F AN l 13),PANCCN, 

ENSIGN ChNNL$(13)t BEGIM5), ENC(5)t FARA^S(7), NSCAL614), 
TITL£(32)t F T ( 5 ) , BCG (3), E P R < 6 5 C ) , F N T S U 3 * 2 ) » 

I T R Y IE 5 C ) , T I T l 3 ) ...... 

ECU 1 VALE NCE < CEN N L S t A I RCh N ) t ( B E G I N , F F X I* L X ) , (cNCtSCEAN) , 
<FCMLN,CATUM t (lTPYU)fERRUH 
unro UM P< prnn = !U5i 1*1 f *8$ 


C1F 


VECTOR 

VECTOR 


VECTOR 

VECTOR 


VECTOR 
VECTOR 
VEC ICR 
VECTOR 
VEC ICR 


DALLES eCC(l) = Hi! , 
VALUES TITLE (1) = IF1, 
162, “ 
tel, 

I C 3 , 

tC4, 

IC6 , 

ice t 
icic 


VALLES 

VALLES 


VALLES 

VALLES 

VALLES 

VALLES 

VALLE 

VALLFS 


1*1 , 18$ 
71-.73I 
$63 t 
$02, 
K2, 
$C5, 
$07, 
$C9, 

$C 12 , 


$E 1 
1 .C-l .41 
1. 5-1.6! 
.67-.S4I 
.58-, 64$ 
.52-.57I 
.48-. 52$ 
.4 l- .46* 


, 3 2- • 3 6 $ 


2.C-2.61 
2.C-2.6S 
1.C-1.4J 
.62-.7CI 
.55-.6C* 

.5C-.5h J 
.46-. 4.9 $_ 

R$G = IFACIANCEI 
FACTCPU ) = 2"~" 

3464.5 ,3076 
1244.3,1017 

C • , c . , c 

Fill * *5, ,8 
FCCNST (1) - 3. v3.v3.f3., 3. ,3. ,3., 3. ,3. ,3., 3 
SAVLAF(l) - 0, , C . , C • t C .,0*tC.fC.tC.fC.tC.tC 

PFLAG = 2 

TIT = ILEVEL RA [ I C $ , I IN TGRL RATIO! 

I .G. NFXFUX . G 


NCE $ 

2725., 133. 8, 1C 16. I, 3484. 5, 1C 16. I, 2166. I, 
8.5,2459.7,21 11.2,17 12.7, 14 65.6, 
7.2f840.65ffc22.65,G.fC.fC.tC.,C., 


► 2 . , 2 « 
► c • , c . 


VCC ICR VALLFS 

ERKU) = $h*l*****ERRCRL. 
VEC1CR VALLES 

LRR(43) - iH 1 1***<*ERRCR2. 
VEC1CR VALLES 

ERR I 85) = fh'l*****ERRCR3. 
VECTOR VALLES 

CRR(127) = $H f M + «**ER«CR4. 
VECTOR VALLES 

ERR ( 16S) = * 1*****EPRCR5 . 

VICTOR VALLES 

F W R < 2 1 1 ) - IH ■ 1* ****ERRCR6. 
CICR VALLFS 


• * 


VtClCR VALLES 

= tH 1 1***** ERRCR7 . 


ERR I2S3 ) 
VECTOR VALLES 


CTCR VALLES 

ERRI255) - IH 1 i * * * ♦ * E R RC.R8 « 
VECTOR V A L L C 5 


.G. NFXFUX .G. 13* 1 * 
C .G. SCFAN . G • 24. •* 
FILRAX -L. I. 1 * 

1 .C. CFANNL .G. 13 
CURRENT MN INVAUC. 1 * 
CURRENT RAX INVALIC. 1 * 
INSUFFICIENT STORAGE.* 
TCF( i ,Il,F l ) .L. C.’* 


E R R C 3 3 7 ) = $h*l*****EFRCR5. TCFI*,ll,F*) .0. F KICTMM 
VECTOR VALLES 

E R R ( 3 7 9 ) = $ H 1 1*****ERRCR1C. CFANNL .G. CNCF AN . 1 * I 

VEC ICR VALUES 

ER R < 4 2 1 I ^ IF 1 ERROR I 1 . EEG 1M ' , I 1 ,F 1 I .L. N A . ■ # i 
VECTOR VALUES 

ERRI462) = ih 1 1<** J MERRCR12. ENU',I1,FM .€. N B . * * % 

VECTCR VALLES 

E R R ( 5 C 5 ) = ih* 1*****ERRCR13. FULSt U I C T F .G. FTS.»* $ 
VECTOR VALLES 

ERRI 54 7 ) » Sh* l*****ERRCR1.4- PULSE UCTF ,L. 5.** $ 
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STtPt 1) 


LP 


{V c P t F I L U S TF AN F I L y A X . 1 * * 
fsC LARK AREA PHCVICEC^ * 
LI SCRIP INANT « L * C. 1 * t 
PT5 ♦ L . 1 • 1 * * 


v e f e [ g f r nct fclnl.'m 


MfCKK HALILS 

= ii -(j*****EFRCRL , > 

VFCTCR VALLES 
F P R < £ 3 l ) - 

V C C T C R nALLES 

hRRU?3> = j|- « i*****6RRCRL7 

V E C R 3 VALLES 

F ,* < < 7 1 ? ) = $j-M*****tPRCPLt 
VfCIC'J V A L L t S 

= th 1 IS- 

V t C K 0 VALLES 

,rj{|7^] r il- * L all ^ >!, *^tRPCt^2C« FFILE *1* 1 • * l 

«iMrPfcMLA 0 FlSTEFn) A 

p,, Cftmig j U yK|?2^),I.UII ,CP££L,CFILe.CUM. uaiTc,l.CMN, 

CN»r' ,m0 .XCCEl ,KCCE2»LMT ,C4l I'VE ,CCL l\ £ , NS i ,<v SL> , 

K£,M,i\t.Kf ,1C (1) .TFL6C , I F * C X , * E $ E S V < 4 I , C F AC T S < .<..9 >, 

i:Tnij(i9)*ciniEii5i»cii$nis)»<.s , '4‘*i : i4fcitC»-i£G* 

CCAivG .CeafvG^HE ,CM*CN$S,UCt-AIV .CI»CCE .C’lECA. LRlCC, 

CFllF.LPEEL.CM>»CS 

eeflCAPlE CAF.M-42?)»nfST,CSTAHT,CSIflWT,NV,NX,KC,l,[F, 

lCF.M:F,l'< : *T.SIAf , T2tS€' s i:T.2,'-<£,»i;i.2,»eAi;F,>,LACElc. 

5MMA<2CJ ,LAIUM24> .ICCC.EI24 1 . I C E A M 2 4 I , Ii'AGEI'jAC) 

M-tM'VLF' NEXT .0. t* f LIVC T I CN tfETLi-lv 
1i<4NSF-r;p ft SIEMIVEXT) 

L !M< . linF'FCC . ) 

L - L + CCHAIV 
CI-/UM = 12 

T l L * X = t 

MV = i.'i 
N k* - t • c 
F ! S = C 

(P-l*i*F*G»5»TCF(F) sl 5*J 

ZEi 3 C«TEEClN-»*PGGINl5)*ENC*-«ENLI5}*FAflAP$*NSCM>. t FlLES» 

«EJL AN^^IM’LiT^FACTLH I L ( , F U f l X . E E G 1 M 1) , E MJ l 1 ) . 0 ANM . f I A * I" A > » U F HI • 

PTStlVPL jPFlLfc , 

h H: \ E V i: R TYPE • C . SFLlLIpf - A N C * FFLAG • E • 2 ♦ PFLAL - - 


PASS - i 

V. HE NEVER FILYA* * L * 1 

TRANSFER U ERR 3 

CR WFENFVcR C h A N M ♦ L * 1 .CP. C F A N N L • • 1 ^ 


1 


1 


T R a A S F L P T C F P P *» 

Ck hrLMrV&R « L * i. • C R * P IN « £ • 7 * 

TRANSFER ft EPR8 

CP UHEMiVfcT M* .G- /. -CR* - i« 

TRANSFER TC FF«c 
,CK v* P b N F V c F7S *L» 1 


TRANSFER TC FFR18 
C R k H E N C V fc R FFILE - L * 1 
TRANSFER TC ERH2C 
NC CF CCNTITICML 
f.F = ENCU) - BEG I N ( 1 ) ♦ 1 
hEMTVFH TCP • L E « C, TRANSFER TC ERRU 
MCLGP LP i FCK * = I » 1 t K -G- b 
iAhhNEVLR ENC(K) - GEGIMK) + 1 
FFENFVcR L NC l K ) - PEGINl*) ♦ T 
r N D t K ) - EEGLMK) + 1 ♦ N E • 


* l ♦ 

• L • 


FTS* TRANSFER TC 
5 . A N C • K • C E • * 


C « TRANSFER TC ERP14 


CCMIKLf 

TFPCLCF LP12* FCR K = i • i t ^ ♦ G « 2 
VH-NFVCP T C F IK) .1* C. 

TRANSFER TC ERR8 

CR WHENEVER TLP(K) # G. t N C { K + 1 ) - P E G l N 1 K 4 1 ) + 1 
# ANC. fcNL 1K+ 1 ) -t EG I N ( K+ I ) + 1 .NE. C 


E H R i 1 
, A i\ C . 
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IRANSFL-R 1C ERK9 
ENC CF CCNCITICNAL 
CCN I IMli 

c i 2£ = F [ L M > *CF A N M * 7 
J - PT$*CHANNl*5 
K K = ASSIGN- (C?C) 

kF F N E V L B KK .L. SI ZF*2 ♦ PTS*2 4 I + 5 f TRANSFER TC EFR7 
KK = ASSIGN. C LAJ^HAC* *FARAN'S*SlZE.**3,SlZE«*-3) 

BLOCK = GtlKM.(e) 

LEVEL = GETN>T.(?LCCK + <t) 

ACCR6S - GE 1 NX 1 • ( G LOCK ) 

KK - ASSIGN* IL&FWAC*fFARA^StFTS*FTSiI*«*3) 

SETClN.lPtAN,FIlMX # ChANNL f 7) 

SEICI^. ICtWFILMK.OANM,?) 

SEir I f*. (FM t F I S tCFANNL ,S ) 

CANCEL. ILAKRAC. ) 

TFRCLGH LFlt FOR K=l,l f K *G - CPANNL 
PANCCN = C* 

PEAT ANC PRIM DATA f^FXr-t> f SCF an , A T HCFN t FCCN5 T t PANCCN 


CFANNL 


■TCT2 m ) .FLAG, 


WHENEVER F > rL > .G. 13 .CR. NFXMX .L. 1, TRANSFER TC E«R1 
WHENEVER SCHAN .G. 2 4 .CR. SCHAN .1. Ct TRANSFER TC ERR2 
khFNEVFR PANCCN ■ L * C • 

PC C N SI (NFXNLX ,2) = 7-9 
CTFERM SE 

PCCNS 1 (KP*M X ,2 ) * PANCCN 
ENC CF CCNrmCNAL 
CFNMS INfXMX) = SCHAN 
PCCNS 1 t KF>M > » 1 ) = FCCNST 
VhFNFVFR SCFAN *E* C, TRANSFER TC L F l 
PRIM FCRNA1 $lF+,T99,H«SPtCTRAL EANC - *.2Ct*J, 

li ILE <SChAN*2-n ,TITLE(SCFAN*2) 

CL N I IME 
PASS = 1 

ZEPcl(TClt..TCT(7fi)»TCIi...TCll(2)*ICT2««*TCT2ll3)*FL^G* 

rC13...TCC3113>.PT...PT(5).LlNES,FNTS...FNTSI2Sn 

WHENEVER pass .E. 1 

WHENEVER CHANNL .G. CNCHANt TRANSFER TC ERP1C 
TERCLCE LF13, ECR K=l,l. K .C. 5 
WHENEVER ENC l K )-EEG I M K > + l .£. C 

WEENEVER K .LE. 3» ( KK= L 1 1 t K K .G . 1 3 »FNT S ( KK ,K ) =» 1) 
TRANSFER TC EF13 
CX WHENEVER BEGIN IK) .1. N A 
TRANSFER TC ERR11 
CR WEENEVER ENC IK) .C. N E 
TRANSFER TC ERR12 
ENC CF CCNCITICNAL 
CCNT IME 

ENC CF CCNCITICNAL 
ELNCTICN K6TCRN 
FLNCTICN RE I CRN 
LINES = LINES + 1 
J = C 

TERCLCE L P 3 « F CR K=L,l, K .G. CEANNL 

TFRCCGE IF22, FCP J=J-U*l f CENNLS(J) .NE. C .CR . J .G. 12 
TCT IK , 1 ) - C . 

(KK = l f ltKK.G.mn«TCT(K,ll = TCTIK.1J 4 
FNT IKK,K,1)/TCF ! 

WHENEVER FLAG .£. C, TRANSFER TC LP2 
TERCLCE LF4, FCH F=2*l, P -G. A 

WEENEVER ENC I F- 1 1 -E EC 1 M F- 11*1 .E. C. TRANSFER TC LPR 
EXCCLTE BCCNCS. 

WEENEVER PASS .E. 1 
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L P ^ 


LP* 


IP3 


STEP! ? ) 


LP2 \ 


INTCRL 


LP 2 


( KKcLEFT , l ,*K .G .R IGF T , I C T t K ,F ) = ICT(K t F> + 

1 FN T IKK ) 

CTFERMSE 

TCI 1 IP-2 ) = C. 
kF‘ EnEVER F . N E * 2 

IKK-LEFT*1,KK*G.P1GH,TCT1 (F-2> - TCTlfF-2) 4 
1 FM ( KK * K, P- 1 I/TOM P-2) I 

CNC CF CCNCI TICNAL 

i X K = L E F 1 * l 5 KK*G*R IGM , T C T <K,f ) = TCTlK.F) 4 
1 IFM(KKtKtF-l)-*cJMFIlESfK,F-in.P.2) 

cm; cr ccnci t icnal 
CCM [Mh 

l>FF NEVER PASS • E • ? »MC t ENClPJ-CEClMiMl . N E * C .AIM. 

1 EM (3 )-eCGIM3) + l .Nt. C 

TCU = TCT 1 <l) / fCT 1 (21 *FACTCMChNMS < J )) 
r C T 2 t K > - T C T 2 ( K ) 4 ITCH - N c A N ( F I L E S r K . * \) « P ♦ 2 •• 

f n r cf t c m: m c n a l 

TFP.CLCF LPtM FCR F=5tl, F .G. t 

WHENEVER EM! ( P - i )-EEGtMF-L 1+1 .£• Ct TRANSFER TC IF* 
GXiECLTF SINFSIW 
UH: NEVER PASS .E. L 

TCT <K , P ) = 1 C T ( K , F ) + C22 
C TFE»*ISE 

TCT1 (F-4 ) = C22 

T C 7 ( K , P ) = TIT(K,F> 4 U22 - NEAMFi L £ S , K * P- il > . F . 2 

cm: cf ccnci t 

C CM IMF 

vPtNCVER PASS .6. 2 . A N C . E NE < 4 ) - 6 EG I M 4 > + l - N E • C .AM. 

1 Fm; (5 I-FEGIMMM . K E • C 

TCTI = TC T1 (1>/TCTI 12 } *FACTCR (CFNMS (J >) 

TCT31K) = T C I 3 ( K ) 4 (TCTI - ► £ A M F I L F S » K p 7 M .P . 2 
EM CF CCMITICNAl 
CCN 1 IMF 

lP=l f L,F.G.5fFHF) = C) 

FLAG = \ 

H.MTICN 

MEMVE* PASS • E . 2, TRANSFER IC.STLCEV 
FILES = FILES 4 I 

V\FEf\r VER FILES .G. FILNAXt TRANSFER TC £ P n 1 S 
J . * C 

TFhCCGF L F 6 t FCP K=l t t, K ,G. CFANM 

Vr P. C L G r L r 2 1 t F C R CFNMSIJ) ^c* C *CO. J .C. 13 

TPRCtCF LF 2 < FCR F=l,l t F . G • 5 

feFCMVER I V F £ « E * 1FKICTM • A N C • FFLAG .£. C 
1 .AM;. F .G, 3, TRANSFER TC L F 2 

hKENEVER F . G . it TRANSFER TC INTGRL 
KFENEVER PNTS (K f P ) . c . C 

MAMFIUI5 f K»F) = C. 

TRANSFER TC L F 2 
bM: CF C CM I TICNAL 

MAMMLES,K,P) - TC T ( K , P4 1) /F N T S < K , F > 

TRANSFER TC LF2 
■whenever LINES .LE. 1 

FEAN IFILEStK.F J = C. 

TRANSFER TC LP2 
- ENC C F CCNCL TICNAL 

MAMFIIES,K,F> = TCT <K,P4 1) /< L IMS - 1) 

CCM IME 

»MNEVhR T 'i P E . E . fPMCTF* . ANC. FFLAG . b . C 
SAVLAR ( K I ” MAN (1 f K ,2 ) 

SAVPAMK) = FEAMltN,3) 

TRANSFER TC LF t 

FR WHENEVER MAN(FILlS,K,2> C. 

RE AMF I LES , K ,6 ) ^ C . 




TRANSFER TC C K 1 

enc cp ccnciticnal 

FEAN(FILES,K,6) = RE AM F I L E S , K »2 ) /F EAM 1 1 L ES • K , 3 )♦ 

1 FAC1CR(CFNNLS(J) ) 

CK1 FhENEVER F EAM F ! LE S , K ,5 J .E. C. 

FEAMF TIES, K, 7) = C. 

TRANSFER 1C LF6 
EM) CF CCNCITICNAL 

FEAMFIL£S,K,7> = F E AN ( F I L E S . K , A ) / F E AN < F I LE S » K , 5 ) * 

1 FACTCB (CFNNLS ( J I I 

LP6 CCMINLE 

PFLAG = 1 
PA'S = 2 
FLNCTICN REHRN 

S7CCEV TFPCLGP L F 7 , FCP K=l,l» K .G. CFANNL 

(P-l»l,P.G.3,Fl=TCT(K,F4l)/FNTS<K,F), 

1 CEVIFILE$»K»P) = SCPT.(Fl) ) 

(P=A,l,P.G.E,Fl=TCI<K,F+l)/(LIN£S-l)t 
l CEV<FILES,K,F) = SCRT.(Fl) ) 

FI = ICI2(K)yMLINES-l) 

CEV(FILES»K,fc) = SCSI. (FI) 

FI = 103 IK) /(LINES-1) 

>'* w ! F ILES.K , 7 ) = SCPT.(Fl) 

L P 7 CCMINLE 

PASS - 1 
FINC1ICN REURA 

STEPU) EXECLIE F LC) 1 . < N SC Al E » 09 f 6 , A , U ) 

ViFENfcVER 1 7 F E . N E . IFKICTM . ANC. ( E N C I 1) - E E G -IK ( 1) + 1 .E. C 

1 .CR. ENC <2 )-eEGIM2)*l .E. C .CP. ENC(3)-6EC-IM3M1 .E. C 

2 .C*. E AC (A > - E EG I N <A) + 1 .E. C 

3 .CR. ENC (5 »-BEGIMS) + l .E. C), TRANSFER TC EPR19 

KK s FfiEESF. (ELCCK, LEVEL, 4CCHES) 

KK = AS S I GN . < C » C ) 

KPENCVtR KK .L. FILFAX*2 + 2, TRANSFER TC E F R 7 
KK = ASSIGN. ILAFRAC. , F 4 R A F S , F I L F A X , F I L F A X ) 

< Is I ♦ 1 , L.G.F UES ,X( I ) = FIN + (FAX - FIN)/(FILE$-I)*I 1-11) 

J = C 

TFPCLGF IFF, FCK K=L,L, K .G. CFANNL 
LP2C IFUCLGF LP2C » FCR J=J+l,l, CFNNLS(J) .NE. C .CP. J .G. 13 

TFKCLGF LP9, FCR F=l,l, P .C. 2 

FFENEVEP F . E . 1 .AND. I E NL ( 2 ) -B EG I N ( 2 ) ♦ 1 . E . C 
1 .CP. ENL l 3 1-6EG IN I 3 ) ♦ 1 .£. C) 

TRANSFER It L F 9 

CP FPeNEVER p .E. 2 .ANC. < EN C ( A )- e DC I M A ) + 1 .£. C 
1 .CP. ENC(5)-BEGIM5M1 .£. C) 

TRANSFER 7 C LF9 
ENC CF CCNCI 1 1 C N A L 

EXL-CLIE FLCT2.(IFAGt,7.,3.C,TCCCC.,lCCC.) 

PRINT FCPFAT HF L , T3 1 ,H* SPECTRAL 8 AN C • , 2C 6 , F • / • , 2C fc ■> I , 

1 TITLE ICFNNLS[J)*2-1>,TTTLEICFNNLS(J)*2)» 

2 T[ I I [F-l ) $2 I, TIT ( IF-l)*2+l) 

THRCLGF LFICt FCP KK=1,L, KK .G. 3 

( 1*1 , I , I .G.F I LES , Y( I) = FE AM I ,K ,F4*j M1C.P . I 3-KK )) 
EXECLIE FLCT3.ieCC(KK) t xm,Vll),FUES) 

IP1C CCMINLE 

EXEClTE FL CT A . ( 8 , F SG ) 

PRIM FCPFAT » I P 9/ 1H9 , TAC ,F- ' CLRRENT 1 * 1 
LP9 CCMINLE 

PRINT FCRFAT i l F 1 , T 3 3 , F • L A F F PANEL*, T93, 

1 F'LAFF PANEL LEVEL INTEGFL'/lF ,T49, 

2 F 'LEVEL INTEGSL* ,TIC8,2(F'PATIC *),T9, 

3 2 ( H 1 C A R K LAFP FANEL I N T E G R L 1 NT E GR L * , T 69 ) , T 13C / 

A F* CRPNT 1 , 1 IF 1 FE AN • ,SA ) , 2 IF * RAT IC * , S 3 ) ,T7C * 

5 ?lh'CtV*,SE)//*J 
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FULL = f'CAMFFILfc t K,«)*PCCNST<J,i>/SAVLAMK) 

F U C R = FFAMFFlLF.K.SUFCCNSnj^l/SAVFAMM 
TPhCLO LFllt FCH KK=i,l, *K ,C, FILES 

PRIM FCHVM 1 1 h v F5.3t*tF7.2'Sl>»T47«2(F<3.2 f <;i),T67t7(F7.2,SLJ4i, 

1 MKK) f ^AMKK t K t l)...MiflMHK,K f 7}, 

2 LEV IKK ,K ,1) •••DEV (KK ,K ,7 ) 

LPll CENTIME 

PRIM FCRf*AT tF'-N E C I: L T A L '//*i 
TFRCLCP L P ? 3 t F'C« K K ■= I » I » KK ,G . FILES 

PRIM FCPMT !2F *F1C.4*$, FAC7CR<CFMVLS(J)]*CEV(KK,K t 4 )*2 

1 / F E A |\ T K K * K , 4 ) 

L P 2 3 CCKUMl 

FiUNr FCRf^AT tP '-FiCTCP 1 ' , T2,F* ) = 1 * F 1 C . 5 * S 2 t F * L W / C N * C P - S T IE H - 1 P 1 / 

1 F It - XiCCV 

1 F i * - X1CVH* 8 - Xl**i, CHMSIJJt 

i FflCTCP ICFNMS ( J) ) 

PRIM FCRHM tF'-RbSLLIS F RL R FLLSE fc U T F CALCL L AT 1 CN S , • / 

1 F * C INTEGRAL CF L * I* F • * S 1C , F 1 C . 4 / 

1 F« I ^ T L- 0 H -fi L CF F A N E L • , SS ,F ]C . A / 

2 F' FEAN LA*F HtlLPT f f$lC»FLC.4/ 

2 F ■ MAN PANEL FEIGFT',S 9^10.4*1, 

3 PF ICL*SAUAF { K ) »FlrtICF*SAvPAN (K) fSAvLAMK )#SAVPAMK ) 

PRINT FCPMT $ F 1 CALCLLATtC FUSE « 1 C T F 1 ♦ S 4 f F I C .4 * F 1 - LAPP*/ 

3 IF , T 3C ? F 1C. 4 ,F ' - PANEL 1 / 

4 F * FUSE MLTF CCNS1AM IS EC • , S 1 » F 1C ,4 . 1- f - LAPP 1 / 

E IF , T3C t F 1C. A ,P 1 - PANEL 1 M* 

6 PF TCLiFMCFiFCCNST IJtl ) tFCCNST < J f 2 ) 

LPE CCMIME 

TVFE = 1 $ 

PFLAG = 2 



1 K= ] , ] ,K. 

G. 1 

?tSAVlAr* (K ) 


FLNC1 1 C ^ 

RE 1 1 R N 

eitF 

PIMM FCHM1 

FPR U ) 

£ R R 1 

F P R l R . 

I * 1 
TRANSFER 

TC 

ELEF 

ERR2 

1 = 4 3 
TRANSFER 

TC 

ELEF 

ERR 3 

I ^ Ht 
TRANSFER 

!C 

ELEF 

ERH4 

I " 12 7 

TRAN SHE K 

TC 

FLEF 

ERRS 

I = 169 

transfer 

TC 

HLEF 

ERR E 

I = 211 
TRANSFER 

TC 

CLEF 

ERR 7 

1 - 2E3 
TRANSFER 

TC 

CLEF 

ERR8 

I = 2 9* 
TRANSFER 

1C 

1 2 L F F I 

ERRS 

1 = 3 3 7 
TRANSFER 

TC 

C LEFT 

ERR 1C 

I = 37S 
TRANSFER 

TC 

CLEF 

ERR 11 

1 = 4 2 L 
TRANSFER 

TC 

ELEF 1 

ERR12 

1=463 

TRANSFER 

TC 

HLEF1 

ERR13 

I ECS 

transfer 

TC 

ELEF 

ERR U 

I - 5*7 
TRANSFER 

TC 

ELEF 

ERR 1 5 

i = ses 

TRANSFER 

TC 

PLEF 
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EBRie i = #. 2 i 



transfer 

TC 

ELEF 

ERR n 

I = £73 
transfer 

TC 

CLEF 

ERR L fi 

1 = 715 
TRANSFER 

TC 

ELEF 

ER R 1 S 

I : 757 

transfer 

TC 

ELEF 

ERR2C 

I = 7 5 S 
TRANSFER 

TC 

ELEF 


ei_EH PRIM FCRPA1 EFPUl.K 

ERPCR . 


LP1S 

LP18 


INTERNAL Fl-NCl ICI^ F 1 FRCC. 

TFRCLGP LP1E, FCR P = l,l. F .G. 5 B , r ,. (n 

PI ( F > = FTIF) ♦ 1 

TFPCLCF LFIM FCP K*l»lt K *C. 13 

► FEMrVtP CHNMSIK) .E. C. IPAMfEF TC LFIS 


1 = 1 + 1 
V,hEKE\/ER F •£• 1 
f M ( FI IP ) • I *P> 
C TFEPU I SE 

FM (f TIM 'I .F) 
ERG CF CCNCIT1CML 
Cf MIME 


= CATLF IK ) 

= CATLF IK) - T C T I 1 » 1 ) 


EKC CF CCNCI T1CNAL 
CCM IME 
FLKCT1CN REURN 
ENC CF FLNC1ICN 


LP14 

L P 1 8 
CLll 

LPU 
Ct 1 2 


IMtRKAL FLKCTICK 6CLNC$. 
F.FENfVER P .E. 2 
LEM = 1 
RIGM = Fill) 

TRANSFER 1C PIN 
EKE CF CCNC t T 1 CM L 


> 1 = C - 

uwr i Gl- LFl^t FC R KK-l?lf KK • FTIF~1) 

*FENE^CR FM(KK f K»P-U *G. ► > 1 • Ml = FMIKK,K»F-U 
CM IMF 

t-HCLGH L F 1 5 i FCP K K = 1 • 1 • H K -C. FT(F-l) 

VFENEVER F>l*FIF-l) *L* FM ( K K »K *F- 1 ) ♦ TRANSFER TC CLT1 

CM IME 
K = I 

EFT = KK t 

FRCCCF LPl^f FCP K K 3 f T (F- 1 ) »- 1 * K K *1* 1 
FhE NEVER l u >l*F tP — 1 > .1. FM(KK»K«F“1)» T.PANSFEF TC CLTc 

.CN I INLE 
:K = PTIF-1) 

! 1 C F 1 - |< K 

• FEKtVER PM (LEFT ,K .F-l J .L. F M ( P 1 G FT > K . P- 1 >- 1C .AKC. P .E. 
KK = LEM 

u = i t i,i.G.PTip-i)M>m= FMU,K,p-n,>xm= n 

EKECLIE L£CAM.I>! < i>Vll)*ITPV»ltFT(F-l)»A(l)*AI2)»A(2)) 


A 
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HTN 


LP1? 


f:X£CLU' MJIS1. 

C R WHtMVFR FM(RTGFT»K f P-L) • L * PM (lEfT f K»F-l )-lC .*M. F .E. 4 
K K = R IGF T 

( T=lti*I -C.PT(F-l) ,UIIJs FMI i*K,F- 1) t >X( I )- I) 

EXeCLTE LSCAR.(>MYYm,ITRY f l »mF-l),MlifM2)*M2)) 
bXGClTF H DJIM. 

cm cf ccmmticml 

F I C C L E = (LEFT 4 RIGHT )/2 

IFF! = MCCU - TrF(F-2)/2. + , E 

RIGHT = MCCLG + TCP(F-2>/2* - .5 

WFEMiVER LFFT » L * 1, LEFT - 1 

WHENEVER RIGHT .G. FTIP-1M RIGHT = FTlF-1) 

PMMKtP-1) = FMMK f F-n + RIGHT - LEFT ♦ 1 
FIMTMN RhRRN 
EM CF FINCTICN 


IMEPH.At FUCTICF Si'HFSM 

WHENEVER FTlF-ll/Z^ * E . PT(F-l), FTIF-I) = FT(F-l) - 1 
Q 22 = PM(ltKfF-l ) 

CCcF ^ 2. 

THRCIGH LP17, FCR KK»2 f l, K K . C ♦ FI(F-l) - 1 

WHFMVER GCFF .E. 4. 

CCCF * 2. 

L IHERWI £H 
CCCF = 4. 

EM CF CCNCITICNAL 

Q77 = i;?2 4 CttF*FM (KKfKtP-1 ) 

CCN 1 IMt 

C 22 = C22 4 FMIFTtf-1 ) f K f P-1 ) 

WHEMVER TYPE . E * FFWICTH* 

C22 - C22/PCCNST(JtF-4) 

CM CF ccm iiicn/il 

C .22 - G 2 2 / 2 • 

HLMTICN R E T l P N 
FM CF FIMTICK 


INTERNAL PLMTICN ACJIST. 

Mil ^ Mil - PM IKK »K »P-1 ) 

FI = M2J.F.2 - 4.*M3)*M1) 

WHENEVER FI . L * C.fTRANSFER TC ERP17 
F2 = {-.0(2) 4 KHMFUl/(2.^[3}) 

F2 ^ (-M2) - SCPT. (FI ) )/(2.*M3 )) 

WHENEVER F 2 .L* F2 .AM. KK . E . RIGHT 
LEFT = F 2 

CR WHEMvkH F? *GC. F3 .ANC. KK . E * P I C FT 
LEFT = F3 

CR WHENEVER F2 . L • F3 
RIGHT = F 2 
CTFFRW IFF 
kICFT = F? 

EM CF CCMI 1 1 C h /N L 
FLMTICN RE11PN 
EM C F F LNCT ICN 
FM CF FLM11CN 
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Appendix E 
PROGRAM THERML. 

The purpose of this program is to calculate a set of statistics from a given set of files. 

This program is very similar to CALIB. with the exceptions that a tape data base is not used 
and all data is assumed to be of one spectral channel (infrared). The plots and table produced 
are identical in format to CALIB. 

The user must specify the beginning and ending reselm numbers of the cold plate, the hot 
plate, and the ambient area. Both the hot plate and the ambient area reselms are corrected 
by subtracting the mean of the cold plate from the previous line. We are assuming that the cold 
plate mean does not change significantly from line to line. 

After this, all the reselms in the cold plate, the hot plate, and the ambient area are aver- 
aged to get a mean cold plate level, a mean hot plate level, and a mean ambient level. Then 
the standard deviation of each reselm from its respective mean is calculated. 

The aforementioned calculations are repeated for each of the files requested. After all 
the files have been processed, they are sorted into ascending numerical order of date. Plots 
are then made of the mean cold plate level, mean hot plate level, and the mean ambient le^el 
versus mission number. This is followed by a plot of the respective standard deviation versus 
mission number and a tabulation of all the calculated statistics, dates, mission numbers and 
comments. 
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STEP (1) - PERFORMED THE 1ST TIME THERML. -IS CALLED OR WHENEVER PROCESS. READS 
M0DE1 = $ START $ 

This section of the program performs the initialization necessary for POINT 
and sets up the default values for the input variables. Data cards are then read 
that may set the maximum number of files to be read, the maximum number of 
points in any one pulse, and the beginning and ending reselm numbers of the 
cold plate, the hot plate, and the ambient area. The input variables are then 
checked for reasonableness. After this, the dynamic storage for the program 
is allocated, and the program returns control to POINT. 

STEP (2) - PERFORMED BEFORE EACH FILE IS PROCESSED. 

First the totalling array, the point counter array, the date array, the line 
counter, and the mission number are zeroed. 

If this is the first pass through the data, the mission number, the date, 
and any special comment are read off data cards. These variables are checked for 
reasonableness as well as the countinued reasonableness of the beginning and ending 
reselm numbers for the pulses. The program returns control to POINT. If this 
is pass two, the program does nothing except return to POINT. 

STEP (3) - PERFORMED BEFORE EACH POINT OF A SCAN LINE IS PROCESSED. 

No action taken. 

STEP (4) - PERFORMED AFTER EACH LINE IN THE FILE IS READ. 

First the line counter is incremented, and the mean of the cold plate is 
determined (This value will be used In the cold plate correction of the next 
line). If this Is the first line bf the file to be processed, nothing else is 
done and' the program returns control to POINT. 

If this is pass one, the level of each reselm within each of the three 
possible pulses is accumulated in the totalling array. If this is the second pass 
through the data, the standard deviation of the reselm levels from the mean 
level within each of the 3 possible pulses is accumulated in the totalling 
arrays. After both pass 1 or pass 2, the point counter array is zeroed, and 
control returns to PQTNT, 
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STEP (5) PERFORMED AFTER AN ENTIRE FILE HAS BEEN PROCESSED 

If this is the first pass through the file, the file counter is incremented 
and then checked to see if it has exceeded the declared maximum number of files 
to be processed. Then, the means are calculated for the cold plate, the hot plate, 
and the ambient area; the mission number and data are saved; the special comment 
is stored; the pass number is set to 2; and control returns to POINT* If this 
is the 2nd pass through the data, the standard deviation is determined for the 
cold plate, the hot plate, and the ambient area; the pass number is set to 1; 
and control returns to POINT. 

STEP (6) PERFORMED WHENEVER PROCESS READS M0DE1 - $FINAL$ 

After the statistics are put into ascending order of date, the plotting- image 
array is set up- The first plot is then made up of the means of the cold plate, 
the hot plate, and the ambient area versus the mission numbers. The 2nd plot is 
made of the standard deviations of the cold plate, the hot plate, and the ambient 
area versus the mission number* In both plots, the mission # of the last file 
processed is largest mission number plotted* The smallest mission number plotted 
is 1 or the mission number of the last file processed minus 61, whichever is larger. 
After this, all of the statistics, dates, mission numbers, and special comments 
are tabulated, and control returns to POINT. 

PTPROC* 

This routine saves the value of each reselm that is within the specified 
boundaries of each pulse. These values are stored in an array according to the 
position relative to the beginning of the pulse (the point counter for each pulse) 
and according to the pulse number (3 possible). Pulse number 1 is considered 
the cold plate, and its uncorrected value is saved. All the remaining pulses 
are cold- plate corrected before being stored. That is, the average cold plate 
level of the previous line is subtracted from each point value before it is stored* 
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INPUT VARIABLES 

WHEN READ VARIABLE DEFAULT/MODE 

STEP(l) FILMAX 1 

INTEGER 

STEP (1) BEGIN(l) . . . 0 

BEGIN (3) INTEGER 

STEP(l) END (1 ) . . . 0 

END(3) INTEGER 

STEP (1) PTS 50 

INTEGER 

STEP (2) MISSON 0 

INTEGER 

STEP(2) FDATE(l) . * . FDATE(3) 0 

INTEGER 


DESCRIPTION 

The maximum number of files that are to be 
processed. Fewer files may be processed, 
hilt not more. This number can not be 
less than 1. 

This is an array of beginning reselm 

numbers for each of 3 possible pulses. 
(Begin (1) - cold plate , Begin (2) - hot 
plate, Begin (3) - ambient area.) 

These numbers must not be less than the 

beginning reselm numbers for the line. 

This is an array of ending reselm numbers 
for each 3 possible pulses. These numbers 
must not be greater than the ending reselm 
number for the line. 

This is the maximum number of reselms in 
any of the three possible pulses. This 
number must not be less than 1. 

This is the mission number that is to be 
associated with this file. This number 
must not be less than 1, 

This is the date (in month, day, year, form) 
that is to be associated with this file. 
FDATE(l) must not be less than 0 or greater 
than 12. 

FDATE(2) must not be less than 0 or greater 
than 31. 


STEP (2) COMMNT(l) . . . 

C0MMNT(8) 


BLANKS 

ALPHANUMERIC 


FDATE (3) must not be less than 0 or greater 
than 99. 

This is any special comment that is to 
be associated with this file. No error 
checking is performed on this array. 
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OUTPUT 

PLOTS 

Only two plots are made by this program* The first plot is of the means 
of the cold plate, the hot plate, and the ambient area versus the mission number: 
$ = mean of the cold plate 
* - mean of the hot plate 
8 = mean of the ambient area 

The second plot is the standard deviations of the cold plate, the hot plate, 
and the ambient area versus the mission number. The same plotting characters 


as before are used. 

After the plots follows a table of all the statistics calculated along with 


the dates, the mission numbers, and the special comments: 


DATE 

MISSION 

COLD PLATE MEAN 
HOT PLATE MEAN 
AMR TENT MEAN 
COLD PLATE DEV 
HOT PLATE DEV 
AMBIENT DEV 


the date associated with this data 

the mission # associated with this data 

the mean level of the cold plate 

the mean level of the hot plate 

the mean level of the ambient area 

the standard deviation of the cold plate 

the standard deviation of the hot plate 

the standard deviation of the ambient area 


DATA CARDS SETUP 

1. FILMAX = 4, BEGIN (1) - 10,60,80, END (1) = 50, 70, 90, 

PTS = 41* 

2. XNBIN = 777, UNIT - 4, FILE = 8, NSA = 1, 100, 1, 1, 200, 

3. MISS0N = 60, FDATE(l) - 10,29,50, 

COMMNT = $MY BIRTHDAY $* 


4. 

FILE =8, NSA = 1, 100, 

1, 1, 200, 

1 * 

5. 

FILE = 9, NSA = 1, 100, 

1, 1, 200, 

1 * 

6. 

MISSON = 62, FDATE(l) 

- 11, 01, 72, 


COMMNT = $WHEN I WAS HIRED $* 


7. 

EILE = 9, NSA = 1, 100, 

1, Is 200, 

1* 

8. 

FILE = 7, NSA = 1, 100, 

1, 1, 200, 

1* 


1 * 
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9. MTS SON = 61, FDATE(l) = 02, 31, 70, 

COMMNT = $ ABSURDITY** 2 $ * 

10. FILE = 7, NS A = 1, 100, 1, 1, 200, 1 * 

11. MODEl = $FINAL$ * 

12. MODEl =* $ RETURN? * 

1. This card specifies that the maximum number of files to be processed is four, 
the beginning reselm # of the cold plate Is 10, and the ending reselm // is 50. 

The beginning reselm # of the hot plate is 60 and the ending reselm // Is 70. The 
beginning reselm // of the ambient area is 80 and the ending reselm // is 90. A 
maximum of 41 reselms are in any one of the pulses. 

2. This is the input data required by PROCESS under POINT. Tape 777 will be 
mounted on unit 4 and file 8 will be accessed. Every line from 1 to 100 will be 
processed and every reselm from 1 to 2 00 will be used. This starts the calculation 
of the means by pass 1. 

3. These cards specify that mission number 61, the date 10, 29, 50, and the 
comment $MY BIRTHDAY? will be associated with this file. 

4. This is more input data required by PROCESS under POINT. INBIN & UNIT do not 
have to be specified since the tape is already mounted. This starts the calculation 
of the standard deviations for pass 2. 

5. See card 2 

6. See card 3 

7. See card 4 

8. See card 2 

9. See card 3 

10. See card 4 

11. This tells POINT to perform STEP (6). This section of the program is 
responsible for plotting and tabulating the previously calculated statistics. 

The last mission number that will be plotted is mission number 61 since this is 
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the mission number of the last file processed. Mission number 1 will be the 
first plotted. Mission number 62 will not get plotted but will be included in the 
tabulated statistics. As can be seen, fewer files may be processed than originally 
declared, and the files may be in any order (being careful of course to have the 
largest desired mission number to be plotted as the last file to be processed). 

12. This tells POINT to return control to the original program that called POINT. 

ERROR COMMENTS 
ERROR 1. FILM AX. L. 1. 

This indicates that the maximum number of files to be processed is less than 

1 . 

ERROR 2. PTS.L.I. 

This indicates that the declared maximum number of reselms of any one pulse 
is less than 1. 

ERROR 3. PTS.L. PULSE WDTH.L.l. 

This indicates that the width of one of the pulses is either less than 1 or 
greater than the declared maximum number of reselms (PTS). 

ERROR 4. INSUFFICIENT STORAGE. 

Since the program uses dynamic storage allocation for its incore data bases* 
the possibility exists that someone will try to process more data than can physically 
fit into the machine. The main in-core data base requires* 

FILMAX*18 + PTS*3 + 9 + FILMAX locations when FILMAX is greater than 
or equal to 6. It requires FILMAX* 18+PTS* 3+1 5 locations when FILMAX is less than 6. 
(See the section on storage requirements,). 

ERRORS . MISSON .L.O. 

This indicates that the mission number to.be associated with the file is less 
than 0. 

ERR0R6. INVALID MONTH. 

This indicates that FDATE(l) is either less than 0 or greater 12. 

ERR0R7 . INVALID DAY. 

This indicates that FDATE(2) is either less than 0 or greater than 31. 

ERROR 8 . INVALID YEAR. 

This indicates that FDATE(3) is either less than 0 or greater than 99. 
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ERR0R9 . BEGIN (X) . L. NA* 

This indicates that the beginning reselm number of pulse "X" is less than 
the beginning reselm number for the lines of the file, 

ERRORIO. END(X) . G.NE. 

This indicates that the ending reselm number of pulse "X" is greater than 
the ending reselm number for the lines of the file. 

ERR0R11. MORE FILES THAN FILMAX. 

This indicates that the program attempted to process more files than 
originally declared with FILMAX. 


SUBROUTINES NEEDED 

The following is a list of external subroutines needed for the correct 
operation of this program: 

ASSIGN * this routine performs all the dynamic allocations 

necessary in the A core box. 

System error processing routine 

this routine provides POINT with a connection to THERML.'s 
point-processing routine PTPROC. (Part of POINT.) 

All four of these routines are used in setting up and 
printing the point plots of the thermal statistics 
calculated. 

this is the program that repetitively calls THERML. with 
different levels to progressively process the data (see 
POINT, writeup). 

This routine, called by POINT., handles all 1/0 
with data tapes (in ERIM format). 

This routine calculates square roots. 

this routine inserts either a binary or floating point 
zero in all of the requested variables or arrays. 

In addition to these external functions, a define package must also be included. 
This package Inserts 3 new operators into the MAD compiler and defines what 
machine instructions they each represent. These three operators are .LI., .SI., 


ERROR. 

LINK. 


PL0T1 . I 
PLOT 2 . ) 
PL0T3. ( 
PL0T4 . ' 


POINT. 


PROCESS. 


SQRT. 

ZERO. 
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and .LIF. and are called the indirection operators. They are used in this program 
in conjunction with the dynamic storage allocation routines for accessing 
dynamically allocated aroays without using the usual subscription routines provided 
by the MAD compiler. Instead, they use the two internal functions S32< and S82. 
defined in the program for all subscripting. 

S32 . is used for accessing a 2 -dimensional array with maximum dimensions 
of (FILMAX, 3). 

S82. is used for accessing a 2 '“dimensional array with maximum dimensions 
of (FILMAX, 8), 

STORAGE REQUIREMENTS 

This program requires 1959 decimal locations of core before any dynamic storage 
allocation takes place. The equation for determining how much storage will be 
allocated with any given set of input cards is as follows: 

FILMAX* 19 + PTS *3 + 9 when FILMAX >_ 6 
FILMAX* 18 + PTS *3+15 when FILMAX < 6 


165 



2?rjm 


FORMERLY WILLOW RUN LABORATORIES. THE UNIVERSITY OF MICHIGAN 


MAC (09 AUG VERSION) PROGRAM LISTING 


CXTEKNAL FUNCTION THERML . < X» Y ) 

s> I: F fc KE NC F S_ t N 

* ULFINt PACKAGE FCR LOAD INDIRECT ( .L I . f.L IF- ) AND 
STORE INDIRECT!. SI.). THE V ARI ABLES BEIN G R EFERENCED 
MIST HE IN THE A CCkE OCX. IF THE AGGRESS BEING REFERENCED 
IS IN T H i" B C OR E B C X THE VARIABLE C CNT AI IS I N G_ J HE .ADDRESS . 
MIST HAVE A NODE CF 5. 

OF FINE L NARY OPERATOR tLl.i PRECEDENCE SAME AS .A6S„ 

M-CCE STRUCTURE 1 - .LI. I 

JMP *+5|AC,*+l --- 

JMP * + 6 i M Q f ♦ + 1 

JMP *+7,LA,*+l ..... — 

CL A* R 

C L I AC 

STO T 

*' T < T 

JMP *-b _ — - 

SLVv T 

.JMP __ *-7 _ 

END 

MODE STRUCTURE 0 - .LI._ 0, SAME SECUENCE AS .LI. 1 
DEFINE UNARY OPERATOR .L1F.* PRECLCt.NCE SAME AS .ABS. 

Mere STRICTURE 1 = .LIF. 0, SAME SECUENCE AS .LI. L 
KUCt STRUCTURE C = • L 1 F • 1, SAME SEQUENCE AS .LI. 1 

DEFINE BINARY OPE RATOR .SI., PRECFDENCb SAME AS = 

MCI:L STRUCTURE 1 = 1 .SI. 1 
JMP *+l«Dr»*+4 

JMP **8.AC,*+1 

JMP **9|M-G,**1 

JMP *4lC,LA,*43 

JMP _ H11,AC»*+1 

JMP *4l2,LA,*+t 

LL A ft _ 

S T {* A 

l IT 7 __ _____ 

STD* " A 
JMP *-2 

S ID* " A ' 

JMP 

SLIr* A 

JMP *- 6 

STL T 

JMP *-lC 

SLv. T 


JMP *-12 

CNC 


MODE 

STRUCTURE 

C = 

c .si. 

0, 

SAf'E 

SECUENCE 

AS 

1 .SI. 

1 

MODE 

STRUCTURE 

1 = 

1 .SI. 

0, 

SAME 

SECUENCE 

AS 

1 .SI. 

1 

MODE 

STRUCTURE 

C = 

c .si. 

1* 

SAME 

SECUENCE 

AS 

L .SI. 

L 


NORMAL > UOC IT INTEGER 


FLOATING POINT TCTU),$DRT., X , Y t F l , F2 , FDAT UM 

DIMENSION BEG I M 3 ) , E NC"(’3 NP A"R AMS < 6 ) * NSC A LEU ) » P f ( 3 ) , BCC ( 3). 
. I CRRU6 3) tFCATE 13) .CUM MM (6) 

ECU I VALENC E ( FUATIK. DATUM ( 1) ) I M EAN » PAR AMS ( 3)) V 

_J LG F V yPARAM SU ) ) , ( FNT , P A P AMS ( 5 )) T 

2 (DATE , PAR AMS (6 ) ) f ( CM M NT . P A U AMS ( 7 I ) « 

3 ( TAG,PARAMS I 8 ) ) 

VECTOR VALLES CCO(l) = $ $ $ i ' 7 Y* & ~ V t 8 * 
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STtPUl 


L P 1 


STEP! 2 ) 


VGCfC'R VALLES I* SGI - if-'EAN LEVEL* 

VECIHK VALLES PSG2 = ^DEVIATION LEVELS 

vector VALLES err ML*!* — !L ** •- J 

V EC T ( • R VALLES H R R l b 4 3 I = ili # it * * * * * E RKCR2 • PTS • L * !•'* 

VECTOR VALLES E K R ( 0 B 5 ) = *h M ** ** * ER* CR3 _ PTS ■ l* • HUL * E ^ * L . . 

VECTOR VALLES E R K < 12 7) - fth ■ 1 *****EKKCR4 . INSUFFICIENT STORAGE* 
VECTOR VALLES E R R ( 1 1 9 ) = tH M ***** t*RLR5 • __P i SS I CN *L. O* 

V L C TOR VALLES ERR 1211) = SH • 1 * *** * EH RCR6 * INVALID PLINTH.** 

VECTOR VALUES ERR l 2 b 3 ) = m M ***** ERRCR7 . 1N V A L IC_C AY .■_*__ 

VtCTCfT VALLES tRR(295) = $HM ***** ERKCRH . INVALID YEAR.** 

VECTOR VALLES ERM337) = $H 1 1 * * * * * E R K C K9 . _B E G I N ( • ♦ ll.H 1 ) . I . 
VEC1UR VALLES ERRT379) = $h * 1 * * *** t RRCR 1 0 . ENLKMIlfHM *0* NU. 
VEC TOR VALLES L R R ( A 2 L > = *h • 1 *** **ERRCR 1 1 . POKE FILE5 THAN F1LMAX 
STATEMENT LABEL STCPtn 

ERASABLE J L is K t ? 2 5 > , C L M T , CKE E L t CF I L £ , CL l *E , TWRj 1 E y 0 C HA h , 

q ^ o S t PR’cG T p COE 1 t P OL L 2 f UN l T T C A L l N fc t C L L I N E t NS A » N S B r 

KS,NA f NB tKP 1 10(1) t IFLAGi I FACKi RESERVE ) ,CFACTR (49 ) » _ 

CTI TL2 < 19 ) tOT ITLF. ( 19 ) ,CLI ST t 19 ) . CSPARE ( A6 ) , OFL AGt 

QDANG »CL ANC »CRP » GNA » CNSS ? CNCHANt CP CD t » CR EC A » OR ECC * 

UFILE .CREEL tCNhRCS ' . ^ 

ERASABLE D A I A ( 4? 3 ) » 1 T E ST » CS TART t OSTAkT r EW tLj_I_P»_ 


h 

i •** 

* 

i 

i 

$_ 

* " 

* t 
* 


» i ' T .F_ ~ ~ ' I--.-.--— — — — 

[ CP t N CP ? NEXT7 s’T’ART2 J REAC1 2,REACL2»RLACP2 t LAREL2, 

E XTRA (20 iDATUP (24) , 1CCDE (24 ) t I CHAM 24 ) t I M AGE (5 70 
INTERNAL FUNCTION S32.(S1»S2) = 3*(Sl-l) + S2 

INTERNAL FUNCTION SC2.(S3.S4) = 0* ( S3 - 1 ) S4 

M-CMiVER NEXT .G . 6, FUNCTION RETURN 

TRANSFER TC STEP ( NEXT) 

LINK. (PTPlVcC. )" 

L = L + CCHAN - - - - 

FTLPAX = 1 
HT C = 5C 

ZERU.(l>tr,IN. ..KEG IN t 3) , EN UVi . £Wl 3 > , P AR AH S . . .P A R AMS (81 , F I L E S , 

NSCALE . . .NSCAIE U) ) 

UtAr, A'lC PRINT DA t A FI L V A X « BEGIN! 11* END I 1 ) * PTS 

WHENEVER FiLMAX .L. 1 . 

TRANSFER TC EPR1 

Ci< In H E N E V t K PTS . L . l_ _ 

TRANSFER TC CRP2 
Ehi: UF CCnDITICNAL 


IIPLUCiH LP1» FCR 1=1*1-* I *tj« 3 

VHtNtVCK END ( l ) -BEG IN 1 1 ) + 1_ . L. I .CR. EAC 1 1 )-BEGIN( I) +1 -G* PTS, 
TRANSFER TC C HR 3 

CCMIMJh 

J = F 1LHAX 

in H t N E V L- i< F ILHA X .L. 6* J = 6 

I = F l L P A X #3 
K = PTS* 3 
KK = ASSIGN. ( C *0 > 

VnEENEvER KK .L. 1*2 + HLPAX» X2 * 

PASS - 1 

KK = ASSIGN. ( THERHL. »PARAHS*FILHAX*J) 

KK = ASS I GN'.l C7P ARA?S * I , I , K , F U P AX , F l LM AX *8 , F I LP AX » 

FI = HNCI1I - BEGIN! n ♦ 1 

FLNCriCN RETLRN 

ZCP.O. I TCI... TTT!^* »PT. ..PTI 3) .FDAT E...FCATEI3) . FL AGiLUit ES.*- 
H USCN ) 

b I: N C V E K PAS S .E . I 


K + 9 * J, TRANSFER TC ERK4 


IK =1.1 *K.G.8»CCPPNT (K) = $ $1 

«CAC AM) PRINT DATA H I S SCN ,...£0 A.IE..! l.U . CQtIPiNT.! 1 1 1 

WHENEVER MSSCN .1. 0 

transfer TO. ER RSl ! •- : 

OR WHENEVER FDATE(l) .L. 0 .CR. FCAIE(l) .G. 12 

TRANSFER TC . E R_tLG 
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Cf< WHENEVER FnAT£(2) .L. 0 .CR. F C AT E ( 2 ) .G. 31 

TRANSFER TC GRR7 

OR WHENEVER FUA TE l 3 ) .L. 0 .CR. Ft A T E ( 3 ) .G . 99 

TRANSFER TC FRR8 

END CF CONDITIONAL 

THRCLGH LP2, F OR KsM, K »G. 3 

WHENEVER UEGIMK) .L. NA 

TRANSFER TC ERR. 9 

CR WHENEVER ENDUU .G. NB 

TRANSFER TC„ ER 8.1.0 _ _ 

END CF CCNUITICNAL 

LP2 CONTINUE 

END CF CCNC I T1CNAL 

FLNC TI ON RE 1 CRN' 

STEIM 3 ) FUNCTION RETURN 

STEPUI LINES = LINES + 1 _ 1 

J = ( TO T -0 • , K = 1 , 1 , K . G • P T II ) , TC T + . L I . ( PNT -$32 . ( K , 1 ) ) / F l ) + . 5 

WhEUF. VrR FLAG_.,F._q 

FLAG = 1 


1 

f 


1 

1 

2 


CUT 

STEPl 6 ) 


1 


STLOtV 

i 


STEpU ) 


LP4 
L P 3 


TRANSFER TC CUT _ _ 

END CF CCiN’DI \ ICNAU 

WHENEVER PASS .6. 1 

<K=l,l,K.G.3, TCMK+1T * TCTIML1 ♦ 

(, TCT = p. T 1=1 ? 1 T 1 . Q.PL jKJLi 

TCT + .11. < FNT-S32. ( ltKIll) 

C I HER W I SF 

( K= 1 1 1 , K . G • 3 » TCT(K+1)“ = TCT ( K * 1 ) + 

i ICJeO . • 1 « I t_l 9 1 • G . PT < K ) f 

TCT + (.LI. ( F NT - 5 3 2 • ( l?K ) ) - 

.. » HF. (KEAN-S3Z. (FILFStK ))> - P.2 )) 

LNU OF CCNOI T I CNAL 

(K=1,L,K*G»3 f PT(KJ j_0 ) 

FUNCTION RETURN 

WHENEVER PASS .E . 2, TR ANS FE R TC STLCEV 

FILES = FILES + I 

WHENEVER FILES .G. FILMAX, TRANSFER TC ERR11 

CK^1»X»K VG 3 ’* (^FAN-S32. (FILES, K) > .SI. 

TC T ( K+ 1 > / ( { ENC ( K ) - BEGIN(K) f 1 M ( L INE S - l ) ) _) 

X ( F 1 LC S ) = t* 1 SSCN 

(CATE-FILES) .SI. FLAT L ( 2 ) + lCO*FOATt(lJ _+ 1G0C0*FCATE ( 3) 

(K=l»l ,K.G. 0 t ( C K P N t - S 8 2 V ( FILES Vk") j .‘SU 'CC> M NT(K)) 

PASS = 2 

function "re Turn 

( K = 1 , 1,K«G#3,TCT = TCT(K+l)/( | L l NE S- 1) * ( E N C ( K ) - B E G I N l K > + 1 > 

(IJEV-S32. ( FILE'S# K ) >” ^ S l\ SCRT . ( TC T ) ) 

PASS = 1 1 

FUNCTION RETURN 

( K = L , 1 f K . G .F ILES, ( TAG-K ) .SI. K) 

T h R C U G H L P 3 » FCR K=1,L, K .G. (FILES-1) 

THROUGH LPA, FCR Is(K+l),l, J .G . FILES 

WHENEVER .LI.(CATE-.Ll.fTAG-K)) .GE. . L I . ( C ATE- • L I ♦ ( T AG- I ) ) 

J = .LI. (TAG-K) 

(TAG-K) .SI. .LI. ( T AG- I } 

( TAG- I ) .SU_J 

EinO CF CONDITIONAL 

CONTINUE _ 

CONTINUE 


NSCALE = 1 

EXeCCTC PLCT1. I NSC ALE, 5,10,6,10) 

F 2 = 1 . __ 

WHENEVER X ( F I LE S J -6 L . • GE. 1., F2 = X(FILES) - 61. 
FI = X ( F 1 LES ) 

EXECUTE PLCT2. (T^AGE ,Fl ,F2i510. ,0.’ ) 

THROUGH LP5, FOR 1=1,1, 1 .0. 3 


168 



FORMERLY WILLOW RUN LABORATORIES, THE UNIVERSITY OF MICHIGAN 



(K«l,l»K.G.FlLGS f Ym = . L I F. < ME AK-S32 . tK , I) >> 

F X t C U T F P L C T 3 , • [ 0 C C (I) ,X (l) ,Y (1) »FIL£S ) __ 

LP 5 COM 1 riliE 

PRINT FCRNAI UH1** 

EXECUTE P'L C T ■ ( 1 0 , K SG 1 ) 

PRIM rCKf'AT $lHC,T3etH^ISSiCN.KU?BER?** 

EXECUTE PLCT2. (INAGI: ,FL,F2 ,6. > 

THRUL’GH LPfct FCR 1 = 1,1, I «0» 3 

<K = J , I ik7G.h : 'lLES,Y{ K) = .LlF. (CEV-S32. Ut l M ) 

EXECUTE: PLCT3. (BCD ( I) ,XU ) ._YU1 [.rJFILES]_ 

LP 6 CCMJMt 

PRIM FORMAT 1 1 H 1 * 5 

EXECUTE FLCT4. llSfKSG?)" 

PRIM FORMAT iIHC.T26,H < ^ISS [CN NUMB ER'** 

' " 'PRINT FORMAT F L H 1 ,117, 2 1 C CL D PC T 1 , T 4 l ) t T ICC/ IH *T16, 

1 2 (H 1 PLATE PLATE A M b l t M 1 , 1 4 0 ) , T 1 00 / 

2 H f DATE MISSCN 1 , 3 l H ' M E M\ M ♦ 3 ( H * D E V 1 , S 5 > / / * * 

THROUGH LP7, FOR 1=1,1, l FILES 

<K=l,l,K„G.3tY(K) = * L I T • (NEAK-S32. ( * L 1. ( T AG- I) , K >) , 

1 y ( K + 3 1 = . L l F * ( D£ V-S32. ( . L U (TA G - I ), K_ JJJ 

(K=l ,1 ,K.g7r7'cC^PNTTkT' = .Li . (CMMNT-S8 2. ( .L WITAG- M ,K )) ) 

K = .11 . (CATR-*LI - TTAG-I )> 

K = (K - K/lCCCC*iCCCC) *100 + K/1CCC0 

PRIM FORMAT flH , t 6 , SI ,F6 , SL , k ( F 7 .2 , SI )/ IH ,115,606**, 

1 K t X t.LI . (TAG-I M tYC 1)'. ••'¥(6) iCDFKK-r ( U .. .CCffMNT( 8) 


LP f 


l 


BLEhl 


ERR 1 


CCMINUE ___ 

PKUT FORMAT CCLC PLATE'/F' * - FCT PL A T E 1 / H • 8 - AMBIENT 1 / 

1 H 1 « * , 

F CM: T ICR KEiLPN 

PR 1 M FORM A T ERR ( I ) ,K 

E f< k L R * 

I = 1 

TRANSFER 'IC BLEH 


tR P2 I = 43 

TRANSFER TC 0 L EH 
fc RR 3 l = 85 

TRANSFER TC BLEH 
tRR A I = 127 

- ; A Si S F E R " T C "F i L ETC 


ERP.5 
t R X6 
ERR 7 
fc R P. 9 
FRR 9 
tRP IQ 
ERRli 
LL cH 


l = 169 

TRANSFER 1C ElLEH 

I = 211 

TRANSFER TC BLEH 

1 = 253 

TRANSFER TC PL EH 
1 * 295 

TRANSFER TC BLEH 
I = 337 _.... 

TRANSFER TC BLEH1 

1 = 379 

TRANSFER TC BLEH 1 

1 = 421 

TRANSFER TC BLEH 
PRIM FORMAT FRRIU 
ERROR* 


I NT EkfJA L FUNCTION .PTP RjDj C_._ 

THROUGH LP 8 , FCR K = l,l, K * G* 3 

WHENEVER IP *LE» ENC(K) .AND. IP *GE* REGIN(K) 
WHEnEVCR CJKCDE .E. I 
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K * FDATLF 

liATUM 1 ) = K 

END OF CPNOlTICrsAL _ 

IMtK) = P T t K I + 1 

UhENEvER K ♦ E . 1 _ 

(PKT-S32. IPT(l)tl) ) .51. CATUP(l) 

OTHERWISE 

(PM-S32. tPTCK) f K) ) .SI. DATUM 1 ) - J 

END CP CCM) I T I.PNA.L 

END CF CCN'CITICNAL 

IP 8 CONTINUE _ 

FLNCTLCN RETURN 

END OF FUNCTION _ _ 

END CF FUNCTION 
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